How to subtract cylinder from a cube, add bevel and make it all smooth
I'd recommend starting with autosmooth along with a weighted normals modifier to help the result after
a bevel modifier.
Here is my cube's modifier stack:

Here is the result: