Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory-Based Ab Initio Molecular Dynamics II: Extensions to the Isobaric-Isoenthalpic and Isobaric-Isothermal Ensembles.

In the previous paper of this series [Ko, H.-Y. et al. J. Chem. Theory Comput. 2020, 16, 3757-3785], we presented a theoretical and algorithmic framework based on a localized representation of the occupied space that exploits the inherent sparsity in the real-space evaluation of the exact exchange (EXX) interaction in finite-gap systems. This was accompanied by a detailed description of exx, a massively parallel hybrid message-passing interface MPI/OpenMP implementation of this approach in Quantum ESPRESSO (QE) that enables linear scaling hybrid density functional theory (DFT)-based ab initio molecular dynamics (AIMD) in the microcanonical/canonical (NVE/NVT) ensembles of condensed-phase systems containing 500-1000 atoms (in fixed orthorhombic cells) with a wall time cost comparable to semi-local DFT. In this work, we extend the current capabilities of exx to enable hybrid DFT-based AIMD simulations of large-scale condensed-phase systems with general and fluctuating cells in the isobaric-isoenthalpic/isobaric-isothermal (NpH/NpT) ensembles. The theoretical extensions to this approach include an analytical derivation of the EXX contribution to the stress tensor for systems in general simulation cells with a computational complexity that scales linearly with system size. The corresponding algorithmic extensions to exx include optimized routines that (i) handle both static and fluctuating simulation cells with non-orthogonal lattice symmetries, (ii) solve Poisson's equation in general/non-orthogonal cells via an automated selection of the auxiliary grid directions in the Natan-Kronik representation of the discrete Laplacian operator, and (iii) evaluate the EXX contribution to the stress tensor. Using this approach, we perform a case study on a variety of condensed-phase systems (including liquid water, a benzene molecular crystal polymorph, and semi-conducting crystalline silicon) and demonstrate that the EXX contributions to the energy and stress tensor simultaneously converge with an appropriate choice of exx parameters. This is followed by a critical assessment of the computational performance of the extended exx module across several different high-performance computing architectures via case studies on (i) the computational complexity due to lattice symmetry during NpT simulations of three different ice polymorphs (i.e., ice Ih, II, and III) and (ii) the strong/weak parallel scaling during large-scale NpT simulations of liquid water. We demonstrate that the robust and highly scalable implementation of this approach in the extended exx module is capable of evaluating the EXX contribution to the stress tensor with negligible cost (<1%) as well as all other EXX-related quantities needed during NpT simulations of liquid water (with a very tight 150 Ry planewave cutoff) in ≈5.2 s ((H2O)128) and ≈6.8 s ((H2O)256) per AIMD step. As such, the extended exx module presented in this work brings us another step closer to routinely performing hybrid DFT-based AIMD simulations of sufficient duration for large-scale condensed-phase systems across a wide range of thermodynamic conditions.

Self-interaction correction in water-ion clusters.

We study the importance of self-interaction errors in density functional approximations for various water-ion clusters. We have employed the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method in conjunction with the local spin-density approximation, Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA), and strongly constrained and appropriately normed (SCAN) meta-GGA to describe binding energies of hydrogen-bonded water-ion clusters, i.e., water-hydronium, water-hydroxide, water-halide, and non-hydrogen-bonded water-alkali clusters. In the hydrogen-bonded water-ion clusters, the building blocks are linked by hydrogen atoms, although the links are much stronger and longer-ranged than the normal hydrogen bonds between water molecules because the monopole on the ion interacts with both permanent and induced dipoles on the water molecules. We find that self-interaction errors overbind the hydrogen-bonded water-ion clusters and that FLOSIC reduces the error and brings the binding energies into closer agreement with higher-level calculations. The non-hydrogen-bonded water-alkali clusters are not significantly affected by self-interaction errors. Self-interaction corrected PBE predicts the lowest mean unsigned error in binding energies (≤50 meV/H2O) for hydrogen-bonded water-ion clusters. Self-interaction errors are also largely dependent on the cluster size, and FLOSIC does not accurately capture the subtle variation in all clusters, indicating the need for further refinement.

Exploring and enhancing the accuracy of interior-scaled Perdew-Zunger self-interaction correction.

The Perdew-Zunger self-interaction correction (PZ-SIC) improves the performance of density functional approximations for the properties that involve significant self-interaction error (SIE), as in stretched bond situations, but overcorrects for equilibrium properties where SIE is insignificant. This overcorrection is often reduced by local scaling self-interaction correction (LSIC) of the PZ-SIC to the local spin density approximation (LSDA). Here, we propose a new scaling factor to use in an LSIC-like approach that satisfies an additional important constraint: the correct coefficient of the atomic number Z in the asymptotic expansion of the exchange-correlation (xc) energy for atoms. LSIC and LSIC+ are scaled by functions of the iso-orbital indicator zσ, which distinguishes one-electron regions from many-electron regions. LSIC+ applied to the LSDA works better for many equilibrium properties than LSDA-LSIC and the Perdew, Burke, and Ernzerhof generalized gradient approximation (GGA), and almost close to the strongly constrained and appropriately normed (SCAN) meta-GGA. LSDA-LSIC and LSDA-LSIC+, however, fail to predict interaction energies involving weaker bonds, in sharp contrast to their earlier successes. It is found that more than one set of localized SIC orbitals can yield a nearly degenerate energetic description of the same multiple covalent bond, suggesting that a consistent chemical interpretation of the localized orbitals requires a new way to choose their Fermi orbital descriptors. To make a locally scaled down SIC to functionals beyond the LSDA requires a gauge transformation of the functional's energy density. The resulting SCAN-sdSIC, evaluated on SCAN-SIC total and localized orbital densities, leads to an acceptable description of many equilibrium properties including the dissociation energies of weak bonds.

The Fermi-Löwdin self-interaction correction for ionization energies of organic molecules.

(Semi)-local density functional approximations (DFAs) suffer from self-interaction error (SIE). When the first ionization energy (IE) is computed as the negative of the highest-occupied orbital (HO) eigenvalue, DFAs notoriously underestimate them compared to quasi-particle calculations. The inaccuracy for the HO is attributed to SIE inherent in DFAs. We assessed the IE based on Perdew-Zunger self-interaction correction on 14 small to moderate-sized organic molecules relevant in organic electronics and polymer donor materials. Although self-interaction corrected DFAs were found to significantly improve the IE relative to the uncorrected DFAs, they overestimate. However, when the self-interaction correction is interiorly scaled using a function of the iso-orbital indicator zσ, only the regions where SIE is significant get a correction. We discuss these approaches and show how these methods significantly improve the description of the HO eigenvalue for the organic molecules.

Simple hydrogenic estimates for the exchange and correlation energies of atoms and atomic ions, with implications for density functional theory.

Exact density functionals for the exchange and correlation energies are approximated in practical calculations for the ground-state electronic structure of a many-electron system. An important exact constraint for the construction of approximations is to recover the correct non-relativistic large-Z expansions for the corresponding energies of neutral atoms with atomic number Z and electron number N = Z, which are correct to the leading order (-0.221Z5/3 and -0.021Z ln Z, respectively) even in the lowest-rung or local density approximation. We find that hydrogenic densities lead to Ex(N, Z) ≈ -0.354N2/3Z (as known before only for Z ≫ N ≫ 1) and Ec ≈ -0.02N ln N. These asymptotic estimates are most correct for atomic ions with large N and Z ≫ N, but we find that they are qualitatively and semi-quantitatively correct even for small N and N ≈ Z. The large-N asymptotic behavior of the energy is pre-figured in small-N atoms and atomic ions, supporting the argument that widely predictive approximate density functionals should be designed to recover the correct asymptotics. It is shown that the exact Kohn-Sham correlation energy, when calculated from the pure ground-state wavefunction, should have no contribution proportional to Z in the Z → ∞ limit for any fixed N.

A step in the direction of resolving the paradox of Perdew-Zunger self-interaction correction. II. Gauge consistency of the energy density at three levels of approximation.

The Perdew-Zunger (PZ) self-interaction correction (SIC) was designed to correct the one-electron limit of any approximate density functional for the exchange-correlation (xc) energy, while yielding no correction to the exact functional. Unfortunately, it spoils the slowly varying (in space) limits of the uncorrected approximate functionals, where those functionals are right by construction. The right limits can be restored by locally scaling down the energy density of the PZ SIC in many-electron regions, but then a spurious correction to the exact functional would be found unless the self-Hartree and exact self-xc terms of the PZ SIC energy density were expressed in the same gauge. Only the local density approximation satisfies the same-gauge condition for the energy density, which explains why the recent local-scaling SIC is found here to work excellently for atoms and molecules only with this basic approximation and not with the more advanced generalized gradient approximations (GGAs) and meta-GGAs, which lose the Hartree gauge via simplifying integrations by parts. The transformation of energy density that achieves the Hartree gauge for the exact xc functional can also be applied to approximate functionals. Doing so leads to a simple scaled-down self-interaction correction that is typically much more accurate than PZ SIC in tests for many molecular properties (including equilibrium bond lengths). The present work unambiguously shows that the largest errors of PZ SIC applied to standard functionals at three levels of approximation can be removed by restoring their correct slowly varying density limits. It also confirms the relevance of these limits to atoms and molecules.

Self-interaction error overbinds water clusters but cancels in structural energy differences.

We gauge the importance of self-interaction errors in density functional approximations (DFAs) for the case of water clusters. To this end, we used the Fermi-Löwdin orbital self-interaction correction method (FLOSIC) to calculate the binding energy of clusters of up to eight water molecules. Three representative DFAs of the local, generalized gradient, and metageneralized gradient families [i.e., local density approximation (LDA), Perdew-Burke-Ernzerhof (PBE), and strongly constrained and appropriately normed (SCAN)] were used. We find that the overbinding of the water clusters in these approximations is not a density-driven error. We show that, while removing self-interaction error does not alter the energetic ordering of the different water isomers with respect to the uncorrected DFAs, the resulting binding energies are corrected toward accurate reference values from higher-level calculations. In particular, self-interaction-corrected SCAN not only retains the correct energetic ordering for water hexamers but also reduces the mean error in the hexamer binding energies to less than 14 meV/[Formula: see text] from about 42 meV/[Formula: see text] for SCAN. By decomposing the total binding energy into many-body components, we find that large errors in the two-body interaction in SCAN are significantly reduced by self-interaction corrections. Higher-order many-body errors are small in both SCAN and self-interaction-corrected SCAN. These results indicate that orbital-by-orbital removal of self-interaction combined with a proper DFA can lead to improved descriptions of water complexes.

Molecular enneanuclear CuII phosphates containing planar hexanuclear and trinuclear sub-units: syntheses, structures, and magnetism.

Highly symmetric enneanuclear copper(ii) phosphates [Cu9(Pz)6(µ-OH)3(µ3-OH)(ArOPO3)4(DMF)3] (PzH = pyrazole, Ar = 2,6-(CHPh2)2-4-R-C6H2; R = Me, 2MeAr; Et, 2EtAr; iPr, 2iPrAr; and Ar = 2,6-iPr2C6H3, 2Dip) comprising nine copper(ii) centers and pyrazole, hydroxide and DMF as ancillary ligands were synthesized by a reaction involving the arylphosphate monoester, 1, copper(i)chloride, pyrazole, and triethylamine in a 4 : 9 : 6 : 14 ratio. All four complexes were characterized by single crystal structural analysis. The complexes contain two distinct structural motifs within the multinuclear copper scaffold: a hexanuclear unit and a trinuclear unit. In the latter, the three Cu(ii) centres are bridged by a µ3-OH. Each pair of Cu(ii) centers in the trinuclear unit are bridged by a pyrazole ligand. The hexanuclear unit is made up of three dinuclear Cu(ii) motifs where the two Cu(ii) centres are bridged by an -OH and a pyrazole ligand. The three dinuclear units are connected to each other by phosphate ligands. The latter also aid the fusion of the trinuclear and the hexanuclear motifs. Magnetic studies reveal a strong antiferromagnetic exchange between the Cu(ii) centres of the dinuclear units in the hexanuclear part and a strong spin frustration in the trinuclear part leading to a degenerate ground state.

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance.

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car-Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson's equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. This is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500-1000 atoms (e.g., (H2O)256) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

A step in the direction of resolving the paradox of Perdew-Zunger self-interaction correction.

Self-interaction (SI) error, which results when exchange-correlation contributions to the total energy are approximated, limits the reliability of many density functional approximations. The Perdew-Zunger SI correction (PZSIC), when applied in conjunction with the local spin density approximation (LSDA), improves the description of many properties, but overall, this improvement is limited. Here, we propose a modification to PZSIC that uses an iso-orbital indicator to identify regions where local SICs should be applied. Using this local-scaling SIC (LSIC) approach with LSDA, we analyze predictions for a wide range of properties including, for atoms, total energies, ionization potentials, and electron affinities and, for molecules, atomization energies, dissociation energy curves, reaction energies, and reaction barrier heights. LSIC preserves the results of PZSIC-LSDA for properties where it is successful and provides dramatic improvements for many of the other properties studied. Atomization energies calculated using LSIC are better than those of the Perdew, Burke, and Ernzerhof generalized gradient approximation (GGA) and close to those obtained with the strongly constrained and appropriately normed meta-GGA. LSIC also restores the uniform gas limit for the exchange energy that is lost in PZSIC-LSDA. Further performance improvements may be obtained by an appropriate combination or modification of the local scaling factor and the particular density functional approximation.

Structural Diversity in Supramolecular Organization of Anionic Phosphate Monoesters: Role of Cations.

Syntheses and structures of anionic arylphosphate monoesters [ArOP(O)2(OH)]- (Ar = 2,6-CHPh2-4-R-C6H2; R = Me/Et/iPr/tBu) with different counter cations are reported. The counter cations were varied systematically: imidazolium cation, 2-methyl imidazolium cation, N-methyl imidazolium cation, N,N'-alkyl substituted imidazolium cation, 1,4-diazabicyclo[2.2.2]octan-1-ium cation, 4,4'-bipyridinium dication, and magnesium(II) dication. The objective was to examine if the supramolecular structure of anionic arylphosphate monoesters could be modulated by varying the cation. It was found that an eight-membered P2O4H2-hydrogen-bonded dimeric motif involving intermolecular H-bonding between the [P(O)(OH)] unit of the anionic phosphate monoester along with the counter cation is formed with 2-methyl imidazolium cation, N-methyl imidazolium cation, N,N'-alkyl substituted imidazolium cation, 1,4-diazabicyclo[2.2.2]octan-1-ium cation, and magnesium(II) dication; both discrete and polymeric H-bonded structures are observed. In the case of imidazolium cations and 1,4-diazabicyclo[2.2.2]octan-1-ium cation, the formation of one-dimensional polymers (single lane/double lane) was observed. On the other hand, two types of phosphate motifs, intermolecular H-bonded dimer and an open-form, were observed in the case of 4,4'-bipyridinium dication.

Stretched or noded orbital densities and self-interaction correction in density functional theory.

Semilocal approximations to the density functional for the exchange-correlation energy of a many-electron system necessarily fail for lobed one-electron densities, including not only the familiar stretched densities but also the less familiar but closely related noded ones. The Perdew-Zunger (PZ) self-interaction correction (SIC) to a semilocal approximation makes that approximation exact for all one-electron ground- or excited-state densities and accurate for stretched bonds. When the minimization of the PZ total energy is made over real localized orbitals, the orbital densities can be noded, leading to energy errors in many-electron systems. Minimization over complex localized orbitals yields nodeless orbital densities, which reduce but typically do not eliminate the SIC errors of atomization energies. Other errors of PZ SIC remain, attributable to the loss of the exact constraints and appropriate norms that the semilocal approximations satisfy, suggesting the need for a generalized SIC. These conclusions are supported by calculations for one-electron densities and for many-electron molecules. While PZ SIC raises and improves the energy barriers of standard generalized gradient approximations (GGAs) and meta-GGAs, it reduces and often worsens the atomization energies of molecules. Thus, PZ SIC raises the energy more as the nodality of the valence localized orbitals increases from atoms to molecules to transition states. PZ SIC is applied here, in particular, to the strongly constrained and appropriately normed (SCAN) meta-GGA, for which the correlation part is already self-interaction-free. This property makes SCAN a natural first candidate for a generalized SIC.

Perdew-Zunger self-interaction correction: How wrong for uniform densities and large-Z atoms?

Semilocal density functionals for the exchange-correlation energy of a many-electron system cannot be exact for all one-electron densities. In 1981, Perdew and Zunger (PZ) subtracted the fully nonlocal self-interaction error orbital-by-orbital, making the corrected functional exact for all collections of separated one-electron densities and making no correction to the exact functional. Although the PZ self-interaction correction (SIC) eliminates many errors of semilocal functionals, it is often worse for equilibrium properties of sp-bonded molecules and solids. Nonempirical semilocal functionals are usually designed to be exact for electron gases of uniform density and, thus, also make 0% error for neutral atoms in the limit of large atomic number Z, but PZ SIC is not so designed. For localized SIC orbitals, we show analytically that the local spin density approximation (LSDA)-SIC correlation energy per electron of the uniform gas in the high-density limit makes an error of -50% in the spin-unpolarized case and -100% in the fully spin-polarized case. Then we extrapolate from the Ne, Ar, Kr, and Xe atoms to estimate the relative errors of the PZ SIC exchange-correlation energies (with localized SIC orbitals) in the limit of large atomic number: about +5.5% for the LSDA-SIC and about -3.5% for nonempirical generalized gradient [Perdew-Burke-Ernzerhof (PBE)-SIC] and meta-generalized gradient strongly constrained and appropriately normed (SCAN)-SIC approximations. The SIC errors are considerably larger than those that have been estimated for LSDA-SIC by approximating the localized SIC orbitals for the uniform gas and may explain the errors of PZ SIC for equilibrium properties, opening the door to a generalized SIC that is more widely accurate.

Modulation of the nuclearity of molecular Mg(ii)-phosphates: solid-state structural change involving coordinating solvents.

Herein, we report the synthesis and molecular structures of various magnesium(ii)-phosphate monoesters. By using a bulky aryl substituted phosphate monoester, ArOPO3H2 (Ar = 2,6-(CHPh2)2-4-tBu-C6H2), we have reproducibly assembled mono-, di-, tetra- (cage and ring), hexa-, and polynuclear magnesium(ii)-phosphate monoesters. Interestingly, the hexanuclear magnesium(ii)-phosphate monoester encapsulates an open-cage dodecanuclear water cluster.

Neutral and anionic phosphate-diesters as molecular templates for the encapsulation of a water dimer.

The study of water clusters in a confined state has received considerable attention. The simplest among such assemblies of water molecules is the water dimer. Stabilizing this unit at ambient conditions in an appropriate host-template, however, has remained a challenge. Herein, we report 2,6-(diphenylmethyl)-4-iso-propyl-phenyl substituted phosphate diesters (both in their neutral and anionic forms) as molecular templates for hosting the water dimer. The robustness of the water dimer within the neutral phosphate-diester ensemble is remarkable. Even when treated with triethylamine the water dimer is not disrupted, while exclusively the P(O)(OH) unit is deprotonated.

Root-growth of boron nitride nanotubes: experiments and ab initio simulations.

We have synthesized boron nitride nanotubes (BNNTs) in an arc in the presence of boron and nitrogen species. We find that BNNTs are often attached to large nanoparticles, suggesting that root-growth is a likely mechanism for their formation. Moreover, the tube-end nanoparticles are composed of boron, without transition metals, indicating that transition metals are not necessary for the arc synthesis of BNNTs. To gain further insight into this process we have studied key mechanisms for root growth of BNNTs on the surface of a liquid boron droplet by ab initio molecular dynamics simulations. We find that nitrogen atoms reside predominantly on the droplet surface where they organize to form boron nitride islands below 2400 K. To minimize contact with the liquid particle underneath, the islands assume non-planar configurations that are likely precursors for the thermal nucleation of cap structures. Once formed, the caps are stable and can easily incorporate nitrogen and boron atoms at their base, resulting in further growth. Our simulations support the root-growth mechanism of BNNTs and provide comprehensive evidence of the active role played by liquid boron.

Structural, electronic, and dynamical properties of liquid water by ab initio molecular dynamics based on SCAN functional within the canonical ensemble.

We perform ab initio molecular dynamics (AIMD) simulation of liquid water in the canonical ensemble at ambient conditions using the strongly constrained and appropriately normed (SCAN) meta-generalized-gradient approximation (GGA) functional approximation and carry out systematic comparisons with the results obtained from the GGA-level Perdew-Burke-Ernzerhof (PBE) functional and Tkatchenko-Scheffler van der Waals (vdW) dispersion correction inclusive PBE functional. We analyze various properties of liquid water including radial distribution functions, oxygen-oxygen-oxygen triplet angular distribution, tetrahedrality, hydrogen bonds, diffusion coefficients, ring statistics, density of states, band gaps, and dipole moments. We find that the SCAN functional is generally more accurate than the other two functionals for liquid water by not only capturing the intermediate-range vdW interactions but also mitigating the overly strong hydrogen bonds prescribed in PBE simulations. We also compare the results of SCAN-based AIMD simulations in the canonical and isothermal-isobaric ensembles. Our results suggest that SCAN provides a reliable description for most structural, electronic, and dynamical properties in liquid water.

Hydroxide diffuses slower than hydronium in water because its solvated structure inhibits correlated proton transfer.

Proton transfer via hydronium and hydroxide ions in water is ubiquitous. It underlies acid-base chemistry, certain enzyme reactions, and even infection by the flu. Despite two centuries of investigation, the mechanism underlying why hydroxide diffuses slower than hydronium in water is still not well understood. Herein, we employ state-of-the-art density-functional-theory-based molecular dynamics-with corrections for non-local van der Waals interactions, and self-interaction in the electronic ground state-to model water and hydrated water ions. At this level of theory, we show that structural diffusion of hydronium preserves the previously recognized concerted behaviour. However, by contrast, proton transfer via hydroxide is less temporally correlated, due to a stabilized hypercoordination solvation structure that discourages proton transfer. Specifically, the latter exhibits non-planar geometry, which agrees with neutron-scattering results. Asymmetry in the temporal correlation of proton transfer leads to hydroxide diffusing slower than hydronium.

Ab initio theory and modeling of water.

Water is of the utmost importance for life and technology. However, a genuinely predictive ab initio model of water has eluded scientists. We demonstrate that a fully ab initio approach, relying on the strongly constrained and appropriately normed (SCAN) density functional, provides such a description of water. SCAN accurately describes the balance among covalent bonds, hydrogen bonds, and van der Waals interactions that dictates the structure and dynamics of liquid water. Notably, SCAN captures the density difference between water and ice Ih at ambient conditions, as well as many important structural, electronic, and dynamic properties of liquid water. These successful predictions of the versatile SCAN functional open the gates to study complex processes in aqueous phase chemistry and the interactions of water with other materials in an efficient, accurate, and predictive, ab initio manner.

In situ Characterization of Nanoparticles Using Rayleigh Scattering.

We report a theoretical analysis showing that Rayleigh scattering could be used to monitor the growth of nanoparticles under arc discharge conditions. We compute the Rayleigh scattering cross sections of the nanoparticles by combining light scattering theory for gas-particle mixtures with calculations of the dynamic electronic polarizability of the nanoparticles. We find that the resolution of the Rayleigh scattering probe is adequate to detect nanoparticles as small as C60 at the expected concentrations of synthesis conditions in the arc periphery. Larger asymmetric nanoparticles would yield brighter signals, making possible to follow the evolution of the growing nanoparticle population from the evolution of the scattered intensity. Observable spectral features include characteristic resonant behaviour, shape-dependent depolarization ratio, and mass-dependent line shape. Direct observation of nanoparticles in the early stages of growth with unobtrusive laser probes should give insight on the particle formation mechanisms and may lead to better-controlled synthesis protocols.