Scaling laws for van der Waals interactions in nanostructured materials

Four-Dimensional Scaling of Dipole Polarizability: From Single-Particle Models to Atoms and Molecules

Scaling laws enable the determination of physicochemical properties of molecules and materials as a function of their size, density, number of electrons or other easily accessible descriptors. Such relations can be counterintuitive and nonlinear, and ultimately yield much needed insight into quantum mechanics of many-particle systems. In this work, we show on the basis of single-particle models, multielectron atoms and molecules that the dipole polarizability of quantum systems is generally proportional to the fourth power of a characteristic length, computed from the ground-state wave function. This four-dimensional (4D) scaling is independent of the ratio of bound-to-bound and bound-to-continuum electronic transitions and applies to many-electron atoms when a correlated length metric is used. Finally, this scaling law is applied to predict the polarizability of molecules by electrostatically coupled atoms-in-molecules approach, obtaining approximately 8% absolute and relative accuracy with respect to hybrid density functional theory (DFT) on the QM7-X data set of organic molecules, providing an efficient and scalable model for the anisotropic polarizability tensors of extended (bio)molecules.

System-Specific Parameter Optimization for Nonpolarizable and Polarizable Force Fields

The accuracy of classical force fields (FFs) has been shown to be limited for the simulation of cation-protein systems despite their importance in understanding the processes of life. Improvements can result from optimizing the parameters of classical FFs or by extending the FF formulation by terms describing charge transfer (CT) and polarization (POL) effects. In this work, we introduce our implementation of the CTPOL model in OpenMM, which extends the classical additive FF formula by adding CT and POL. Furthermore, we present an open-source parametrization tool, called FFAFFURR, that enables the (system-specific) parametrization of OPLS-AA and CTPOL models. The performance of our workflow was evaluated by its ability to reproduce quantum chemistry energies and by molecular dynamics simulations of a zinc-finger protein.

Mixtures Recomposition by Neural Nets: A Multidisciplinary Overview

Artificial Neural Networks (ANNs) are transforming how we understand chemical mixtures, providing an expressive view of the chemical space and multiscale processes. Their hybridization with physical knowledge can bridge the gap between predictivity and understanding of the underlying processes. This overview explores recent progress in ANNs, particularly their potential in the 'recomposition' of chemical mixtures. Graph-based representations reveal patterns among mixture components, and deep learning models excel in capturing complexity and symmetries when compared to traditional Quantitative Structure-Property Relationship models. Key components, such as Hamiltonian networks and convolution operations, play a central role in representing multiscale mixtures. The integration of ANNs with Chemical Reaction Networks and Physics-Informed Neural Networks for inverse chemical kinetic problems is also examined. The combination of sensors with ANNs shows promise in optical and biomimetic applications. A common ground is identified in the context of statistical physics, where ANN-based methods iteratively adapt their models by blending their initial states with training data. The concept of mixture recomposition unveils a reciprocal inspiration between ANNs and reactive mixtures, highlighting learning behaviors influenced by the training environment.

Fluidization and Application of Carbon Nano Agglomerations

Carbon nanomaterials are unique with excellent functionality and diverse structures. However, agglomerated structures are commonly formed because of small-size effects and surface effects. Their hierarchical assembly into micro particles enables carbon nanomaterials to break the boundaries of classical Geldart particle classification before stable fluidization under gas-solid interactions. Currently, there are few systematic reports regarding the structural evolution and fluidization mechanism of carbon nano agglomerations. Based on existing research on carbon nanomaterials, this article reviews the fluidized structure control and fluidization principles of prototypical carbon nanotubes (CNTs) as well as their nanocomposites. The controlled agglomerate fluidization technology leads to the successful mass production of agglomerated and aligned CNTs. In addition, the self-similar agglomeration of individual ultralong CNTs and nanocomposites with silicon as model systems further exemplify the important role of surface structure and particle-fluid interactions. These emerging nano agglomerations have endowed classical fluidization technology with more innovations in advanced applications like energy storage, biomedical, and electronics. This review aims to provide insights into the connections between fluidization and carbon nanomaterials by highlighting their hierarchical structural evolution and the principle of agglomerated fluidization, expecting to showcase the vitality and connotation of fluidization science and technology in the new era.

Self-Assembly via Condensation of Imine or Its N-Substituted Derivatives

ConspectusCompared to traditionally used irreversible chemical reactions, dynamic covalent chemistry (DCC) including imine formation represents a more advanced technique in the preparation of molecules with complex structures and topologies, whose syntheses require the formation of many bonds. By allowing the occurrence of error checking and self-correcting, it is likely that the target molecules with high enough thermodynamic stability could be self-assembled in high or even quantitative yield. Two questions are raised herein. First, it becomes a central problem in self-assembly that how to endow a target product with high enough thermodynamic stability so that it can be produced as the major or the only product within the self-assembly library. Second, the reversible nature of dynamic bonds jeopardizes the intrinsic stability of the products. More specifically, the imine bond which represents the mostly used dynamic covalent bond, is apt to undergo hydrolysis in the presence of water. Developing new approaches to make imine more robust and compatible with water is thus of importance. In this account, we summarized the progress made in our group in the field of self-assembly based on C═N bond formation. In organic solvent where an imine bond is relatively robust, we focus on studying how to enhance the thermodynamic stability of a target molecule by introducing intramolecular forces. These noncovalent interactions either release enthalpy to favor the formation of the target molecule or preorganize the building blocks into specific conformations that mimic the product, so that the entropy loss of the formation of the latter is thus suppressed. In water, which often leads to imine hydrolysis, we developed two strategies to enhance the water-compatibility. By taking advantage of multivalency, namely, multiple bonds are often more robust than a single bond, self-assembly via condensation of imine was performed successfully in water, a solvent that is considered as forbidden zone of imine. Another approach is to replace typical imine with its more robust and water compatible derivatives, namely, either hydrazone or oxime, whose C═N bonds are generally less electrophilic compared to typical imine. With the water-compatible dynamic bonds in hand, a variety topological nontrivial molecules such as catenanes and knots was self-assembled successfully in aqueous media, driven by hydrophobic effect. When the self-assembled molecules in the form of rings and cages were designed for supramolecular purposes, water-compatibility endows a merit that allows the hosts to take advantage of hydrophobic effect to drive host-guest recognition, enabling various tasks to be accomplished, such as separation of guest isomers with similar physical properties, recognition of highly hydrated anions, as well as stabilization of guest dimers.

Enhanced capacitive pressure sensing performance by charge generation from filler movement in thin and flexible PVDF-GNP composite films

This study introduces an approach to overcome the limitations of conventional pressure sensors by developing a thin and lightweight composite film specifically tailored for flexible capacitive pressure sensors, with a particular emphasis on the medium and high pressure range. To accomplish this, we have engineered a composite film by combining polyvinylidene fluoride (PVDF) and graphite nanoplatelets (GNP) derived from expanded graphite (Ex-G). A uniform sized GNPs with an average lateral size of 2.55av and an average thickness of 33.74 av with narrow size distribution was obtained with a gas-induced expansion of expandable graphite (EXP-G) combined with tip sonication in solvent. By this precisely controlled GNP within the composite film, a remarkable improvement in sensor sensitivity has been achieved, surpassing 4.18 MPa-1 within the pressure range of 0.1 to 1.6 MPa. This enhancement can be attributed to the generation of electric charge from the movement of GNP in the polymer matrix. Additionally, stability testing has demonstrated the reliable operation of the composite film over 1000 cycles. Notably, the composite film exhibits exceptional continuous pressure sensing capabilities with a rapid response time of approximately 100 milliseconds. Experimental validation using a 3 × 3 sensor array has confirmed the accurate detection of specific contact points, thus highlighting the potential of the composite film in selective pressure sensing. These findings signify an advancement in the field of flexible capacitive pressure sensors that offer enhanced sensitivity, consistent operation, rapid response time, and the unique ability to selectively sense pressure.

Simultaneous quantification of Young's modulus and dispersion forces with nanoscale spatial resolution

Many advances in polymers and layered materials rely on a precise understanding of the local interactions between adjacent molecular or atomic layers. Quantifying dispersion forces at the nanoscale is particularly challenging with existing methods often time consuming, destructive, relying on surface averaging or requiring bespoke equipment. Here, we present a non-invasive method able to quantify the local mechanical and dispersion properties of a given sample with nanometer lateral precision. The method, based on atomic force microscopy (AFM), uses the frequency shift of a vibrating AFM cantilever in combination with established contact mechanics models to simultaneously derive the Hamaker constant and the effective Young's modulus at a given sample location. The derived Hamaker constant and Young's modulus represent an average over a small (typically <100) number of molecules or atoms. The oscillation amplitude of the vibrating AFM probe is used to select the length-scale of the features to analyse, with small vibrations able to resolve the contribution of sub-nanometric defects and large ones exploring effectively homogeneous areas. The accuracy of the method is validated on a range of 2D materials in air and water as well as on polymer thin films. We also provide the first experimental measurements of the Hamaker constant of HBN, MoT2, WSe2and polymer films, verifying theoretical predictions and computer simulations. The simplicity and robustness of the method, implemented with a commercial AFM, may support a broad range of technological applications in the growing field of polymers and nanostructured materials where a fine control of the van der Waals interactions is crucial to tune their properties.

Stereodynamic Strategies to Induce and Enrich Chirality of Atropisomers at a Late Stage

Enantiomers, where chirality arises from restricted rotation around a single bond, are atropisomers. Due to the unique nature of the origins of their chirality, synthetic strategies to access these compounds in an enantioselective manner differ from those used to prepare enantioenriched compounds containing point chirality arising from an unsymmetrically substituted carbon center. In particular stereodynamic transformations, such as dynamic kinetic resolutions, thermodynamic dynamic resolutions, and deracemizations, which rely on the ability to racemize or interconvert enantiomers, are a promising set of transformations to prepare optically pure compounds in the late stage of a synthetic sequence. Translation of these synthetic approaches from compounds with point chirality to atropisomers requires an expanded toolbox for epimerization/racemization and provides an opportunity to develop a new conceptual framework for the enantioselective synthesis of these compounds.

Electrochemically Versatile Graphite Nanoplatelets Prepared by a Straightforward, Highly Efficient, and Scalable Route

Nanostructured carbon materials with tailor-made structures (e.g., morphology, topological defect, dopant, and surface area) are of significant interest for a variety of applications. However, the preparation method selected for obtaining these tailor-made structures determines the area of application, precluding their use in other technological areas of interest. Currently, there is a lack of simple and low-cost methodologies versatile enough for obtaining freestanding carbon nanostructures that can be used in either energy storage or chemical detection. Here, a novel methodology for the development of a versatile electrochemically active platform based on freestanding graphite nanoplatelets (GNP) has been developed by exploiting the interiors of hollow carbon nanofibers (CNF) comprising nanographene stacks using dry ball-milling. Even though ball-milling could be considered as a universal method for any carbonaceous material, often, it is not as simple (one step, no purification, and no solvents), efficient (just GNP without tubular structures), and quick (just 20 min) as the sustainable method developed in this work, free of surfactants and stabilizer agents. We demonstrate that the freestanding GNP developed in this work (with an average thickness of 3.2 nm), due to the selective edge functionalization with the minimal disruption of the basal plane, can act either as a supercapacitor or as a chemical sensor, showing both a dramatic improvement in the charge storage ability of more than 30 times and an enhanced detection of electrochemically active molecules such as ascorbic acid with a 236 mV potential shift with respect to CNF in both cases. As shown here, GNP stand as an excellent versatile alternative compared to the standard commercially available carbon-based materials. Overall, our approach paves the way for the discovery of new nanocarbon-based electrochemical active platforms with a wide electrochemical applicability.

The Halogen Bond in Weakly Bonded Complexes and the Consequences for Aromaticity and Spin-Orbit Coupling

The halogen bond complexes CF3X⋯Y and C2F3X⋯Y, with Y = furan, thiophene, selenophene and X = Cl, Br, I, have been studied by using DFT and CCSD(T) in order to understand which factors govern the interaction between the halogen atom X and the aromatic ring. We found that PBE0-dDsC/QZ4P gives an adequate description of the interaction energies in these complexes, compared to CCSD(T) and experimental results. The interaction between the halogen atom X and the π-bonds in perpendicular orientation is stronger than the interaction with the in-plane lone pairs of the heteroatom of the aromatic cycle. The strength of the interaction follows the trend Cl < Br < I; the chalcogenide in the aromatic ring nor the hybridization of the C−X bond play a decisive role. The energy decomposition analysis shows that the interaction energy is dominated by all three contributions, viz., the electrostatic, orbital, and dispersion interactions: not one factor dominates the interaction energy. The aromaticity of the ring is undisturbed upon halogen bond formation: the π-ring current remains equally strong and diatropic in the complex as it is for the free aromatic ring. However, the spin-orbit coupling between the singlet and triplet π→π* states is increased upon halogen bond formation and a faster intersystem crossing between these states is therefore expected.

Prediction of the electronic structure of single-walled GeS nanotubes

The structure and electronic properties of puckered GeS nanotubes have been investigated using first-principles density functional theory calculation. Our results show that both the armchair and zigzag GeS nanotubes are semiconductor materials with an adjustable band gap. The band gap increases gradually with increasing the tube diameter, and slowly converges to the monolayer limit. On the application of strain, the GeS nanotubes provide interesting strain-induced band gap variation. When the compressive strain reached 20%, zigzag GeS nanotubes are completely transformed into armchair GeS nanotubes. In addition, the elastic properties of the relatively stable armchair GeS nanotubes have been studied, the Young's modulus of the armchair (11, 11), (13, 13) and (15, 15) nanotubes were calculated to be 227.488 GPa, 211.888 GPa and 213.920 GPa, respectively. Our work confirms that compared with carbon nanotubes, two-dimensional materials with a puckered structure are easier to realize phase transition by stress.

Small Size, Big Impact: Recent Progress in Bottom-Up Synthesized Nanographenes for Optoelectronic and Energy Applications

Bottom-up synthesized graphene nanostructures, including 0D graphene quantum dots and 1D graphene nanoribbons, have recently emerged as promising candidates for efficient, green optoelectronic, and energy storage applications. The versatility in their molecular structures offers a large and novel library of nanographenes with excellent and adjustable optical, electronic, and catalytic properties. In this minireview, recent progress on the fundamental understanding of the properties of different graphene nanostructures, and their state-of-the-art applications in optoelectronics and energy storage are summarized. The properties of pristine nanographenes, including high emissivity and intriguing blinking effect in graphene quantum dots, superior charge transport properties in graphene nanoribbons, and edge-specific electrochemistry in various graphene nanostructures, are highlighted. Furthermore, it is shown that emerging nanographene-2D material-based van der Waals heterostructures provide an exciting opportunity for efficient green optoelectronics with tunable characteristics. Finally, challenges and opportunities of the field are highlighted by offering guidelines for future combined efforts in the synthesis, assembly, spectroscopic, and electrical studies as well as (nano)fabrication to boost the progress toward advanced device applications.

Shape-Controlled Self-Assembly of Light-Powered Microrobots into Ordered Microchains for Cells Transport and Water Remediation

Nature presents the collective behavior of living organisms aiming to accomplish complex tasks, inspiring the development of cooperative micro/nanorobots. Herein, the spontaneous assembly of hematite-based microrobots with different shapes is presented. Autonomous motile light-driven hematite/Pt microrobots with cubic and walnut-like shapes are prepared by hydrothermal synthesis, followed by the deposition of a Pt layer to design Janus structures. Both microrobots show a fuel-free motion ability under light irradiation. Because of the asymmetric orientation of the magnetic dipole moment in the crystal, cubic hematite/Pt microrobots can self-assemble into ordered microchains, contrary to the random aggregation observed for walnut-like microrobots. The microchains exhibit different synchronized motions under light irradiation depending on the mutual orientation of the individual microrobots during the assembly, which allows them to accomplish multiple tasks, including capturing, picking up, and transporting microscale objects, such as yeast cells and suspended matter in water extracted from personal care products, as well as degrading polymeric materials. Such light-powered self-assembled microchains demonstrate an innovative cooperative behavior for small-scale multitasking artificial robotic systems, holding great potential toward cargo capture, transport, and delivery, and wastewater remediation.

Covalent RGD-graphene-phthalocyanine nanocomposite for fluorescence imaging-guided dual active/passive tumor-targeted combinatorial phototherapy

Poor tumor selectivity, low stability and quenched fluorescence are the main challenges to be overcome for nanomedicine, and are mainly caused by the dissociation of the nanostructure and aggregation of chromophores in the biological environment. Herein, covalently connected nanoparticles RGD-graphene-phthalocyanine (RGD-GO-SiPc) were constructed based on RGD peptide, silicon phthalocyanine (SiPc) and graphene oxide (GO) via a conjugation reaction for fluorescence imaging-guided cancer-targeted combinatorial phototherapy. The prepared RGD-GO-SiPc exhibited supreme biological stability, high-contrast fluorescence imaging, significantly enhanced NIR absorption, high photothermal conversion efficiency (25.6%), greatly improved cancer-targeting capability, and synergistic photodynamic (PDT) and photothermal therapy (PTT) efficacy along with low toxicity. Both in vitro and in vivo biological studies showed that RGD-GO-SiPc is a kind of promising multifunctional nanomedicine for fluorescence imaging-guided combined photothermal and photodynamic therapy with dual active/passive tumor-targeting properties.

Engineering ferromagnetism in Ni(OH)2 nanosheets using tunable uniaxial pressure in graphene oxide/reduced graphene oxide

The interlayer spaces in two dimensional (2D) layered materials such as graphene, metal oxides and metal chalcogenides can be used in a number of roles that include the trapping of gases, for ion transfer and for water purification applications. In such spaces, "inner" pressure occurs on guest species enclosed between the layers and its variation can, in principal, be used for precisely controlling particular guest properties. In this study, a mixture of two 2D materials including graphene oxide (GO) and nickel hydroxide (Ni(OH)₂), was employed to yield an anisotropic GO-Ni(OH)₂ hybrid 2D sheet. The inner pressure associated with this material was able to be tuned by reduction of the GO (to yield rGO) and this in turn was shown to affect the magnetic behaviour of Ni(OH)₂. The ferromagnetic transition temperature (T_c) for Ni(OH)₂ decreases as the interlayer distance became shorter, which is opposite to the behaviour observed for the application of hydrostatic pressure to the hybrid sheet. The uniaxial pressure affecting the interlayer of the 2D material, and generated by the reduction of GO to rGO, has the potential to not only influence the behaviour of a range of magnetic materials, but also individual properties of other types of functional materials.

Systematic optimization of a fragment-based force field against experimental pure-liquid properties considering large compound families: application to oxygen and nitrogen compounds

The CombiFF approach is a workflow for the automated refinement of force-field parameters against experimental condensed-phase data, considering entire classes of organic molecules constructed using a fragment library via combinatorial isomer enumeration. One peculiarity of this approach is that it relies on an electronegativity-equalization scheme to account for induction effects within molecules, with values of the atomic hardness and electronegativity as electrostatic parameters, rather than the partial charges themselves. In a previous article [M. P. Oliveira, M. Andrey, S. R. Rieder, L. Kern, D. F. Hahn, S. Riniker, B. A. C. Horta and P. H. Hünenberger, J. Chem. Theory. Comput. 2020, 16, 7525], CombiFF was introduced and applied to calibrate a GROMOS-compatible united-atom force field for the saturated acyclic (halo-)alkane family. Here, this scheme is employed for the construction of a corresponding force field for saturated acyclic compounds encompassing eight common chemical functional groups involving oxygen and/or nitrogen atoms, namely: ether, aldehyde, ketone, ester, alcohol, carboxylic acid, amine, and amide. Monofunctional as well as homo-polyfunctional compounds are considered. A total of 1712 experimental liquid densities ρliq and vaporization enthalpies ΔHvap concerning 1175 molecules are used for the calibration (339 molecules) and validation (836 molecules) of the 102 non-bonded interaction parameters of the force field. Using initial parameter values based on the GROMOS 2016H66 parameter set, convergence is reached after five iterations. Given access to one processor per simulated system, this operation only requires a few days of wall-clock computing time. After optimization, the root-mean-square deviations from experiment are 29.9 (22.4) kg m-3 for ρliq and 4.1 (5.5) kJ mol-1 for ΔHvap for the calibration (validation) set. Thus, a very good level of agreement with experiment is achieved in terms of these two properties, although the errors are inhomogeneously distributed across the different chemical functional groups.

Templated Assembly of Nanoparticles into Continuous Arrays

The templated assembly of nanoparticles has been limited so far to yield only discontinuous nanoparticle clusters confined within lithographically patterned cavities. Here, we explored the templated assembly of nanoparticles into continuous 2D structures, using lithographically patterned templates with topographical features sized as the assembled nanoparticles. We found that these features act as nucleation centers, whose exact arrangement determines four possible assembly regimes (i) rotated, (ii) disordered, (iii) closely packed, and (iv) unpacked. These regimes produce structures strikingly different from their geometry, orientation, long-range and short-range orders, and packing density. Interestingly, for templates with relatively distant nucleation centers, these four regimes are replaced with three new ones, which produce large monocrystalline domains that are either (i) uniformly rotated, (ii) uniformly aligned, or (iii) nonuniformly rotated relative to the nucleation lattice. We rationalized our experimental data using a mathematical model, which examines all the alignment possibilities between the nucleation centers and the ideal hexagonal assembly. Our finding provides a new approach for the à la carte obtainment of various nanoscale structures unachievable by natural self-assembly and opens a route for the fabrication of numerous functional nanodevices and nanosystems that could not be realized so far by the standard bottom-up approach.

On the forbidden graphene's ZO (out-of-plane optic) phononic band-analog vibrational modes in fullerenes

The study of nanostructures' vibrational properties is at the core of nanoscience research. They are known to represent a fingerprint of the system as well as to hint the underlying nature of chemical bonds. In this work, we focus on addressing how the vibrational density of states (VDOS) of the carbon fullerene family (C_n: n = 20 → 720 atoms) evolves from the molecular to the bulk material (graphene) behavior using density functional theory. We find that the fullerene's VDOS smoothly converges to the graphene characteristic line-shape, with the only noticeable discrepancy in the frequency range of the out-of-plane optic (ZO) phonon band. From a comparison of both systems we obtain as main results that: (1) The pentagonal faces in the fullerenes impede the existence of the analog of the high frequency graphene's ZO phonons, (2) which in the context of phonons could be interpreted as a compression (by 43%) of the ZO phonon band by decreasing its maximum allowed radial-optic vibration frequency. And 3) as a result, the deviation of fullerene's VDOS relative to graphene may hold important thermodynamical implications, such as larger heat capacities compared to graphene at room-temperature. These results provide insights that can be extrapolated to other nanostructures containing pentagonal rings or pentagonal defects.

Enhancing the CO2-to-CO Conversion from 2D Silver Nanoprisms via Superstructure Assembly

The electrochemical reduction of CO₂ in a highly selective and efficient manner is a crucial step toward its reuse for the production of chemicals and fuels. Nanostructured Ag catalysts have been found to be effective candidates for the conversion of CO₂-to-CO. However, the ambiguous determination of the intrinsic CO₂ activity and the maximization of the density of exposed active sites have greatly limited the use of Ag toward the realization of practical electrocatalytic devices. Here, we report a superstructure design strategy prepared by the self-assembly of two-dimensional Ag nanoprisms for maximizing the exposure of active edge ribs. The vertically stacked Ag nanoprisms allow the exposure of >95% of the edge sites, resulting in an enhanced selectivity and activity toward the production of CO from CO₂ with an overpotential of 152 mV. The Ag superstructures also demonstrate a selectivity of over 90% for 100 h together with a current retention of ≈94% at -600 mV versus the reversible hydrogen electrode and a partial energy efficiency for CO production of 70.5%. Our electrochemical measurements on individual Ag nanoprisms with various edge-to-basal plane ratios and the Ag superstructures led to the identification of the edge ribs as the active sites thanks to the ≈400 mV decrease in the onset potential compared to that of the Ag (111) basal planes and a turnover frequency of 9.2 × 10^-3 ± 1.9 × 10^-3 s^-1 at 0 V overpotential.

Electron confinement meet electron delocalization: non-additivity and finite-size effects in the polarizabilities and dispersion coefficients of the fullerenes

In this work, we used finite-field derivative techniques and density functional theory (DFT) to compute the static isotropic polarizability series (αl with l = 1, 2, 3) for the C60-C84 fullerenes and quantitatively assess the intrinsic non-additivity in these fundamental response properties. By comparing against classical models of the fullerenes as conducting spherical shells (or solid spheres) of uniform electron density, a detailed critical analysis of the derived effective scaling laws (α1 ∼ N1.2, α2 ∼ N2.0, α3 ∼ N2.7) demonstrates that the electronic structure of finite-sized fullerenes-a unique dichotomy of electron confinement and delocalization effects due to their quasi-spherical cage-like structures and encapsulated void spaces-simultaneously limits and enhances their quantum mechanical response to electric field perturbations. Corresponding frequency-dependent polarizabilities were obtained by inputting the αl series into the hollow sphere model (within the modified single frequency approximation), and used to compute the molecular dispersion coefficients (Cn with n = 6, 8, 9, 10) needed to describe the non-trivial van der Waals (vdW) interactions in fullerene-based systems. Using first-order perturbation theory in conjunction with >140 000 DFT calculations, we also computed the non-negligible zero-point vibrational contributions to α1 in C60 and C70, thereby enabling a more accurate and direct comparison between theory and experiment for these quintessential nanostructures.

Size-dependent polarizabilities and van der Waals dispersion coefficients of fullerenes from large-scale complex polarization propagator calculations

While the anomalous non-additive size-dependencies of static dipole polarizabilities and van der Waals C₆ dispersion coefficients of carbon fullerenes are well established, the widespread reported scalings for the latter (ranging from N^2.2 to N^2.8) call for a comprehensive first-principles investigation. With a highly efficient implementation of the linear complex polarization propagator, we have performed Hartree-Fock and Kohn-Sham density functional theory calculations of the frequency-dependent polarizabilities for fullerenes consisting of up to 540 carbon atoms. Our results for the static polarizabilities and C₆ coefficients show scalings of N^1.2 and N^2.2, respectively, thereby deviating significantly from the previously reported values obtained with the use of semi-classical/empirical methods. Arguably, our reported values are the most accurate to date as they represent the first ab initio or first-principles treatment of fullerenes up to a convincing system size.

Predictive QM/MM Modeling of Modulations in Protein-Protein Binding by Lysine Methylation

Lysine methylation is a key regulator of protein-protein binding. The amine group of lysine can accept up to three methyl groups, and experiments show that protein-protein binding free energies are sensitive to the extent of methylation. These sensitivities have been rationalized in terms of chemical and structural features present in the binding pockets of methyllysine binding domains. However, understanding their specific roles requires an energetic analysis. Here we propose a theoretical framework to combine quantum and molecular mechanics methods, and compute the effect of methylation on protein-protein binding free energies. The advantages of this approach are that it derives contributions from all local non-trivial effects of methylation on induction, polarizability and dispersion directly from self-consistent electron densities, and at the same time determines contributions from well-characterized hydration effects using a computationally efficient classical mean field method. Limitations of the approach are discussed, and we note that predicted free energies of fourteen out of the sixteen cases agree with experiment. Critical assessment of these cases leads to the following overarching principles that drive methylation-state recognition by protein domains. Methylation typically reduces the pairwise interaction between proteins. This biases binding toward lower methylated states. Simultaneously, however, methylation also makes it easier to partially dehydrate proteins and place them in protein-protein complexes. This latter effect biases binding in favor of higher methylated states. The overall effect of methylation on protein-protein binding depends ultimately on the balance between these two effects, which is observed to be tuned via several combinations of local features.

Catalysis on Pristine 2D Materials via Dispersion and Electrostatic Interactions

Shape complementarity between catalyst and transition state structure is one of the cornerstones of chemical catalysis. Likewise, noncovalent interactions play a major role in catalysis. It has been predicted computationally and recently confirmed experimentally [Kroeger, A. A.; Hooper, J. F.; Karton, A. ChemPhysChem, 2020, 21, 1675-1681] that pristine graphene can efficiently catalyze chemical processes via π-interactions and shape complementarity. Here we show that other two-dimensional materials with different electronic structures and chemical compositions (h-BN and graphane) can also catalyze chemical processes that proceed via planar transition state structures. These include the bowl-to-bowl inversions in corannulene and sumanene and the rotation about the C-C bond in substituted biphenyls. This catalytic activity is achieved through shape complementarity between planar nanomaterial and planar transition state structure, enabling disproportionate stabilization of the transition state structures over the nonplanar reactants and products. A DFT-based energy decomposition analysis shows that this catalytic activity is mainly driven by dispersion and electrostatic forces, which together outweigh the Pauli repulsion term. These findings enrich and expand the concept of catalysis by pristine 2D materials.

Pristine Graphene as a Racemization Catalyst for Axially Chiral BINOL

Despite versatile applications of functionalized graphene in catalysis, applications of pure, unfunctionalized graphene in catalysis are in their infancy. This work uses both computational and experimental approaches to show that single-layer graphene can efficiently catalyze the racemization of axially chiral BINOL in solution. Using double-hybrid density functional theory (DHDFT) we calculate the uncatalyzed and catalyzed Gibbs free reaction barrier heights in a number of representative solvents of varying polarity: benzene, diphenyl ether, dimethylformamide (DMF), and water. These calculations show that (i) graphene can achieve significant catalytic efficiencies (▵▵G^≠ _cat ) varying between 47.2 (in diphenyl ether) and 60.7 (in DMF) kJ mol^-1 . An energy decomposition analysis reveals that this catalytic activity is driven by electrostatic and dispersion interactions. Based on these computational results, we explore the graphene-catalyzed racemization of axially chiral BINOL experimentally and show that single-layer graphene can efficiently catalyze this process. Whilst the uncatalyzed racemization requires high temperatures of over 200 °C, a pristine single-layer graphene catalyst makes it accessible at 60 °C.

Screening nature of the van der Waals density functional method: a review and analysis of the many-body physics foundation

We review the screening nature and many-body physics foundation of the van der Waals density functional (vdW-DF) method [Berland Ket al2015Rep. Prog. Phys.78066501], a systematic approach to construct truly nonlocal exchange-correlation energy density functionals. To that end we define and focus on a class of consistent vdW-DF versions that adhere to the Lindhard screening logic of the full method formulation. The consistent-exchange vdW-DF-cx version [Berland K and Hyldgaard P 2014Phys. Rev. B89035412] and its spin extension [Thonhauser Tet al2015Phys. Rev. Lett.115136402] represent the first examples of this class; in general, consistent vdW-DFs reflect a concerted expansion of a formal recast of the adiabatic-connection formula [Hyldgaard Pet al2014Phys. Rev. B90075148], an exponential summation of contributions to the local-field response, and the Dyson equation. We argue that the screening emphasis is essential because the exchange-correlation energy reflects an effective electrodynamics set by a long-range interaction. Two consequences are that (1) there are, in principle, no wiggle room in how one balances exchange and correlation, for example, in vdW-DF-cx, and that (2) consistent vdW-DFs have a formal structure that allows them to incorporate vertex-correction effects, at least in the case of levels that experience recoil-less interactions (for example, near the Fermi surface). We explore the extent to which the strictly nonempirical vdW-DF-cx formulation can serve as a systematic extension of the constraint-based semilocal functionals. For validation, we provide a complete survey of vdW-DF-cx performance for broad molecular processes, for the full set of 55 benchmarks in GMTKN55 [Goerigk Let al2017Phys. Chem. Chem. Phys.1932184] and comparing to the quantum-chemistry calculations that are summarized in that paper. We also provide new vdW-DF-cx results for metal surface energies and work functions that we compare to experiment. Finally, we use the screening insight to separate the vdW-DF nonlocal-correlation term into pure-vdW-interaction and local-field-susceptibility effects and present tools to compute and map the binding signatures.

Interactions of ions across carbon nanotubes

The interactions between a pair of Li⁺ ions across a semiconducting (8,0)CNT and a conducting (5,5)CNT has been investigated by density functional theory. The direct Coulomb interaction between the ions is almost completely screened. The band structure of the CNTs is not affected by the Li⁺ ions, but the Fermi level is raised to accommodate the extra electrons. Because of the unique band structure of CNTs this results in an effective attraction between the ions, which is greater for the (8,0)CNT. In contrast, a Cl^- ion inside a CNT forms a chemical bond which modifies the band structure. Again, the electrostatic field of the ion is almost completely screened outside of the tube. Nevertheless, the adsorption of a Li⁺ ion outside is favored by a Cl^- ion inside. This apparent attraction is mainly caused by a lowering of the work function of the CNT by the presence of the Cl^-.

Quantum mechanics of proteins in explicit water: The role of plasmon-like solute-solvent interactions

Quantum-mechanical van der Waals dispersion interactions play an essential role in intraprotein and protein-water interactions-the two main factors affecting the structure and dynamics of proteins in water. Typically, these interactions are only treated phenomenologically, via pairwise potential terms in classical force fields. Here, we use an explicit quantum-mechanical approach of density-functional tight-binding combined with the many-body dispersion formalism and demonstrate the relevance of many-body van der Waals forces both to protein energetics and to protein-water interactions. In contrast to commonly used pairwise approaches, many-body effects substantially decrease the relative stability of native states in the absence of water. Upon solvation, the protein-water dispersion interaction counteracts this effect and stabilizes native conformations and transition states. These observations arise from the highly delocalized and collective character of the interactions, suggesting a remarkable persistence of electron correlation through aqueous environments and providing the basis for long-range interaction mechanisms in biomolecular systems.

Impact of nuclear vibrations on van der Waals and Casimir interactions at zero and finite temperature

Recent advances in measuring van der Waals (vdW) interactions have probed forces on molecules at nanometric separations from metal surfaces and demonstrated the importance of infrared nonlocal polarization response and temperature effects, yet predictive theories for these systems remain lacking. We present a theoretical framework for computing vdW interactions among molecular structures, accounting for geometry, short-range electronic delocalization, dissipation, and collective nuclear vibrations (phonons) at atomic scales, along with long-range electromagnetic interactions in arbitrary macroscopic environments. We primarily consider experimentally relevant low-dimensional carbon allotropes, including fullerenes, carbyne, and graphene, and find that phonons couple strongly with long-range electromagnetic fields depending on molecular dimensionality and dissipation, especially at nanometric scales, creating delocalized phonon polaritons that substantially modify infrared molecular response. These polaritons, in turn, alter vdW interaction energies between molecular and macroscopic structures, producing nonmonotonic power laws and nontrivial temperature variations at nanometric separations feasible in current experiments.

Mimicking the Structure and Function of Ant Bridges in a Reconfigurable Microswarm for Electronic Applications

In nature, social insects are capable of self-organizing into various sophisticated and functional structures through local communications, which facilitate their cooperative accomplishment of complex tasks that are beyond the capabilities of individuals. Emulating this collective behavior in artificial robotic systems promises benefits in various engineering fields and has been partially realized through elaborate algorithm and physical designs. However, developing a swarm robotic system with group-level functionality at small scales remains a challenge. Herein, a microswarm system that mimics the structure and function of an ant bridge is realized by employing functionalized magnetite nanoparticles, which are paramagnetic and electrically conductive, as the building blocks. Through the application of a programmed oscillating magnetic field, the building blocks are reconfigured into a ribbon-like microswarm, which can perform reversible elongation with a high aspect ratio and, thus, is capable of constructing a conductive pathway for electrons between two disconnected electrodes with the bodies of functionalized nanoparticles. Furthermore, the microswarm is demonstrated to serve as a microswitch, repair broken microcircuits, and constitute flexible circuits, exhibiting a promising future for the practical applications in the electronics field.

A generally applicable atomic-charge dependent London dispersion correction

The so-called D4 model is presented for the accurate computation of London dispersion interactions in density functional theory approximations (DFT-D4) and generally for atomistic modeling methods. In this successor to the DFT-D3 model, the atomic coordination-dependent dipole polarizabilities are scaled based on atomic partial charges which can be taken from various sources. For this purpose, a new charge-dependent parameter-economic scaling function is designed. Classical charges are obtained from an atomic electronegativity equilibration procedure for which efficient analytical derivatives with respect to nuclear positions are developed. A numerical Casimir-Polder integration of the atom-in-molecule dynamic polarizabilities then yields charge- and geometry-dependent dipole-dipole dispersion coefficients. Similar to the D3 model, the dynamic polarizabilities are precomputed by time-dependent DFT and all elements up to radon (Z = 86) are covered. The two-body dispersion energy expression has the usual sum-over-atom-pairs form and includes dipole-dipole as well as dipole-quadrupole interactions. For a benchmark set of 1225 molecular dipole-dipole dispersion coefficients, the D4 model achieves an unprecedented accuracy with a mean relative deviation of 3.8% compared to 4.7% for D3. In addition to the two-body part, three-body effects are described by an Axilrod-Teller-Muto term. A common many-body dispersion expansion was extensively tested, and an energy correction based on D4 polarizabilities is found to be advantageous for larger systems. Becke-Johnson-type damping parameters for DFT-D4 are determined for more than 60 common density functionals. For various standard energy benchmark sets, DFT-D4 slightly but consistently outperforms DFT-D3. Especially for metal containing systems, the introduced charge dependence of the dispersion coefficients improves thermochemical properties. We suggest (DFT-)D4 as a physically improved and more sophisticated dispersion model in place of DFT-D3 for DFT calculations as well as other low-cost approaches like semi-empirical models.

Ubiquitous energy conversion of two-dimensional molecular crystals

Two-dimensional (2D) atomic crystals have triggered significant excitement due to their rich physics as well as potential industrial applications. The possibility of a molecular counterpart with scalable processability and superior performance is intriguing from both fundamental and applied perspectives. Here, we present the freestanding 2D molecular charge-transfer bis(ethylenedithio)tetrathiafulvalene-C₆₀ crystals prepared by a modified Langmuir-Blodgett method, with precisely controlled few-layer thickness and centimeter-scale lateral dimension. The interconversion of intrinsic excited process, the long-range ordering and anisotropic stacking arrangement of the molecular layered crystals generate external stimuli responsive behaviors and anisotropic spin-charge conversion with magnetic energy conversion ability, as well as a superior UV photosensitivity. Moreover, the 2D freestanding crystals demonstrate superior magneto-electrical properties. These results suggest that a new class of 2D atomically thin molecular crystals with novel electronic, optical and magnetic properties have great potential for spintronic, energy and sensor applications.

Micromechanical exfoliation of graphene on the atomistic scale

Mechanical exfoliation techniques are widely used to create high quality graphene samples for analytical use. Increasingly, mechanical methods are used to create large quantities of graphene, yet there is surprisingly little molecular insight into the mechanisms involved. We study the exfoliation of graphene with sticky tape using molecular dynamics. This is made possible by using a recently developed molecular dynamics forcefield, GraFF, to represent graphene's dispersion interactions. For nano-sized flakes we observe two different mechanisms depending on the polymer-adhesive used. A peeling mechanism which mixes shearing and normal mode exfoliation promotes synthesis of graphene rather than many-layered graphite. Armed with this new chemical insight we discuss the experimental methods that could preferentially produce graphene by mechanical exfoliation. We also introduce a mathematical model describing the repeated exfoliation of graphite.

Influence of Pore Size on the van der Waals Interaction in Two-Dimensional Molecules and Materials

Despite the importance of porous two-dimensional (2D) molecules and materials in advanced technological applications, the question of how the void space in these systems affects the van der Waals (vdW) scaling landscape has been largely unanswered. Analytical and numerical models presented herein demonstrate that the mere presence of a pore leads to markedly different vdW scaling across nonasymptotic distances, with certain relative pore sizes yielding effective power laws ranging from simple monotonic decay to the formation of minima, extended plateaus, and even maxima. These models are in remarkable agreement with first-principles approaches for the 2D building blocks of covalent organic frameworks (COFs), and reveal that COF macrocycle dimers and periodic bilayers exhibit unique vdW scaling behavior that is quite distinct from their nonporous analogs. These findings extend across a range of distances relevant to the nanoscale, and represent a hitherto unexplored avenue towards governing the self-assembly of complex nanostructures from porous 2D molecules and materials.

Are octahedral clusters missing on the carbon energy landscape?

We report a new class of carbon nanostructures at a lower sub-nano end of the size scale with a surprising stability, as compared to the well-known carbon fullerenes. The octahedral carbon clusters contain tetragonal rings, which, in spite of a common belief, prove to be an energy efficient means of plying graphene sheets to make three-dimensional spheroid shapes, similar to fullerenes. The two families of structures are shown to be competitive at small sizes (∼20 atoms) at room temperature, and for higher temperatures, at both small and large sizes (>200 atoms). Our calculations demonstrate that both vibrational and electronic spectra of these cluster families are similar, which thus might cloud their experimental identification. However, there is a sufficiently strong shift in vibrational frequencies below 160 and in the range of 600-800 cm^-1, which should help to identify different types of carbon clusters experimentally. We propose octahedral clusters and other structures containing tetragonal rings as viable structural elements and building units in inorganic chemistry and materials science of carbon along with fullerenes.

Theory and practice of modeling van der Waals interactions in electronic-structure calculations

The accurate description of long-range electron correlation, most prominently including van der Waals (vdW) dispersion interactions, represents a particularly challenging task in the modeling of molecules and materials. vdW forces arise from the interaction of quantum-mechanical fluctuations in the electronic charge density. Within (semi-)local density functional approximations or Hartree-Fock theory such interactions are neglected altogether. Non-covalent vdW interactions, however, are ubiquitous in nature and play a key role for the understanding and accurate description of the stability, dynamics, structure, and response properties in a plethora of systems. During the last decade, many promising methods have been developed for modeling vdW interactions in electronic-structure calculations. These methods include vdW-inclusive Density Functional Theory and correlated post-Hartree-Fock approaches. Here, we focus on the methods within the framework of Density Functional Theory, including non-local van der Waals density functionals, interatomic dispersion models within many-body and pairwise formulation, and random phase approximation-based approaches. This review aims to guide the reader through the theoretical foundations of these methods in a tutorial-style manner and, in particular, highlight practical aspects such as the applicability and the advantages and shortcomings of current vdW-inclusive approaches. In addition, we give an overview of complementary experimental approaches, and discuss tools for the qualitative understanding of non-covalent interactions as well as energy decomposition techniques. Besides representing a reference for the current state-of-the-art, this work is thus also designed as a concise and detailed introduction to vdW-inclusive electronic structure calculations for a general and broad audience.

A TPD-based determination of the graphite interlayer cohesion energy

Temperature Programmed Desorption (TPD) spectroscopy was used to determine the binding energies of polycyclic aromatic hydrocarbons C _n H _m (22 ≤ n ≤ 60) with highly oriented pyrolytic graphite. These energies were then used to estimate the dispersive graphite interlayer cohesion by means of a refined extrapolation method proposed by Björk et al. This yields a cohesion energy of 44.0 ± 3.8 meV per carbon atom. We discuss some limits of the TPD-based approach and contrast our values with previous determinations of the interlayer cohesion energy of graphite.

Dissipative Self-Assembly of Anisotropic Nanoparticle Chains with Combined Electrodynamic and Electrostatic Interactions

Dissipative self-assembly of colloidal nanoparticles offers the prospect of creating reconfigurable artificial materials and systems, yet the phenomenon only occurs far from thermodynamic equilibrium. Therefore, it is usually difficult to predict and control. Here, a dissipative colloidal solution system, where anisotropic chains with different interparticle separations in two perpendicular directions transiently arise among largely disordered silver nanoparticles illuminated by a laser beam, is reported. The optical field creates a nonequilibrium dissipative state, where a disorder-to-order transition occurs driven by anisotropic electrodynamic interactions coupled with electrostatic interactions. Investigation of the temporal dynamics and spatial arrangements of the nanoparticle system shows that the optical binding strength and entropy of the system are two crucial parameters for the formation of the anisotropic chains and responsible for adaptive behaviors, such as self-replication of dimer units. Formation of anisotropic nanoparticle chains is also observed among colloidal nanoparticles made from other metal (e.g., Au), polymer (e.g., polystyrene), ceramic (e.g., CeO₂ ), and hybrid materials (e.g., SiO₂ @Au core-shell), suggesting that light-driven self-organization will provide a wide range of opportunities to discover new dissipative structures under thermal fluctuations and build novel anisotropic materials with nanoscale order.