Pseudomagnetic fields in graphene nanobubbles of constrained geometry: A molecular dynamics study

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.