Unraveling H2 chemisorption and physisorption on metal decorated graphene using quantum Monte Carlo

Molecular hydrogen has the potential to significantly reduce the use of carbon dioxide emitting energy processes. However, hydrogen gas storage is a major bottleneck for its large-scale use as current storage methods are energy intensive. Among different storage methods, physisorbing molecular hydrogen at ambient pressure and temperatures is a promising alternative-particularly in light of the advancements in tunable lightweight nanomaterials and high throughput screening methods. Nonetheless, understanding hydrogen adsorption in well-defined nanomaterials remains experimentally challenging and reference information is scarce despite the proliferation of works predicting hydrogen adsorption. We focus on Li, Na, Ca, and K, decorated graphene sheets as substrates for molecular hydrogen adsorption, and compute the most accurate adsorption energies available to date using quantum diffusion Monte Carlo (DMC). Building on our previous insights at the density functional theory (DFT) level, we find that a weak covalent chemisorption of molecular hydrogen, known as Kubas interaction, is feasible on Ca decorated graphene according to DMC, in agreement with DFT. This finding is in contrast to previous DMC predictions of the 4H2/Ca+ gas cluster (without graphene) where chemisorption is not favored. However, we find that the adsorption energy of hydrogen on metal decorated graphene according to a widely used DFT method is not fully consistent with DMC. The reference adsorption energies reported herein can be used to find better work-horse methods for application in large-scale modeling of hydrogen adsorption. Furthermore, the implications of this work affect strategies for finding suitable hydrogen storage materials and high-throughput methods.

Thermodynamic and structural insights into the repurposing of drugs that bind to SARS-CoV-2 main protease

Although researchers have been working tirelessly since the COVID-19 outbreak, so far only three drugs - remdesivir, ronapreve and molnupiravir - have been approved for use in some countries which directly target the SARS-CoV-2 virus. Given the slow pace and substantial costs of new drug discovery and development, together with the urgency of the matter, repurposing of existing drugs for the ongoing disease is an attractive proposition. In a recent study, a high-throughput X-ray crystallographic screen was performed for a selection of drugs which have been approved or are in clinical trials. Thirty-seven compounds have been identified from drug libraries all of which bind to the SARS-CoV-2 main protease (3CL^pro). In the current study, we use molecular dynamics simulation and an ensemble-based free energy approach, namely, enhanced sampling of molecular dynamics with approximation of continuum solvent (ESMACS), to investigate a subset of the aforementioned compounds. The drugs studied here are highly diverse, interacting with different binding sites and/or subsites of 3CL^pro. The predicted free energies are compared with experimental results wherever they are available and they are found to be in excellent agreement. Our study also provides detailed energetic insights into the nature of the associated drug-protein binding, in turn shedding light on the design and discovery of potential drugs.

Enhancing the Accuracy of Ab Initio Molecular Dynamics by Fine Tuning of Effective Two-Body Interactions: Acetonitrile as a Test Case

Grimme's dispersion-corrected density functional theory (DFT-D) methods have emerged among the most practical approaches to perform accurate quantum mechanical calculations on molecular systems ranging from small clusters to microscopic and mesoscopic samples, i.e., including hundreds or thousands of molecules. Moreover, DFT-D functionals can be easily integrated into popular ab initio molecular dynamics (MD) software packages to carry out first-principles condensed-phase simulations at an affordable computational cost. Here, starting from the well-established D3 version of the dispersion-correction term, we present a simple protocol to improve the accurate description of the intermolecular interactions of molecular clusters of growing size, considering acetonitrile as a test case. Optimization of the interaction energy was performed with reference to diffusion quantum Monte Carlo calculations, successfully reaching the same inherent accuracy of the latter (statistical error of ∼0.1 kcal/mol per molecule). The refined DFT-D3 model was then used to perform ab initio MD simulations of liquid acetonitrile, again showing significant improvements toward available experimental data with respect to the default correction.

Pandemic drugs at pandemic speed: infrastructure for accelerating COVID-19 drug discovery with hybrid machine learning- and physics-based simulations on high-performance computers

The race to meet the challenges of the global pandemic has served as a reminder that the existing drug discovery process is expensive, inefficient and slow. There is a major bottleneck screening the vast number of potential small molecules to shortlist lead compounds for antiviral drug development. New opportunities to accelerate drug discovery lie at the interface between machine learning methods, in this case, developed for linear accelerators, and physics-based methods. The two in silico methods, each have their own advantages and limitations which, interestingly, complement each other. Here, we present an innovative infrastructural development that combines both approaches to accelerate drug discovery. The scale of the potential resulting workflow is such that it is dependent on supercomputing to achieve extremely high throughput. We have demonstrated the viability of this workflow for the study of inhibitors for four COVID-19 target proteins and our ability to perform the required large-scale calculations to identify lead antiviral compounds through repurposing on a variety of supercomputers.

Unusual reversibility in molecular break-up of PAHs: the case of pentacene dehydrogenation on Ir(111)

In this work, we characterise the adsorption of pentacene molecules on Ir(111) and their dissociation behaviour as a function of temperature.

Correction: Unusual reversibility in molecular break-up of PAHs: the case of pentacene dehydrogenation on Ir(111)

Correction for ‘Unusual reversibility in molecular break-up of PAHs: the case of pentacene dehydrogenation on Ir(111)’ by Davide Curcio et al., Chem. Sci., 2021, DOI: 10.1039/d0sc03734f.

Electronic correlations and transport in iron at Earth's core conditions

The transport properties of iron under Earth’s inner core conditions are essential input for the geophysical modelling but are poorly constrained experimentally. Here we show that the thermal and electrical conductivity of iron at those conditions remains high even if the electron-electron-scattering (EES) is properly taken into account. This result is obtained by ab initio simulations taking into account consistently both thermal disorder and electronic correlations. Thermal disorder suppresses the non-Fermi-liquid behavior of the body-centered cubic iron phase, hence, reducing the EES; the total calculated thermal conductivity of this phase is 220 Wm⁻¹ K⁻¹ with the EES reduction not exceeding 20%. The EES and electron-lattice scattering are intertwined resulting in breaking of the Matthiessen’s rule with increasing EES. In the hexagonal close-packed iron the EES is also not increased by thermal disorder and remains weak. Our main finding thus holds for the both likely iron phases in the inner core.

The heat and electrical conductivity of Earth’s core matter represent key input quantities for geophysical models of the Earth’s core evolution and geodynamo. Here, the authors show how the scattering due to interactions between electrons has a relatively weak impact on the electrical and thermal conductivities of iron at core conditions.

Translucency of Graphene to van der Waals Forces Applies to Atoms/Molecules with Different Polar Character

Graphene has been proposed to be either fully transparent to van der Waals interactions to the extent of allowing switching between hydrophobic and hydrophilic behavior, or partially transparent (translucent), yet there has been considerable debate on this topic, which is still ongoing. In a combined experimental and theoretical study we investigate the effects of different metal substrates on the adsorption energy of atomic (argon) and molecular (carbon monoxide) adsorbates on high-quality epitaxial graphene. We demonstrate that while the adsorption energy is certainly affected by the chemical composition of the supporting substrate and by the corrugation of the carbon lattice, the van der Waals interactions between adsorbates and the metal surfaces are partially screened by graphene. Our results indicate that the concept of graphene translucency, already introduced in the case of water droplets, is found to hold more generally also in the case of single polar molecules and atoms, which are apolar.

Physisorption of Water on Graphene: Subchemical Accuracy from Many-Body Electronic Structure Methods

Wet carbon interfaces are ubiquitous in the natural world and exhibit anomalous properties, which could be exploited by emerging technologies. However, progress is limited by lack of understanding at the molecular level. Remarkably, even for the most fundamental system (a single water molecule interacting with graphene), there is no consensus on the nature of the interaction. We tackle this by performing an extensive set of complementary state-of-the-art computer simulations on some of the world's largest supercomputers. From this effort a consensus on the water-graphene interaction strength has been obtained. Our results have significant impact for the physical understanding, as they indicate that the interaction is weaker than predicted previously. They also pave the way for more accurate and reliable studies of liquid water at carbon interfaces.

Performance of van der Waals Corrected Functionals for Guest Adsorption in the M2(dobdc) Metal-Organic Frameworks

Small-molecule binding in metal-organic frameworks (MOFs) can be accurately studied both experimentally and computationally, provided the proper tools are employed. Herein, we compare and contrast properties associated with guest binding by means of density functional theory (DFT) calculations using nine different functionals for the M2(dobdc) (dobdc4- = 2,5-dioxido,1,4-benzenedicarboxylate) series, where M = Mg, Mn, Fe, Co, Ni, Cu, and Zn. Additionally, we perform Quantum Monte Carlo (QMC) calculations for one system to determine if this method can be used to assess the performance of DFT. We also make comparisons with previously published experimental results for carbon dioxide and water and present new methane neutron powder diffraction (NPD) data for further comparison. All of the functionals are able to predict the experimental variation in the binding energy from one metal to the next; however, the interpretation of the performance of the functionals depends on which value is taken as the reference. On the one hand, if we compare against experimental values, we would conclude that the optB86b-vdW and optB88-vdW functionals systematically overestimate the binding strength, while the second generation of van der Waals (vdW) nonlocal functionals (vdw-DF2 and rev-vdW-DF2) correct for this providing a good description of binding energies. On the other hand, if the QMC calculation is taken as the reference then all of the nonlocal functionals yield results that fall just outside the error of the higher-level calculation. The empirically corrected vdW functionals are in reasonable agreement with experimental heat of adsorptions but under bind when compared with QMC, while Perdew-Burke-Ernzerhof fails by more than 20 kJ/mol regardless of which reference is employed. All of the functionals, with the exception of vdW-DF2, predict reasonable framework and guest binding geometries when compared with NPD measurements. The newest of the functionals considered, rev-vdW-DF2, should be used in place of vdW-DF2, as it yields improved bond distances with similar quality binding energies.

A comparison between quantum chemistry and quantum Monte Carlo techniques for the adsorption of water on the (001) LiH surface

We present a comprehensive benchmark study of the adsorption energy of a single water molecule on the (001) LiH surface using periodic coupled cluster and quantum Monte Carlo theories. We benchmark and compare different implementations of quantum chemical wave function based theories in order to verify the reliability of the predicted adsorption energies and the employed approximations. Furthermore we compare the predicted adsorption energies to those obtained employing widely used van der Waals density-functionals. Our findings show that quantum chemical approaches are becoming a robust and reliable tool for condensed phase electronic structure calculations, providing an additional tool that can also help in potentially improving currently available van der Waals density-functionals.

Nanoscale π-π stacked molecules are bound by collective charge fluctuations

Non-covalent π−π interactions are central to chemical and biological processes, yet the full understanding of their origin that would unite the simplicity of empirical approaches with the accuracy of quantum calculations is still missing. Here we employ a quantum-mechanical Hamiltonian model for van der Waals interactions, to demonstrate that intermolecular electron correlation in large supramolecular complexes at equilibrium distances is appropriately described by collective charge fluctuations. We visualize these fluctuations and provide connections both to orbital-based approaches to electron correlation, as well as to the simple London pairwise picture. The reported binding energies of ten supramolecular complexes obtained from the quantum-mechanical fluctuation model joined with density functional calculations are within 5% of the reference energies calculated with the diffusion quantum Monte-Carlo method. Our analysis suggests that π−π stacking in supramolecular complexes can be characterized by strong contributions to the binding energy from delocalized, collective charge fluctuations—in contrast to complexes with other types of bonding.

Attractive, non-covalent interactions between aromatic rings—termed π−π stacking—is common in chemistry but difficult to model. Here the authors report a quantum-mechanical model to show the importance of collective charge fluctuations for understanding pi-stacked supramolecular systems.

Toward Accurate Adsorption Energetics on Clay Surfaces

Clay minerals are ubiquitous in nature, and the manner in which they interact with their surroundings has important industrial and environmental implications. Consequently, a molecular-level understanding of the adsorption of molecules on clay surfaces is crucial. In this regard computer simulations play an important role, yet the accuracy of widely used empirical force fields (FF) and density functional theory (DFT) exchange-correlation functionals is often unclear in adsorption systems dominated by weak interactions. Herein we present results from quantum Monte Carlo (QMC) for water and methanol adsorption on the prototypical clay kaolinite. To the best of our knowledge, this is the first time QMC has been used to investigate adsorption at a complex, natural surface such as a clay. As well as being valuable in their own right, the QMC benchmarks obtained provide reference data against which the performance of cheaper DFT methods can be tested. Indeed using various DFT exchange-correlation functionals yields a very broad range of adsorption energies, and it is unclear a priori which evaluation is better. QMC reveals that in the systems considered here it is essential to account for van der Waals (vdW) dispersion forces since this alters both the absolute and relative adsorption energies of water and methanol. We show, via FF simulations, that incorrect relative energies can lead to significant changes in the interfacial densities of water and methanol solutions at the kaolinite interface. Despite the clear improvements offered by the vdW-corrected and the vdW-inclusive functionals, absolute adsorption energies are often overestimated, suggesting that the treatment of vdW forces in DFT is not yet a solved problem.

Saturation of electrical resistivity of solid iron at Earth's core conditions

We report on the temperature dependence of the electrical resistivity of solid iron at high pressure, up to and including conditions likely to be found at the centre of the Earth. We have extended some of the calculations of the resistivities of pure solid iron we recently performed at Earth's core conditions (Pozzo et al. in Earth Planet Sci Lett 393:159-164, 2014) to lower temperature. We show that at low temperature the resistivity increases linearly with temperature, and saturates at high temperature. This saturation effect is well known as the Mott-Ioffe-Regel limit in metals, but has been largely ignored to estimate the resistivity of iron at Earth's core conditions. Recent experiments (Gomi et al. in Phys Earth Planet Int 224:88-103, 2013) coupled new high pressure data and saturation to predict the resitivity of iron and iron alloys at Earth's core conditions, and reported values up to three times lower than previous estimates, confirming recent first principles calculations (de Koker et al. in Proc Natl Acad Sci 109:4070-4073, 2012; Pozzo et al. in Nature 485:355-358, 2012, Phys Rev B 87:014110-10, 2013, Earth Planet Sci Lett 393:159-164, 2014; Davies et al. in Nat Geosci 8:678-685, 2015). The present results support the saturation effect idea.

Molecular Lifting, Twisting, and Curling during Metal-Assisted Polycyclic Hydrocarbon Dehydrogenation

The atomistic understanding of the dissociation mechanisms for large molecules adsorbed on surfaces is still a challenge in heterogeneous catalysis. This is especially true for polycyclic aromatic hydrocarbons, which represent an important class of organic compounds used to produce novel graphene-based architectures. Here, we show that coronene molecules adsorbed on Ir(111) undergo major conformational changes during dissociation. They first tilt upward with respect to the surface, still keeping their planar configuration, and subsequently experience a rotation, which changes the molecular axis orientation. Upon lifting, the internal C-C strain is initially relieved; as the dehydrogenation proceeds, the molecules experience a progressive increase in the average interatomic distance and gradually settle to form dome-shaped nanographene flakes. Our results provide important insight into the complex mechanism of molecular breakup, which could have implications in the synthesis of new carbon-based nanostructured materials.

Energy benchmarks for methane-water systems from quantum Monte Carlo and second-order Møller-Plesset calculations

The quantum Monte Carlo (QMC) technique is used to generate accurate energy benchmarks for methane-water clusters containing a single methane monomer and up to 20 water monomers. The benchmarks for each type of cluster are computed for a set of geometries drawn from molecular dynamics simulations. The accuracy of QMC is expected to be comparable with that of coupled-cluster calculations, and this is confirmed by comparisons for the CH4-H2O dimer. The benchmarks are used to assess the accuracy of the second-order Møller-Plesset (MP2) approximation close to the complete basis-set limit. A recently developed embedded many-body technique is shown to give an efficient procedure for computing basis-set converged MP2 energies for the large clusters. It is found that MP2 values for the methane binding energies and the cohesive energies of the water clusters without methane are in close agreement with the QMC benchmarks, but the agreement is aided by partial cancelation between 2-body and beyond-2-body errors of MP2. The embedding approach allows MP2 to be applied without loss of accuracy to the methane hydrate crystal, and it is shown that the resulting methane binding energy and the cohesive energy of the water lattice agree almost exactly with recently reported QMC values.

Water on BN doped benzene: a hard test for exchange-correlation functionals and the impact of exact exchange on weak binding

Density functional theory (DFT) studies of weakly interacting complexes have recently focused on the importance of van der Waals dispersion forces, whereas the role of exchange has received far less attention. Here, by exploiting the subtle binding between water and a boron and nitrogen doped benzene derivative (1,2-azaborine) we show how exact exchange can alter the binding conformation within a complex. Benchmark values have been calculated for three orientations of the water monomer on 1,2-azaborine from explicitly correlated quantum chemical methods, and we have also used diffusion quantum Monte Carlo. For a host of popular DFT exchange-correlation functionals we show that the lack of exact exchange leads to the wrong lowest energy orientation of water on 1,2-azaborine. As such, we suggest that a high proportion of exact exchange and the associated improvement in the electronic structure could be needed for the accurate prediction of physisorption sites on doped surfaces and in complex organic molecules. Meanwhile to predict correct absolute interaction energies an accurate description of exchange needs to be augmented by dispersion inclusive functionals, and certain non-local van der Waals functionals (optB88- and optB86b-vdW) perform very well for absolute interaction energies. Through a comparison with water on benzene and borazine (B3N3H6) we show that these results could have implications for the interaction of water with doped graphene surfaces, and suggest a possible way of tuning the interaction energy.

Communication: Water on hexagonal boron nitride from diffusion Monte Carlo

Despite a recent flurry of experimental and simulation studies, an accurate estimate of the interaction strength of water molecules with hexagonal boron nitride is lacking. Here, we report quantum Monte Carlo results for the adsorption of a water monomer on a periodic hexagonal boron nitride sheet, which yield a water monomer interaction energy of -84 ± 5 meV. We use the results to evaluate the performance of several widely used density functional theory (DFT) exchange correlation functionals and find that they all deviate substantially. Differences in interaction energies between different adsorption sites are however better reproduced by DFT.

Communication: On the stability of ice 0, ice i, and I(h)

Using ab initio methods, we examine the stability of ice 0, a recently proposed tetragonal form of ice implicated in the homogeneous freezing of water [J. Russo, F. Romano, and H. Tanaka, Nat. Mater. 13, 670 (2014)]. Vibrational frequencies are computed across the complete Brillouin Zone using Density Functional Theory (DFT), to confirm mechanical stability and quantify the free energy of ice 0 relative to ice I(h). The robustness of this result is tested via dispersion corrected semi-local and hybrid DFT, and Quantum Monte-Carlo calculation of lattice energies. Results indicate that popular molecular models only slightly overestimate the stability of ice zero. In addition, we study all possible realisations of proton disorder within the ice zero unit cell, and identify the ground state as ferroelectric. Comparisons are made to other low density metastable forms of ice, suggesting that the ice i structure [C. J. Fennel and J. D. Gezelter, J. Chem. Theory Comput. 1, 662 (2005)] may be equally relevant to ice formation.

Structure of magnesium selenate enneahydrate, MgSeO4·9H2O, from 5 to 250 K using neutron time-of-flight Laue diffraction

The complete structure of MgSeO4·9H2O has been refined from neutron single-crystal diffraction data obtained at 5, 100, 175 and 250 K. It is monoclinic, space group P2₁/c, Z = 4, with unit-cell parameters a = 7.222 (2), b = 10.484 (3), c = 17.327 (4) Å, β = 109.57 (2)°, and V = 1236.1 (6) Å(3) [ρ(calc) = 1770 (1) kg m(-3)] at 5 K. The structure consists of isolated [Mg(H2O)6](2+) octahedra, [SeO4](2-) tetrahedra and three interstitial lattice water molecules, all on sites of symmetry 1. The positions of the H atoms agree well with those inferred on the basis of geometrical considerations in the prior X-ray powder diffraction structure determination: no evidence of orientational disorder of the water molecules is apparent in the temperature range studied. Six of the nine water molecules are hydrogen bonded to one another to form a unique centrosymmetric dodecamer, (H2O)12. Raman spectra have been acquired in the range 170-4000 cm(-1) at 259 and 78 K; ab initio calculations, using density functional theory, have been carried out in order to aid in the analysis of the Raman spectrum as well as providing additional insights into the geometry and thermodynamics of the hydrogen bonds. Complementary information concerning the thermal expansion, crystal morphology and the solubility are also presented.

On the room-temperature phase diagram of high pressure hydrogen: an ab initio molecular dynamics perspective and a diffusion Monte Carlo study

The finite-temperature phase diagram of hydrogen in the region of phase IV and its neighborhood was studied using the ab initio molecular dynamics (MD) and the ab initio path-integral molecular dynamics (PIMD). The electronic structures were analyzed using the density-functional theory (DFT), the random-phase approximation, and the diffusion Monte Carlo (DMC) methods. Taking the state-of-the-art DMC results as benchmark, comparisons of the energy differences between structures generated from the MD and PIMD simulations, with molecular and dissociated hydrogens, respectively, in the weak molecular layers of phase IV, indicate that standard functionals in DFT tend to underestimate the dissociation barrier of the weak molecular layers in this mixed phase. Because of this underestimation, inclusion of the quantum nuclear effects (QNEs) in PIMD using electronic structures generated with these functionals leads to artificially dissociated hydrogen layers in phase IV and an error compensation between the neglect of QNEs and the deficiencies of these functionals in standard ab initio MD simulations exists. This analysis partly rationalizes why earlier ab initio MD simulations complement so well the experimental observations. The temperature and pressure dependencies for the stability of phase IV were also studied in the end and compared with earlier results.

Benchmarking the performance of density functional theory and point charge force fields in their description of sI methane hydrate against diffusion Monte Carlo

High quality reference data from diffusion Monte Carlo calculations are presented for bulk sI methane hydrate, a complex crystal exhibiting both hydrogen-bond and dispersion dominated interactions. The performance of some commonly used exchange-correlation functionals and all-atom point charge force fields is evaluated. Our results show that none of the exchange-correlation functionals tested are sufficient to describe both the energetics and the structure of methane hydrate accurately, while the point charge force fields perform badly in their description of the cohesive energy but fair well for the dissociation energetics. By comparing to ice Ih, we show that a good prediction of the volume and cohesive energies for the hydrate relies primarily on an accurate description of the hydrogen bonded water framework, but that to correctly predict stability of the hydrate with respect to dissociation to ice Ih and methane gas, accuracy in the water-methane interaction is also required. Our results highlight the difficulty that density functional theory faces in describing both the hydrogen bonded water framework and the dispersion bound methane.

Structure, hydrogen bonding and thermal expansion of ammonium carbonate monohydrate

We have determined the crystal structure of ammonium carbonate monohydrate, (NH4)2CO3·H2O, using Laue single-crystal diffraction methods with pulsed neutron radiation. The crystal is orthorhombic, space group Pnma (Z = 4), with unit-cell dimensions a = 12.047 (3), b = 4.453 (1), c = 11.023 (3) Å and V = 591.3 (3) Å(3) [ρcalc = 1281.8 (7) kg m(-3)] at 10 K. The single-crystal data collected at 10 and 100 K are complemented by X-ray powder diffraction data measured from 245 to 273 K, Raman spectra measured from 80 to 263 K and an athermal zero-pressure calculation of the electronic structure and phonon spectrum carried out using density functional theory (DFT). We find no evidence of a phase transition between 10 and 273 K; above 273 K, however, the title compound transforms first to ammonium sesquicarbonate monohydrate and subsequently to ammonium bicarbonate. The crystallographic and spectroscopic data and the calculations reveal a quite strongly hydrogen-bonded structure (EHB ≃ 30-40 kJ mol(-1)), on the basis of H...O bond lengths and the topology of the electron density at the bond critical points, in which there is no free rotation of the ammonium cation at any temperature. The barrier to free rotation of the ammonium ions is estimated from the observed librational frequency to be ∼ 36 kJ mol(-1). The c-axis exhibits negative thermal expansion, but the thermal expansion behaviour of the a and b axes is ormal.

Cooperative interplay of van der Waals forces and quantum nuclear effects on adsorption: H at graphene and at coronene

The energetic barriers that atoms and molecules often experience when binding to surfaces are incredibly important to a myriad of chemical and physical processes. However, these barriers are difficult to describe accurately with current computer simulation approaches. Two prominent contemporary challenges faced by simulation are the role of van der Waals forces and nuclear quantum effects. Here we examine the widely studied model systems of hydrogen on graphene and coronene using a van der Waals inclusive density functional theory approach together with path integral molecular dynamics at 50 K. We find that both van der Waals and quantum nuclear effects work together in a cooperative manner to dramatically reduce the barriers for hydrogen atoms to adsorb. This suggests that the low temperature hydrogenation of graphene is easier than previously thought and in more general terms that the combined roles of van der Waals and quantum tunnelling can lead to qualitative changes in adsorption.

Anharmonicity, mechanical instability, and thermodynamic properties of the Cr-Re σ-phase

Using density-functional theory in combination with the direct force method and molecular dynamics we investigate the vibrational properties of a binary Cr-Re σ-phase. In the harmonic approximation, we have computed phonon dispersion curves and density of states, evidencing structural and chemical effects. We found that the σ-phase is mechanically unstable in some configurations, for example, when all crystallographic sites are occupied by Re atoms. By using a molecular-dynamics-based method, we have analysed the anharmonicity in the system and found negligible effects (~0.5 kJ/mol) on the Helmholtz energy of the binary Cr-Re σ-phase up to 2000 K (~0.8T(m)). Finally, we show that the vibrational contribution has significant consequences on the disordering of the σ-phase at high temperature.

Hard Numbers for Large Molecules: Toward Exact Energetics for Supramolecular Systems

Noncovalent interactions are ubiquitous in molecular and condensed-phase environments, and hence a reliable theoretical description of these fundamental interactions could pave the way toward a more complete understanding of the microscopic underpinnings for a diverse set of systems in chemistry and biology. In this work, we demonstrate that recent algorithmic advances coupled to the availability of large-scale computational resources make the stochastic quantum Monte Carlo approach to solving the Schrödinger equation an optimal contender for attaining "chemical accuracy" (1 kcal/mol) in the binding energies of supramolecular complexes of chemical relevance. To illustrate this point, we considered a select set of seven host-guest complexes, representing the spectrum of noncovalent interactions, including dispersion or van der Waals forces, π-π stacking, hydrogen bonding, hydrophobic interactions, and electrostatic (ion-dipole) attraction. A detailed analysis of the interaction energies reveals that a complete theoretical description necessitates treatment of terms well beyond the standard London and Axilrod-Teller contributions to the van der Waals dispersion energy.

Fine tuning of graphene-metal adhesion by surface alloying

We show that bimetallic surface alloying provides a viable route for governing the interaction between graphene and metal through the selective choice of the elemental composition of the surface alloy. This concept is illustrated by an experimental and theoretical characterization of the properties of graphene on a model PtRu surface alloy on Ru(0001), with a concentration of Pt atoms in the first layer between 0 and 50%. The progressive increase of the Pt content determines the gradual detachment of graphene from the substrate, which results from the modification of the carbon orbital hybridization promoted by Pt. Alloying is also found to affect the morphology of graphene, which is strongly corrugated on bare Ru, but becomes flat at a Pt coverage of 50%. The method here proposed can be readily extended to several supports, thus opening the way to the conformal growth of graphene on metals and to a full tunability of the graphene-substrate interaction.

Local electronic structure and density of edge and facet atoms at Rh nanoclusters self-assembled on a graphene template

The chemical and physical properties of nanoclusters largely depend on their sizes and shapes. This is partly due to finite size effects influencing the local electronic structure of the nanocluster atoms which are located on the nanofacets and on their edges. Here we present a thorough study on graphene-supported Rh nanocluster assemblies and their geometry-dependent electronic structure obtained by combining high-energy resolution core level photoelectron spectroscopy, scanning tunneling microscopy, and density functional theory. We demonstrate the possibility to finely control the morphology and the degree of structural order of Rh clusters grown in register with the template surface of graphene/Ir(111). By comparing measured and calculated core electron binding energies, we identify edge, facet, and bulk atoms of the nanoclusters. We describe how small interatomic distance changes occur while varying the nanocluster size, substantially modifying the properties of surface atoms. The properties of under-coordinated Rh atoms are discussed in view of their importance in heterogeneous catalysis and magnetism.

The kinetics of homogeneous melting beyond the limit of superheating

Molecular dynamics simulation is used to study the time-scales involved in the homogeneous melting of a superheated crystal. The interaction model used is an embedded-atom model for Fe developed in previous work, and the melting process is simulated in the microcanonical (N, V, E) ensemble. We study periodically repeated systems containing from 96 to 7776 atoms, and the initial system is always the perfect crystal without free surfaces or other defects. For each chosen total energy E and number of atoms N, we perform several hundred statistically independent simulations, with each simulation lasting for between 500 ps and 10 ns, in order to gather statistics for the waiting time τ(w) before melting occurs. We find that the probability distribution of τ(w) is roughly exponential, and that the mean value <τ(w)> depends strongly on the excess of the initial steady temperature of the crystal above the superheating limit identified by other researchers. The mean <τ(w)> also depends strongly on system size in a way that we have quantified. For very small systems of ~100 atoms, we observe a persistent alternation between the solid and liquid states, and we explain why this happens. Our results allow us to draw conclusions about the reliability of the recently proposed Z method for determining the melting properties of simulated materials and to suggest ways of correcting for the errors of the method.