  Analysis of the strain-induced pseudomagnetic fields (PMFs) generated in
graphene nanobulges under three different substrate scenarios shows that, in
addition to the shape, the graphene-substrate interaction can crucially
determine the overall distribution and magnitude of strain and those fields, in
and outside the bulge. We utilize a combination of classical molecular
dynamics, continuum mechanics, and tight-binding electronic structure
calculations as an unbiased means of studying pressure-induced deformations and
the resulting PMF in graphene nanobubbles of various geometries. The interplay
among substrate aperture geometry, lattice orientation, internal gas pressure,
and substrate type is analyzed in view of strain-engineered graphene
nanostructures capable of confining and/or guiding electrons at low energies.
Except in highly anisotropic geometries, the magnitude of the PMF is generally
significant only near the boundaries of the aperture and rapidly decays towards
the center because under gas pressure at the scales considered here there is
considerable bending at the edges and the central region displays nearly
isotropic strain. When the deflection lead to sharp bends at the edges,
curvature and the tilting of the $p_z$ orbitals cannot be ignored and
contributes substantially to the total field. The strong and localized nature
of the PMF at the boundaries and its polarity-changing profile can be exploited
to trap electrons inside the bubble or of guiding them in channel-like
geometries defined by edges. However, we establish that slippage of graphene
against the substrate is an important factor in determining the degree of
concentration of PMFs in or around the bulge since it can lead to considerable
softening of the strain gradients there. The nature of the substrate emerges
thus as a decisive factor determining the effectiveness of nanoscale PMFs
tailoring in graphene.
