Search

Surfaces are able to control physical-chemical processes in multi-component solution systems and, as such, find application in a wide range of technological devices. Understanding the structure, dynamics and thermodynamics of non-ideal solutions at surfaces, however, is particularly challenging. Here, we use Constant Chemical Potential Molecular Dynamics simulations to gain insight into aqueous NaCl solutions in contact with graphite surfaces at high concentrations and under the effect of applied surface charges: conditions where mean-field theories describing interfaces cannot (typically) be reliably applied. We discover an asymmetric effect of surface charge on the electric double layer structure and resulting thermodynamic properties, which can be explained by considering the affinity of the surface for cations and anions and the cooperative adsorption of ions that occurs at higher concentrations. We characterise how the sign of the surface charge affects ion densities and water structure in the double layer and how the capacitance of the interface - a function of the electric potential drop across the double layer - is largely insensitive to the bulk solution concentration. Notably, we find that negatively charged graphite surfaces induce an increase in the size and concentration of extended liquid-like ion clusters confined to the double layer. Finally, we discuss how concentration and surface charge affect the activity coefficients of ions and water at the interface, demonstrating how electric fields in this region should be explicitly considered when characterising the thermodynamics of both solute and solvent at the solid/liquid interface.

Molecule- and particle-based simulations provide the tools to test, in microscopic detail, the validity of classical nucleation theory. In this endeavour, determining nucleation mechanisms and rates for phase separation requires an appropriately defined reaction coordinate to describe the transformation of an out-of-equilibrium parent phase, for which myriad options are available to the simulator. In this article, we describe the application of the variational approach to Markov processes (VAMP) to quantify the suitability of reaction coordinates to study crystallisation from supersaturated colloid suspensions. Our analysis indicates that collective variables (CVs) that correlate with the number of particles in the condensed phase, the system potential energy and approximate configurational entropy often feature as the most appropriate order parameters to quantitatively describe the crystallisation process. We apply time-lagged independent component analysis to reduce high-dimensional reaction coordinates constructed from these CVs to build Markov State Models (MSMs), which indicate that two barriers separate a supersaturated fluid phase from crystals in the simulated environment. The MSMs provide consistent estimates for crystal nucleation rates, regardless of the dimensionality of the order parameter space adopted; however, the two-step mechanism is only consistently evident from spectral clustering of the MSMs in higher dimensions. As the method is general and easily transferable, the variational approach we adopt could provide a useful framework to study controls for crystal nucleation.

We present the coupling of two frameworks -- the pseudo-open boundary simulation method known as constant potential Molecular Dynamics simulations (C$\mu$MD), combined with QMMD calculations -- to describe the properties of graphene electrodes in contact with electrolytes. The resulting C$\mu$QMMD model was then applied to three ionic solutions (LiCl, NaCl and KCl in water) at bulk solution concentrations ranging from 0.5 M up to 6 M in contact with a charged graphene electrode. The new approach we are describing here provides a simulation protocol to control the concentration of the electrolyte solutions while including the effects of a fully polarizable electrode surface. Thanks to this coupling, we are able to accurately model both the electrode and solution side of the double layer and provide a thorough analysis of the properties of electrolytes at charged interfaces, such as the screening ability of the electrolyte and the electrostatic potential profile. We also report the calculation of the integral electrochemical double layer capacitance in the whole range of concentrations analysed for each ionic species, while the QM simulations provide access to the differential and integral quantum capacitance. We highlight how subtle features, such as the adsorption of potassium at the interface or the tendency of the ions to form clusters, emerge from our simulations, contribute to explaining the ability of graphene to store charge and suggest implications for desalination., Comment: 28 pages, 10 figures

Twisted van der Waals multilayers are widely regarded as a rich platform to access novel electronic phases, thanks to the multiple degrees of freedom such as layer thickness and twist angle that allow control of their electronic and chemical properties. Here, we propose that the stacking domains that form naturally due to the relative twist between successive layers act as an additional "knob" for controlling the behavior of these systems, and report the emergence and engineering of stacking domain-dependent surface chemistry in twisted few-layer graphene. Using mid-infrared near-field optical microscopy and atomic force microscopy, we observe a selective adhesion of metallic nanoparticles and liquid water at the domains with rhombohedral stacking configurations of minimally twisted double bi- and tri-layer graphene. Furthermore, we demonstrate that the manipulation of nanoparticles located at certain stacking domains can locally reconfigure the moir\'e superlattice in their vicinity at the {\mu}m-scale. In addition, we report first-principles simulations of the energetics of adhesion of metal atoms and water molecules on the stacking domains in an attempt to elucidate the origin of the observed selective adhesion. Our findings establish a new approach to controlling moir\'e-assisted chemistry and nanoengineering., Comment: 11 pages, 3 figures

Moir\'e superlattices in twisted two-dimensional materials have generated tremendous excitement as a platform for achieving quantum properties on demand. However, the moir\'e pattern is highly sensitive to the interlayer atomic registry, and current assembly techniques suffer from imprecise control of the average twist angle, spatial inhomogeneity in the local twist angle, and distortions due to random strain. Here, we demonstrate a new way to manipulate the moir\'e patterns in hetero- and homo-bilayers through in-plane bending of monolayer ribbons, using the tip of an atomic force microscope. This technique achieves continuous variation of twist angles with improved twist-angle homogeneity and reduced random strain, resulting in moir\'e patterns with highly tunable wavelength and ultra-low disorder. Our results pave the way for detailed studies of ultra-low disorder moir\'e systems and the realization of precise strain-engineered devices.

NaCl crystal nucleation from metastable solutions has long been considered to occur according to a single-step mechanism. Recent experimental observations suggest that NaCl crystals emerge from disordered intermediates which is seemingly at odds with this established view. Here, we performed biased and unbiased molecular dynamics simulations to analyse and characterise the pathways to crystalline phases from solutions far into the metastable region. We find that large liquid-like NaCl clusters emerge as the solution concentration is increased and two-step nucleation pathways are more dominant than one-step pathways to phase separation. Analyses of cluster size populations and the ion pair association constant show that these clusters are transient, unlike the thermodynamically stable prenulceation cluster solute species that were suggested in other mineralising systems. A Markov State Model was developed to analyse the mechanisms and timescales for nucleation from unbiased molecular dynamics trajectories in a reaction coordinate space characterising the dense regions in clusters and crystalline order. This allowed calculation of the committor probabilities for the system to relax to the solution or crystal states, and to estimate the rate of nucleation, which shows excellent agreement with literature values. From a fundamental nucleation perspective, our work highlights the need to extend the attribute 'critical' to an ensemble of clusters which can display a broad range of structures and include sizeable disordered domains depending upon the reaction conditions.Moreover, our recent simulation studies demonstrated that carbon surfaces catalyse the formation of liquid-like NaCl networks which, combined with the observations here, suggests that alternative pathways beyond the single-step mechanism can be exploited to control the crystallisation of NaCl.

Solid/liquid interfaces control various processes of technological relevance in the process industry and many fundamental physicochemical phenomena. This work examines the link between the atomistic description of mass transfer at solid/liquid interfaces out of equilibrium and the constitutive mass transfer equations typically used to model these processes at the mesoscale. In our analysis, we discuss the microscopic inconsistencies apparent in simplified models of mass transfer whenever non-idealities dominate the liquid phase in the proximity of solid/liquid interfaces. Using C{\mu}MD - a molecular simulation technique to investigate out-of-equilibrium, concentration-driven processes with pseudo open-boundary conditions - we outline a strategy to capture and quantify non-idealities induced by specific interactions between the solid surface and molecules in the fluid phase. We demonstrate our approach by studying electrolyte solutions in technologically important, multi-component systems in which introducing a solid/liquid interface induces substantial deviations in the composition and electrochemical properties compared to the fluid phase bulk. To this aim, we analyse NaCl(aq) solutions in contact with i) the graphite basal surface and ii) {100} NaCl(s) crystalline surfaces. We uncover the tendency of the sodium cation to preferentially adsorb at the surfaces considered, which induces local violations of electroneutrality that leads to the emergence of an electric potential in the fluid phase at the solid/liquid interface., Comment: 12 pages, 4 figures

Graphite is a ubiquitous electrode material with particular promise for use in e.g., energy storage and desalination devices, but very little is known about the properties of the graphite-electrolyte double layer at technologically relevant concentrations. Here, the (electrified) graphite-NaCl(aq) interface was examined using constant chemical potential molecular dynamics simulations; this approach avoids ion depletion (due to surface adsorption) and maintains a constant concentration of ions beyond the surface. Specific Na+ adsorption at the graphite basal surface causes charging of the interface in the absence of an applied potential. At moderate bulk concentrations, this leads to accumulation of counter-ions in a diffuse layer to balance the effective surface charge, consistent with established models of the electrical double layer (DL). Beyond 0.6 M, however, a combination of over-screening and ion crowding in the DL results in alternating compact layers of ion density perpendicular to the interface. The transition to this regime is marked by an increasing DL size and anomalous negative shifts to the potential of zero charge with incremental changes to the bulk concentration. Our observations are supported by changes to the position of the differential capacitance minimum measured by electrochemical impedance spectroscopy. Furthermore, a striking level of agreement between the differential capacitance from simulations and experiments allows us to critically assess the accepted norm that electrochemical capacitance measurements report simply on the density of states of the graphite material. Finally, ion crowding at the highest concentrations (beyond 5 M) leads to the formation of liquid-like NaCl clusters confined to highly non-ideal regions of the double layer, where ion diffusion is up to five times slower than in the bulk.

Graphene-based heterostructures display a variety of phenomena that are strongly tunable by electrostatic local gates. Monolayer graphene (MLG) exhibits tunable surface plasmon polaritons, as revealed by scanning nano-infrared experiments. In bilayer graphene (BLG), an electronic gap is induced by a perpendicular displacement field. Gapped BLG is predicted to display unusual effects such as plasmon amplification and domain wall plasmons with significantly larger lifetime than MLG. Furthermore, a variety of correlated electronic phases highly sensitive to displacement fields have been observed in twisted graphene structures. However, applying perpendicular displacement fields in nano-infrared experiments has only recently become possible (Ref. 1). In this work, we fully characterize two approaches to realizing nano-optics compatible top-gates: bilayer $\text{MoS}_2$ and MLG. We perform nano-infrared imaging on both types of structures and evaluate their strengths and weaknesses. Our work paves the way for comprehensive near-field experiments of correlated phenomena and plasmonic effects in graphene-based heterostructures.

The emerging field of twistronics, which harnesses the twist angle between two-dimensional materials, represents a promising route for the design of quantum materials, as the twist-angle-induced superlattices offer means to control topology and strong correlations. At the small twist limit, and particularly under strain, as atomic relaxation prevails, the emergent moir\'e superlattice encodes elusive insights into the local interlayer interaction. Here we introduce moir\'e metrology as a combined experiment-theory framework to probe the stacking energy landscape of bilayer structures at the 0.1 meV/atom scale, outperforming the gold-standard of quantum chemistry. Through studying the shapes of moir\'e domains with numerous nano-imaging techniques, and correlating with multi-scale modelling, we assess and refine first-principle models for the interlayer interaction. We document the prowess of moir\'e metrology for three representative twisted systems: bilayer graphene, double bilayer graphene and H-stacked $MoSe_2/WSe_2$. Moir\'e metrology establishes sought after experimental benchmarks for interlayer interaction, thus enabling accurate modelling of twisted multilayers., Comment: Main-text: 15 pages, 3 figures; Changes: Fig. 1a revised, expanded reference list, added supplementary information

Broken symmetries induce strong nonlinear optical responses in materials and at interfaces. Twist angle can give complete control over the presence or lack of inversion symmetry at a crystal interface, and is thus an appealing knob for tuning nonlinear optical systems. In contrast to conventional nonlinear crystals with rigid lattices, the weak interlayer coupling in van der Waals (vdW) heterostructures allows for arbitrary selection of twist angle, making nanomechanical manipulation of fundamental interfacial symmetry possible within a single device. Here we report highly tunable second harmonic generation (SHG) from nanomechanically rotatable stacks of bulk hexagonal boron nitride (BN) crystals, and introduce the term twistoptics to describe studies of optical properties in dynamically twistable vdW systems. We observe SHG intensity modulated by a factor of more than 50, polarization patterns determined by moir\'e interface symmetry, and enhanced conversion efficiency for bulk crystals by stacking multiple pieces of BN joined by symmetry-broken interfaces. Our study provides a foundation for compact twistoptics architectures aimed at efficient, scalable, and tunable frequency-conversion, and demonstrates SHG as a robust probe of buried vdW interfaces., Comment: 23 Pages, 17 Figures

Moir\'e superlattices in van der Waals (vdW) heterostructures have given rise to a number of emergent electronic phenomena due to the interplay between atomic structure and electron correlations. A lack of a simple way to characterize moir\'e superlattices has impeded progress in the field. In this work we outline a simple, room-temperature, ambient method to visualize real-space moir\'e superlattices with sub-5 nm spatial resolution in a variety of twisted vdW heterostructures including but not limited to conducting graphene, insulating boron nitride and semiconducting transition metal dichalcogenides. Our method utilizes piezoresponse force microscopy, an atomic force microscope modality which locally measures electromechanical surface deformation. We find that all moir\'e superlattices, regardless of whether the constituent layers have inversion symmetry, exhibit a mechanical response to out-of-plane electric fields. This response is closely tied to flexoelectricity wherein electric polarization and electromechanical response is induced through strain gradients present within moir\'e superlattices. Moir\'e superlattices of 2D materials thus represent an interlinked network of polarized domain walls in a non-polar background matrix., Comment: 21 pages, 4 figures; updated text

The formation of interfacial moir\'e patterns from angular and/or lattice mismatch has become a powerful approach to engineer a range of quantum phenomena in van der Waals heterostructures. For long-lived and valley-polarized interlayer excitons in transition-metal dichalcogenide (TMDC) heterobilayers, signatures of quantum confinement by the moir\'e landscape have been reported in recent experimental studies. Such moir\'e confinement has offered the exciting possibility to tailor new excitonic systems, such as ordered arrays of zero-dimensional (0D) quantum emitters and their coupling into topological superlattices. A remarkable nature of the moir\'e potential is its dramatic response to strain, where a small uniaxial strain can tune the array of quantum-dot-like 0D traps into parallel stripes of one-dimensional (1D) quantum wires. Here, we present direct evidence for the 1D moir\'e potentials from real space imaging and the corresponding 1D moir\'e excitons from photoluminescence (PL) emission in MoSe2/WSe2 heterobilayers. Whereas the 0D moir\'e excitons display quantum emitter-like sharp PL peaks with circular polarization, the PL emission from 1D moir\'e excitons has linear polarization and two orders of magnitude higher intensity. The results presented here establish strain engineering as a powerful new method to tailor moir\'e potentials as well as their optical and electronic responses on demand.

In heterostructures consisting of atomically thin crystals layered on top of one another, lattice mismatch or rotation between the layers results in long-wavelength moir\'e superlattices. These moir\'e patterns can drive significant band structure reconstruction of the composite material, leading to a wide range of emergent phenomena including superconductivity, magnetism, fractional Chern insulating states, and moir\'e excitons. Here, we investigate monolayer graphene encapsulated between two crystals of boron nitride (BN), where the rotational alignment between all three components can be varied. We find that band gaps in the graphene arising from perfect rotational alignment with both BN layers can be modified substantially depending on whether the relative orientation of the two BN layers is 0 or 60 degrees, suggesting a tunable transition between the absence or presence of inversion symmetry in the heterostructure. Small deviations ($<1^{\circ}$) from perfect alignment of all three layers leads to coexisting long-wavelength moir\'e potentials, resulting in a highly reconstructed graphene band structure featuring multiple secondary Dirac points. Our results demonstrate that the interplay between multiple moir\'e patterns can be utilized to controllably modify the electronic properties of the composite heterostructure., Comment: 14 Pages, 12 Figures