Sequential Bending and Twisting around C-C Single Bonds by Mechanical Lifting of a Pre-Adsorbed Polymer

Explicit Modelling of Spectral Bandshapes by a Mixed Quantum-Classical Approach: Solvent Order and Temperature Effects in the Optical Spectra of Distyrylbenzene

The absorption and emission spectral shapes of a flexible organic probe, the distyrylbenzene (DSB) dye, are simulated accounting for the effect of different environments of increasing complexity, ranging from a homogeneous, low-molecular- weight solvent, to a long-chain alkane, and, eventually, a channel-forming organic matrix. Each embedding is treated explicitly, adopting a mixed quantum-classical approach, the Adiabatic Molecular Dynamics - generalized vertical Hessian (Ad-MD|gVH) model, which allows a direct simulation of the environment-induced constraining effects on the vibronic spectral shapes. In such a theoretical framework, the stiff modes of the dye are described at a quantum level within the harmonic approximation, including Duschinsky mixing effects, while flexible degrees of freedom of the solute (e. g. torsions) and those of the solvent are treated classically by means of molecular dynamics sampling. Such a setup is shown to reproduce the distinct effects exerted by the different environments in varied thermodynamic conditions. Besides allowing for a first-principles rationale on the supramolecular mechanism leading to the experimental spectral features, this result represents the first successful application of the Ad-MD|gVH method to complex embeddings and supports its potential application to other heterogeneous environments, such as for instance, pigment-protein complexes or organic dyes adsorbed into metal-organic frameworks.

Dynamic effects on the nonlinear optical properties of donor acceptor stenhouse adducts: insights from combined MD + QM simulations

The second-order nonlinear optical (NLO) responses of a donor-acceptor stenhouse adduct (DASA) are investigated by using a computational approach combining molecular dynamics simulations and density functional theory (DFT) calculations. Specific force fields for the open and closed photoswitching forms are first parameterized and validated according to the Joyce protocol, in order to finely reproduce the geometrical features and potential energy surfaces of both isomers in chloroform solution. Then, DFT calculations are performed on structural snapshots extracted at regular time steps of the MD trajectories to address the influence of the thermalized conformational dynamics on the NLO responses related to hyper-Rayleigh scattering (HRS) experiments. We show that accounting for the structural dynamics largely enhances the HRS hyperpolarizability (βHRS) compared to DFT calculations considering solely equilibrium geometries, and greatly improves the agreement with experimental measurements. Furthermore, we show that the NLO responses of the NLO-active open form are correlated with the bond order alternation along the triene bridge connecting the donor and acceptor moieties, which is rationalized using simple essential state models.

Shear-activation of mechanochemical reactions through molecular deformation

Mechanical stress can directly activate chemical reactions by reducing the reaction energy barrier. A possible mechanism of such mechanochemical activation is structural deformation of the reactant species. However, the effect of deformation on the reaction energetics is unclear, especially, for shear stress-driven reactions. Here, we investigated shear stress-driven oligomerization reactions of cyclohexene on silica using a combination of reactive molecular dynamics simulations and ball-on-flat tribometer experiments. Both simulations and experiments captured an exponential increase in reaction yield with shear stress. Elemental analysis of ball-on-flat reaction products revealed the presence of oxygen in the polymers, a trend corroborated by the simulations, highlighting the critical role of surface oxygen atoms in oligomerization reactions. Structural analysis of the reacting molecules in simulations indicated the reactants were deformed just before a reaction occurred. Quantitative evidence of shear-induced deformation was established by comparing bond lengths in cyclohexene molecules in equilibrium and prior to reactions. Nudged elastic band calculations showed that the deformation had a small effect on the transition state energy but notably increased the reactant state energy, ultimately leading to a reduction in the energy barrier. Finally, a quantitative relationship was developed between molecular deformation and energy barrier reduction by mechanical stress.

Accurate Quantum-Mechanically Derived Force-Fields through a Fragment-Based Approach: Balancing Specificity and Transferability in the Prediction of Self-Assembly in Soft Matter

The wide range of time/length scales covered by self-assembly in soft matter makes molecular dynamics (MD) the ideal candidate for simulating such a supramolecular phenomenon at an atomistic level. However, the reliability of MD outcomes heavily relies on the accuracy of the adopted force-field (FF). The spontaneous re-ordering in liquid crystalline materials stands as a clear example of such collective self-assembling processes, driven by a subtle and delicate balance between supramolecular interactions and single-molecule flexibility. General-purpose transferable FFs often dramatically fail to reproduce such complex phenomena, for example, the error on the transition temperatures being larger than 100 K. Conversely, quantum-mechanically derived force-fields (QMD-FFs), specifically tailored for the target system, were recently shown (J. Phys. Chem. Lett.2022,13, 243) to allow for the required accuracy as they not only well reproduced transition temperatures but also yielded a quantitative agreement with the experiment on a wealth of structural, dynamic, and thermodynamic properties. The main drawback of this strategy stands in the computational burden connected to the numerous quantum mechanical (QM) calculations usually required for a target-specific parameterization, which has undoubtedly hampered the routine application of QMD-FFs. In this work, we propose a fragment-based strategy to extend the applicability of QMD-FFs, in which the amount of QM calculations is significantly reduced, being a single-molecule-optimized geometry and its Hessian matrix the only QM information required. To validate this route, a new FF is assembled for a large mesogen, exploiting the parameters obtained for two smaller liquid crystalline molecules, in this and previous work. Lengthy MD simulations are carried out with the new transferred QMD-FF, observing again a spontaneous re-orientation in the correct range of temperatures, with good agreement with the available experimental measures. The present results strongly suggest that a partial transfer of QMD-FF parameters can be invoked without a significant loss of accuracy, thus paving the way to exploit the method's intrinsic predictive capabilities in the simulation of novel soft materials.

Flexible Superlubricity Unveiled in Sidewinding Motion of Individual Polymeric Chains

A combination of low temperature atomic force microcopy and molecular dynamic simulations is used to demonstrate that soft designer molecules realize a sidewinding motion when dragged over a gold surface. Exploiting their longitudinal flexibility, pyrenylene chains are indeed able to lower diffusion energy barriers via on-surface directional locking and molecular strain. The resulting ultralow friction reaches values on the order of tens of pN reported so far only for rigid chains sliding on an incommensurate surface. Therefore, we demonstrate how molecular flexibility can be harnessed to realize complex nanomotion while retaining a superlubric character. This is in contrast with the paradigm guiding the design of most superlubric nanocontacts (mismatched rigid contacting surfaces).

Critical peeling of tethered nanoribbons

The peeling of an immobile adsorbed membrane is a well-known problem in engineering and macroscopic tribology. In the classic setup, picking up at one extreme and pulling off result in a peeling force that is a decreasing function of the pickup angle. As one end is lifted, the detachment front retracts to meet the immobile tail. At the nanoscale, interesting situations arise with the peeling of graphene nanoribbons (GNRs) on gold, as realized, e.g., by atomic force microscopy. The nanosized system shows a constant-force steady peeling regime, where the tip lifting h produces no retraction of the ribbon detachment point, and just an advancement h of the free tail end. This is opposite to the classic case, where the detachment point retracts and the tail end stands still. Here we characterise, by analytical modeling and numerical simulations, a third, experimentally relevant setup where the nanoribbon, albeit structurally lubric, does not have a freely moving tail end, which is instead elastically tethered. Surprisingly, novel nontrivial scaling exponents appear that regulate the peeling evolution. As the detachment front retracts and the tethered tail is stretched, power laws of h characterize the shrinking of the adhered length, the growth of peeling force and the peeling angle. These exponents precede the final total detachment as a critical point, where the entire ribbon eventually hangs suspended between the tip and the tethering spring. These analytical predictions are confirmed by realistic MD simulations, retaining the full atomistic description, also confirming their survival at finite experimental temperatures.

Topographic signatures and manipulations of Fe atoms, CO molecules and NaCl islands on superconducting Pb(111)

Topological superconductivity emerging in one- or two-dimensional hybrid materials is predicted as a key ingredient for quantum computing. However, not only the design of complex heterostructures is primordial for future applications but also the characterization of their electronic and structural properties at the atomic scale using the most advanced scanning probe microscopy techniques with functionalized tips. We report on the topographic signatures observed by scanning tunneling microscopy (STM) of carbon monoxide (CO) molecules, iron (Fe) atoms and sodium chloride (NaCl) islands deposited on superconducting Pb(111). For the CO adsorption a comparison with the Pb(110) substrate is demonstrated. We show a general propensity of these adsorbates to diffuse at low temperature under gentle scanning conditions. Our findings provide new insights into high-resolution probe microscopy imaging with terminated tips, decoupling atoms and molecules by NaCl islands or tip-induced lateral manipulation of iron atoms on top of the prototypical Pb(111) superconducting surface.

Practical Guide to Single-Protein AFM Nanomechanical Spectroscopy Mapping: Insights and Pitfalls As Unraveled by All-Atom MD Simulations on Immunoglobulin G

Atomic force microscopy is an invaluable characterization tool in almost every biophysics laboratory. However, obtaining atomic/sub-nanometer resolution on single proteins has thus far remained elusive-a feat long achieved on hard substrates. In this regard, nanomechanical spectroscopy mapping may provide a viable approach to overcome this limitation. By complementing topography with mechanical properties measured locally, one may thus enhance spatial resolution at the single-protein level. In this work, we perform all-atom molecular dynamics simulations of the indentation process on a single immunoglobulin G (IgG) adsorbed on a graphene slab. Our simulations reveal three different stages as a function of strain: a noncontact regime-where the mechanical response is linked to the presence of the water environment- followed by an elastic response and a final plastic deformation regime. In the noncontact regime, we are able to identify hydrophobic/hydrophilic patches over the protein. This regime provides the most local mechanical information that allows one to discern different regions with similar height/topography and leads to the best spatial resolution. In the elastic regime, we conclude that the Young modulus is a well-defined property only within mechanically decoupled domains. This is caused by the fact that the elastic deformation is associated with a global reorganization of the domain. Differences in the mechanical response are large enough to clearly resolve domains within a single protein, such as the three subunits forming the IgG. Two events, unfolding or protein slipping, are observed in the plastic regime. Our simulations allow us to characterize these two processes and to provide a strategy to identify them in the force curves. Finally, we elaborate on possible challenges that could hamper the interpretation of such experiments/simulations and how to overcome them. All in all, our simulations provide a detailed picture of nanomechanical spectroscopy mapping on single proteins, showing its potential and the challenges that need to be overcome to unlock its full potential.