Dispersion Correction Research Articles

Machine learning molecular dynamics simulations of liquid methanol

As the simplest hydrogen-bonded alcohol, liquid methanol has attracted intensive experimental and theoretical interest. However, theoretical investigations on this system have primarily relied on empirical intermolecular force fields or ab initio molecular dynamics with semilocal density functionals. Inspired by recent studies on bulk water using increasingly accurate machine learning force fields, we report a new machine learning force field for liquid methanol with a hybrid functional revPBE0 plus dispersion correction. Molecular dynamics simulations on this machine learning force field are orders of magnitude faster than ab initio molecular dynamics simulations, yielding the radial distribution functions, self-diffusion coefficients, and hydrogen bond network properties with very small statistical errors. The resulting structural and dynamical properties are compared well with the experimental data, demonstrating the superior accuracy of this machine learning force field. This work represents a successful step toward a first-principles description of this benchmark system and showcases the general applicability of the machine learning force field in studying liquid systems.

Multilevel quantum mechanical calculations show the role of promoter molecules in the dehydration of methanol to dimethyl ether in H-ZSM-5.

Methyl carboxylate esters promote the formation of dimethyl ether (DME) from the dehydration of methanol in H-ZSM-5 zeolite. We employ a multilevel quantum method to explore the possible associative and dissociative mechanisms in the presence, and absence, of six methyl ester promoters. This hybrid method combines density functional theory, with dispersion corrections (DFT-D3), for the full periodic system, with second-order Møller-Plesset perturbation theory (MP2) for small clusters representing the reaction site, and coupled cluster with single, double, and perturbative triple substitution (CCSD(T)) for the reacting molecules. The calculated adsorption enthalpy of methanol, and reaction enthalpies of the dehydration of methanol to DME within H-ZSM-5, agree with experiment to within chemical accuracy (∼4 kJ mol-1). For the promoters, a reaction pathway via an associative mechanism gives lower overall reaction enthalpies and barriers compared to the reaction with methanol only. Each stage of this mechanism is explored and related to experimental data. We provide evidence that suggests the promoter's adsorption to the Brønsted acid site is the most important factor dictating its efficiency.

Support effects on conical intersections of Jahn-Teller fluxional metal clusters on the sub-nanoscale.

The concept of fluxionality has been invoked to explain the enhanced catalytic properties of atomically precise metal clusters of subnanometer size. Cu3 isolated in the gas phase is a classical case of a fluxional metal cluster where a conical intersection leads to a Jahn-Teller (JT) distortion resulting in a potential energy landscape with close-lying multiminima and, ultimately, fluxional behavior. In spite of the role of conical intersections in the (photo)stability and (photo)catalytic properties of surface-supported atomic metal clusters, they have been largely unexplored. In this work, by applying a high-level multi-reference ab initio method aided with dispersion corrections, we analyze support effects on the conical intersection of Cu3 considering benzene as a model support molecule of carbon-based surfaces. We verify that the region around the conical intersection and the associated Jahn-Teller (JT) distortion is very slightly perturbed by the support when the Cu3 cluster approaches it in a parallel orientation: Two electronic states remain degenerate for a structure with C3 symmetry consistent with the D3h symmetry of unsupported Cu3 at the conical intersection. It extends over a one-dimensional seam that characterizes a physisorption minimum of the Cu3-benzene complex. The fluxionality of the Cu3 cluster, reflected in large fluctuations of relaxed Cu-Cu distances as a function of the active JT mode, is kept unperturbed upon complexation with benzene as well. In stark contrast, for the energetically favored perpendicular orientation of the Cu3 plane to the benzene ring plane, the conical intersection (CI) is located 12 100 cm-1 (∼1.5 eV) above the chemisorption minimum, with the fluxionality being kept at the CI's nearby and lost at the chemisorption well. The first excited state at the perpendicular orientation has a deep well (>4000 cm-1), being energetically closer to the CI. The transition dipole moment between ground and excited states has a significant magnitude, suggesting that the excited state can be observed through direct photo-excitation from the ground state. Besides demonstrating that the identity of an isolated Jahn-Teller metal cluster can be preserved against support effects at a physisorption state and lifted out at a chemisorption state, our results indicate that a correlation exists between conical intersection topography and fluxionality in the metal cluster's Cu-Cu motifs.

Element screening of metal sites in Fe-based Prussian blue framework materials for ammonium ion battery applications: a first-principles study.

Prussian blue framework materials are expected to be the next generation of electrode materials for commercial batteries because their three-dimensional framework structures facilitate the rapid transport and storage of ions and a variety of redox processes. This work compared the calculations of the model before and after the dispersion correction, and the model considering the effect of van der Waals force was more stable. In addition, the distances between H, C and N atoms were within the range of van der Waals force. Thus it was confirmed that NH4+ was adsorbed on the Ax site in the Prussian blue framework material (AxMa[Mb(CN)6]) by van der Waals interaction, and the charge transfer was mainly achieved by the interaction between the H atom in NH4+ and the N atom in the Prussian blue framework. On this basis, the properties of NH4+ batteries were theoretically screened for the Fe-based Prussian blue analogues (PBAs) with different Ma elements (Ma = Co, Cu, Fe, Mg, Mn, Ni, V or Zn). Considering the regulating effect of different metal elements on the electronic structures of PBAs, MgFe and ZnFe PBAs as the electrode materials of NH4+ batteries are expected to show excellent electrochemical energy storage performance in organic electrolytes.

Big data benchmarking: how do DFT methods across the rungs of Jacob's ladder perform for a dataset of 122k CCSD(T) total atomization energies?

Total atomization energies (TAEs) are a central quantity in density functional theory (DFT) benchmark studies. However, so far TAE databases obtained from experiment or high-level ab initio wavefunction theory included up to hundreds of TAEs. Here, we use the GDB-9 database of 133k CCSD(T) TAEs generated by Curtiss and co-workers [B. Narayanan, P. C. Redfern, R. S. Assary and L. A. Curtiss, Chem. Sci., 2019, 10, 7449] to evaluate the performance of 14 representative DFT methods across the rungs of Jacob's ladder (namely, PBE, BLYP, B97-D, M06-L, τ-HCTH, PBE0, B3LYP, B3PW91, ωB97X-D, τ-HCTHh, PW6B95, M06, M06-2X, and MN15). We first use the A25[PBE] diagnostic for nondynamical correlation to eliminate systems that potentially include significant multireference effects, for which the CCSD(T) TAEs might not be sufficiently reliable. The resulting database (denoted by GDB9-nonMR) includes 122k species. Of the considered functionals, B3LYP attains the best performance relative to the G4(MP2) reference TAEs, with a mean absolute deviation (MAD) of 4.09 kcal mol-1. This first-generation hybrid functional, in which the three mixing coefficients were fitted against a small set of TAEs, is one of the few functionals that are not systematically biased towards overestimating the G4(MP2) TAEs, as demonstrated by a mean-signed deviation (MSD) of 0.45 kcal mol-1. The relatively good performance of B3LYP is followed by the heavily parameterized M06-L meta-GGA functional, which attains a MAD of 6.24 kcal mol-1. The PW6B95, M06, M06-2X, and MN15 functionals tend to systematically overestimate the G4(MP2) TAEs and attain MADs ranging between 18.69 (M06) and 28.54 (MN15) kcal mol-1. However, PW6B95 and M06-2X exhibit particularly narrow error distributions. Thus, scaling their TAEs by an empirical scaling factor reduces their MADs to merely 3.38 (PW6B95) and 2.85 (M06-2X) kcal mol-1. Empirical dispersion corrections (e.g., D3 and D4) are attractive, and therefore, their inclusion worsens the performance of methods that systematically overestimate the TAEs.

Effect of density functional approximations on the calculated Jahn–Teller distortion in bis(terpyridine)manganese(III) and related compounds

ContextBis(terpyridine)manganese(III) exhibits Jahn–Teller distortion due to the inequivalent occupation of the degenerate eg orbitals of this high-spin d4 pseudo octahedral complex. Due to the spatially constrained nature of the terpyridine ligand, the central Mn-N bonds will always be shorter than the Mn-N terminal bonds, making it more difficult to distinguish between compression and elongation Jahn–Teller structures for bis(terpyridine)manganese(III). Density functional theory (DFT) calculations were utilized as a tool to evaluate the type of Jahn–Teller distortion in the high-spin d4 bis(terpyridine)manganese(III). The nature of the Jahn–Teller distortion calculated does depend upon the choice of density functional approximation (DFA) with the B3LYP, M06, and OLYP-D3 DFAs giving compression and the PW6B95D3, MN15, and MN15-D3 DFAs giving elongation in gas-phase calculations. All solvent-phase calculations yield an elongated structure for the bis(terpyridine)manganese(III) compound, which is yet to be structurally characterized experimentally. However, both gas and solvent OLYP-D3 calculations result in a compressed structure for the only experimentally isolated and characterized bis(terpyridine)manganese(III) complex, specifically the complex with terpyridine = 4′-(4-methylphenyl)-2,2′:6′,2′′-terpyridine. This alignment with the experimentally observed compression Jahn–Teller structure enhances the credibility of OLYP-D3 calculations in reproducing the observed geometries. The compressed Jahn–Teller geometries were near D2d symmetry with the z-axis for compression defined along the Mn-N central bonds. Elongation Jahn–Teller distortion is not possible along the Mn-N central bonds, due to their spatially constrained nature. Thus, elongation occur along one pair of opposite Mn-N terminal bonds that are longer than the other pair of opposite terminal bonds, with shorter central bonds. The highest symmetry of the elongation Jahn–Teller distortion geometry of bis(terpyridine)manganese(III) is C2v. Criteria to distinguish between a compression and elongation Jahn–Teller geometry for bis(terpyridine)manganese(III) are identified. The nature of the singly occupied eg molecular orbital, exhibiting anti-bonding interaction with the nitrogen-p MOs involved, dictates the type of Jahn–Teller distortion that occurs. The low-energy occupied bonding t2g molecular orbitals establish bonds with and undergo mixing with the ligand molecular orbitals. The OLYP-D3 functional is recommended for calculating bis(terpyridine)manganese(III) and related compounds due to its consistent generation of metal–ligand bonds slightly longer than observed in experiments, in line with the required behavior. Additionally, OLYP-D3 offers a realistic electronic structure for Jahn–Teller distorted bis(terpyridine)manganese(III), correctly identifying alpha eg molecular orbitals as the highest occupied molecular orbital and lowest unoccupied molecular orbital in agreement with experimental electrochemical studies. Furthermore, OLYP-D3 accurately reproduces the experimental compression geometry for the only structurally known bis(terpyridine)manganese(III) compound, instilling confidence in its reliability for such calculations.MethodsDFT geometry optimization and frequency calculations were done on the two different modes of Jahn–Teller distortion of bis(terpyridine)manganese(III), using the OLYP, B3LYP, M06, PW6B95D3, and MN15 functionals, with and without the Grimme’s D3 dispersion correction, and the 6-311G(d,p) or def2TZVPP basis set, as implemented in Gaussian 16. All optimizations were in the gas phase and also in the solvent phase with CH3CN as implicit solvent using IEFPCM.Graphical DFT calculations were utilized to determine the Jahn–Teller effect on the geometry of high-spin d4 bis(terpyridine)manganese(III) complex containing two structurally constrained tridentate ligands.

Intramolecular excimers of open forms of 2H-benzopyran, 2H- and 3H-naphthopyrans in solution: TD-DFT/DFT analysis.

We hypothesized that in the open forms of diphenyl-substituted photochromic compounds, immediately after the photoinduced cleavage of the C-O bond, one of the phenyl rings forms a stack with an aromatic system at the other end of the alkyl linker. The formation of these intramolecular stacking excimers is made possible by a known increase in the rotational mobility of the linker double bonds in the excited state. The flexibilities of linkers are analyzed in terms of changes in their dihedral angles. After internal conversion, stacking is partially broken but it still prevents the recovery of the C-O bond. This explains the significantly greater stability of the open forms of diphenyl-substituted benzopyran and naphthopyran compared to dimethyl-substituted ones. Open structures with residual stacking are intermediate and by thermal motion, they transform into stable trans-trans and cis-trans forms. However, in our opinion, they significantly affect the equilibrium of closed and open isomers in solution. In particular, the calculated stacking energy of 2H-naphthopyran turned out to be higher than that of 3H-naphthopyran, which partly explains the significantly greater stability of the open form of the former. For the DFT/TD-DFT calculations performed in this work, the following functionals were selected to provide reliable stacking for the open forms of all three studied compounds: APFD, M052X, M06HF, as well as B3LYP with D3(BJ) Grimme dispersion correction. The 6-311++G(d,p) basis set and cyclohexane solvent modeled using IEFPCM were also used.

Design and synthesis of a new coordination polymer of zinc (II): Crystal structure and evaluation for its adsorption studies in removing perilous organic dye pollutants in aqueous medium

The contamination of aqueous environments by hazardous organic dye pollutants has become a pressing global concern, necessitating innovative strategies for effective water purification. In this research, we present the design, synthesis, and comprehensive evaluation of a new zinc(II) coordination polymer, formulated as [Zn(HIPA)(HIZ)2]n⋅H2O (CP-1), specifically tailored for adsorption studies aimed at the removal of hazardous organic dye pollutants from aqueous mediums. The CP-1 was effectively synthesized using a hydrothermal synthesis approach and incorporates two distinct organic ligands, namely, 5‑hydroxy isophthalic acid, and imidazole, along with a water molecule within the crystal lattice. The structural and morphological attributes of the synthesized coordination polymer (CP-1) are elucidated through a comprehensive characterization process, including SCXRD, PXRD, FTIR, UV–VIS spectroscopy, TGA, SEM, and BET analysis. Single-crystal diffraction of CP-1 confirms that Zn is coordinated in a tetrahedral environment. The CP-1 was systematically examined for dye adsorption of four hazardous organic dyes namely methylene blue (MB), methyl orange (MO), crystal violet (CV) and congo red (CR). Batch adsorption experiments were conducted, probing the influence of variables such as contact time, pH, temperature, the effect of dosage and initial dye concentration on the adsorption efficiency of CP-1. The CP-1 shows the best adsorption efficiency for MB followed by MO, CV, and CR. The CP-1 exhibits high selectivity and efficient MB and MO dyes adsorption at room temperature and neutral pH. Moreover, three different kinetic models were utilized to attain a deeper comprehension of the reaction mechanism involved in the adsorption of MB and MO. Amongst these models, the pseudo-second-order kinetic model demonstrated the best fit for both cases. Furthermore, the density functional theory with the dispersion correction (Grimme's DFT-D3) was utilized to support the dye's adsorption mechanism with the CP-1. To study this, the dispersion corrected interaction energies were computed for the dye-complex systems at B3LYP-D3/6–31G(d,p)/LANL2DZ level of theory. The efficient adsorption performance, coupled with its structural and chemical stability, suggests that CP-1 holds great potential for practical applications in water treatment processes aimed at addressing the pressing challenge of organic dyes pollutant contamination. Further research and exploration of its adsorption properties can unveil its full potential for future applications as a dye adsorbent.

Influence of Steric and Dispersion Interactions on the Thermochemistry of Crowded (Fluoro)alkyl Compounds.

ConspectusAlkanes play a pivotal role in industrial, environmental, and biological processes. They are characterized by their carbon-carbon single-bond structure, remarkable stability, and conformational diversity. Fluorination of such compounds imparts unique physicochemical properties that often enhance pharmacokinetic profiles, metabolic stability, and receptor interactions while keeping beneficial properties. However, such per- and polyfluoroalkyl substances (PFAS) show a persistent presence in the environment and potential adverse health effects, which propelled them to the forefront of global environmental and health discussions. Alkyl compounds are also prototypical for stereoelectronic (SE) effects that are widely applied in chemistry. Substituents are typically described as electron-density-donating/withdrawing and/or responsible for sterically interacting with reagents or strategic groups in the molecule. That alkane branching can result in higher stability compared to less-branched isomers has been investigated in detail also by testing quantum chemical methods, in particular density functional theory (DFT). Alkane branching results in spatially compact structures with close intramolecular contacts so that at a specific size the detailed balance of attractive London dispersion and covalent versus repulsive Pauli exchange interactions shifts to new, chemically unfragile situations. This may lead to dissociation at room temperature and opens the central question: what is the smallest crowed alkane that cannot be made synthetically? In this Account, we try to shed light on the interplay among the various (free) energy components for crowded (fluoro)alkane dissociation. In this context, homolytic cleavage of the central C-C bond in a series of model alkanes of increasing size with tert-butyl (tBu), adamantyl (Ad), and [1.1.1]propellanyl (Prop) substituents is investigated. Reference energies are calculated at the PNO-LCCSD(T)-F12b level and used to benchmark the performance of contemporary DFT functionals. In line with previous conclusions, the application of dispersion corrections to density functionals is mandatory. For crowed structures, the accurate description of the midrange correlation effects, specifically repulsive van der Waals interactions, is crucial, and we observed that the density-dependent VV10 correction is superior to D4 in this context, although the asymptotic region is better described by the latter. The best available dispersion-inclusive functionals show systematic and reasonably small residual errors and can be safely applied to large systems (>100 atoms), for which coupled cluster methods with large basis sets are not computationally feasible anymore. For qualitatively correct predictions of synthetic accessibility under equilibrium conditions (free energy), the inclusion of thermostatistical (entropy) contributions is also essential. According to our results, tetra-tert-butylmethane (C17tBu) is the largest and most crowded system with a positive dissociation free energy and should be synthesizable. The difference between hydrogenated and perfluorinated systems originates from the increase in the steric repulsion of spatially close substituents, which is not compensated to the same extent by attractive orbital and dispersion interactions. A sometimes-assumed similar steric demand for fluorine and hydrogen atoms is not corroborated by our investigations on crowded systems. Perfluorination is found to substantially decrease thermal stability, rendering perfluorinated hexamethylethane (C8tBuF) the last potentially stable representative.

Combining the D3 dispersion correction with the neuroevolution machine-learned potential

Machine-learned potentials (MLPs) have become a popular approach of modeling interatomic interactions in atomistic simulations, but to keep the computational cost under control, a relatively short cutoff must be imposed, which put serious restrictions on the capability of the MLPs for modeling relatively long-ranged dispersion interactions. In this paper, we propose to combine the neuroevolution potential (NEP) with the popular D3 correction to achieve a unified NEP-D3 model that can simultaneously model relatively short-ranged bonded interactions and relatively long-ranged dispersion interactions. We show that improved descriptions of the binding and sliding energies in bilayer graphene can be obtained by the NEP-D3 approach compared to the pure NEP approach. We implement the D3 part into the gpumd package such that it can be used out of the box for many exchange-correlation functionals. As a realistic application, we show that dispersion interactions result in approximately a 10% reduction in thermal conductivity for three typical metal-organic frameworks.

Structures and stability of K+ cation solvated in Arn clusters

The solvation of K+ cation plays an important role in various phenomena such as biological procedures, geological time, and archaeological properties. Monte Carlo (MC) simulation and DFT method are employed to study the structural and energetic characteristics of the K + Arn (n = 1–14) clusters. The potential model (PM) and the Basin-Hopping (BH) method are the foundation of the MC simulation. The pairwise PM (PW-PM) is improved by introducing the N-body interactions via the polarizable potential model (PPM). On the other side, the DFT functional M05–2X, combined with the 6–311++G(3d2f,2p) basis set, and the Grimme dispersion correction GD3 was used to deeply investigate the geometrical properties and the relative stability of the K + Arn clusters. Starting from n = 12, a structural transition from square antiprism (SA) to icosahedron (ICOS) form is detected. Additionally, the PPM allows us to examine the largest sizes (n = 15–54). Herein, the first ICOS layers are found for n = 12 and 54 cluster sizes, respectively. The binding energy and the second energy difference as a function of cluster size are used to evaluate the relative stability of K + Arn clusters. The obtained data are in concordance with the available results in the literature.

Strain-induced effects on the optoelectronic properties of ZrSe2/HfSe2 heterostructures.

Two-dimensional semiconductor materials have received much attention in recent years due to their wide variety of applications in the field of nano-optoelectronic devices. In this project, we applied stresses ranging from -6 to +6% to the ZrSe2/HfSe2 heterostructure and systematically investigated its electrical and optical properties. It is discovered that stress can effectively modulate the forbidden bandwidth of the ZrSe2/HfSe2 heterojunction; whereas, under compressive stress, the forbidden bandwidth of the material decreases further until the bandgap is zero, leading to the material's transformation from semiconductor to metal. The forbidden band gap of the ZrSe2/HfSe2 heterojunction increases with increasing horizontal biaxial tensile strain. We discovered that the light absorption performance of this heterostructure is significantly better than that of its similar monomolecular layer and that its light absorption intensity can reach an order of magnitude of 104. Under compressive and tensile stresses, the ZrSe2/HfSe2 heterojunctions exhibit different degrees of red or blue shift. The results indicate that constructing ZrSe2/HfSe2 heterojunctions and applying horizontal biaxial stresses to them can significantly modulate the optoelectronic properties of the materials. ZrSe2/HfSe2 heterojunction is a new type of high-performance photogenerated carrier transport device with a wide range of applications. The calculations in this study are carried out the first principles approach of density functional theory, as implemented in the CASTEP module of Materials-Studio2019. The researchers used an ultrasoft reaction potential to calculate the interactions between the ion core and the electrons and applied the Perdew-Burke-Ernzerhof (PBE) and the generalized gradient approximation (GGA) to perform the calculations. The Monkhorst-Pack technique was employed to create the k-point samples utilized for integration on the Brillouin zone, and the k-point grid was uniformly 6 × 6 × 1. In addition, in order to avoid interactions between the atomic layers affecting the properties and stability of the material, such interactions were prevented by adding a 30 Å vacuum layer. Using a plane-wave energy cutoff of 500 eV and the convergence accuracy of the iterative process was set to 1 × 10-5 eV to ensure the accuracy of the computational results, and in addition. The maximum stress in the lattice was limited to less than 0.05 GPa or the interaction force between neighboring atoms was lower than 0.03 eV/Å. For the calculation of the properties of the optical properties, a k-point grid of 18 × 18 × 1 is used for optimization, and the polarization direction of the material is not taken into account, considering that the material is isotropic. This study proposes to apply the Tkatchenko-Scheffler (TS) dispersion correction method in DFT-D to appropriately represent the interlayer van der Waals interaction forces to solve inaccuracies in the computation of van der Waals interactions via density functional theory.

Dispersion-corrected r2SCAN based double-hybrid functionals

The regularized and restored semi-local meta-generalized gradient approximation (meta-GGA) exchange–correlation functional r2SCAN [Furness et al., J. Phys. Chem. Lett. 11, 8208–8215 (2020)] is used to create adiabatic-connection-derived global double-hybrid functionals employing spin-opposite-scaled MP2. The 0-DH, CIDH, QIDH, and 0–2 type double-hybrid functionals are assessed as a starting point for further modification. Variants with 50% and 69% Hartree–Fock exchange (HFX) are empirically optimized (Pr2SCAN50 and Pr2SCAN69), and the effect of MP2-regularization (κPr2SCAN50) and range-separated HFX (ωPr2SCAN50) is evaluated. All optimized functionals are combined with the state-of-the-art London dispersion corrections D4 and NL. The resulting functionals are assessed comprehensively for their performance on main-group and metal-organic thermochemistry on 90 different benchmark sets containing 25 800 data points. These include the extensive GMTKN55 database, additional sets for main-group chemistry, and multiple sets for transition-metal complexes, including the ROST61, the MOR41, and the MOBH35 sets. As the main target of this study is the development of a broadly applicable, robust functional with low empiricism, special focus is put on variants with moderate amounts of HFX (50%), which are compared to the so far successful PWPB95-D4 (50% HFX, 20% MP2 correlation) functional. The overall best variant, ωPr2SCAN50-D4, performs well on main-group and metal-organic thermochemistry, followed by Pr2SCAN69-D4 that offers a slight edge for metal-organic thermochemistry and by the low HFX global double-hybrid Pr2SCAN50-D4 that performs robustly across all tested sets. All four optimized functionals, Pr2SCAN69-D4, Pr2SCAN50-D4, κPr2SCAN50-D4, and ωPr2SCAN50-D4, outperform the PWPB95-D4 functional.

Investigating the molecular crystals of L-Alanine, DL-Alanine, β-Alanine, and Alanine hydrogen chloride: Experimental and DFT analysis of structural and optoelectronic properties

In this study, DFT calculations were employed to assess and contrast the structural, electronic, and optical characteristics of five alanine crystal forms. The Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) with the Tkatchenko-Scheffler dispersion correction scheme was employed and compared to the PBE functional and the local density approximation (LDA) without dispersion correction. Experimentally, alanine crystal powders were characterized through X-ray diffraction (XRD) and UV–VIS absorption spectrum. Additionally, UV–VIS absorption measurements were performed on L-α-alanine, DL-alanine, β-alanine, and alanine hydrochloride solvated in water. Significant deviations (>10%) were observed between the experimental crystal lattice parameters and those predicted by the GGA and LDA methods. Incorporating dispersion correction into the GGA functional improved prediction accuracy. Additionally, time-dependent DFT (TD-DFT) computations using the polarizable continuum model were conducted to analyze the UV/VIS absorption spectrum of α-, β-, DL-, and hydrogen chloride alanine in water and interpret experimental data. Different alanine crystals exhibited substantial anisotropic electronic and optical properties, including effective masses, absorption, and dielectric function. The L-α-HCl structure displayed the smallest band gap (4.67 eV), suggesting its suitability for applications such as L-α/L-α-HCl/L-α quantum well structures. The main effective masses of electrons and holes were also determined for each crystal for a better understanding of their electronic transport properties.

Analytical estimation of sea-level rise impacts on the freshwater lenses of elliptical islands with sloping shorelines

The vulnerability of island freshwater lenses to rising sea levels has been assessed previously using sharp-interface analytical solutions that consider islands as circular or two-dimensional strips, or the lenses are treated as entirely underlain by seawater. However, many islands do not fit these basic conditions. The current study develops an analytical solution for the freshwater-saltwater interface in elliptical islands where the interface intercepts the aquifer base, and explores (for the first time for island lenses) the applicability of a dispersion correction to the analytical solution. An analysis of the lens response to sea-level rise (SLR) and the accompanying land-surface inundation (LSI) shows that for lenses entirely underlain by seawater, the primary impact of SLR is likely the associated LSI. Counterintuitively, lenses in contact with the aquifer basement may experience an increase in the freshwater volume with SLR, although this requires that LSI is limited (e.g., the shoreline is steep) and the water table is sufficiently deep (below the land surface) to allow it to rise commensurately with SLR. The analysis also shows that as the aspect ratio of the elliptical island increases (i.e., larger values approach strip islands and circular islands have a value of 1), the lens volume decreases (for a given island area). Thus, circular islands have the maximum freshwater storage per island area. Specifically, for typical island aquifer conditions, freshwater storage decreased over 50% as a/b increased from 1 to 10. The results obtained from this study and the new analytical solution are expected to assist in the rapid assessment of the freshwater resources of elliptical-like islands.

Benchmark Data Set of Crystalline Organic Semiconductors.

This work reports a Benchmark Data set of Crystalline Organic Semiconductors to test calculations of the structural and electronic properties of these materials in the solid state. The data set contains 67 crystals consisting of mostly rigid molecules with a single dominant conformer, covering the majority of known structural types. The experimental crystal structure is available for the entire data set, whereas zero-temperature unit cell volume can be reliably estimated for a subset of 28 crystals. Using this subset, we benchmark r2SCAN-D3 and PBE-D3 density functionals. Then, for the entire data set, we benchmark approximate density functional theory (DFT) methods, including GFN1-xTB and DFTB3(3ob-3-1), with various dispersion corrections against r2SCAN-D3. Our results show that r2SCAN-D3 geometries are accurate within a few percent, which is comparable to the statistical uncertainty of experimental data at a fixed temperature, but the unit cell volume is systematically underestimated by 2% on average. The several times faster PBE-D3 provides an unbiased estimate of the volume for all systems except for molecules with highly polar bonds, for which the volume is substantially overestimated in correlation with the underestimation of atomic charges. Considered approximate DFT methods are orders of magnitude faster and provide qualitatively correct but overcompressed crystal structures unless the dispersion corrections are fitted by unit cell volume.

Improved Description of Intra- and Intermolecular Interactions through Dispersion-Corrected Second-Order Møller-Plesset Perturbation Theory.

ConspectusThe quantum chemical modeling of organic crystals and other molecular condensed-phase problems requires computationally affordable electronic structure methods which can simultaneously describe intramolecular conformational energies and intermolecular interactions accurately. To achieve this, we have developed a spin-component-scaled, dispersion-corrected second-order Møller-Plesset perturbation theory (SCS-MP2D) model. SCS-MP2D augments canonical MP2 with a dispersion correction which removes the uncoupled Hartree-Fock dispersion energy present in canonical MP2 and replaces it with a more reliable coupled Kohn-Sham treatment, all evaluated within the framework of Grimme's D3 dispersion model. The spin-component scaling is then used to improve the description of the residual (nondispersion) portion of the correlation energy.The SCS-MP2D model improves upon earlier corrected MP2 models in a few ways. Compared to the highly successful dispersion-corrected MP2C model, which is based solely on intermolecular perturbation theory, the SCS-MP2D dispersion correction improves the description of both inter- and intramolecular interactions. The dispersion correction can also be evaluated with trivial computational cost, and nuclear analytic gradients are computed readily to enable geometry optimizations. In contrast to earlier spin-component scaling MP2 models, the optimal spin-component scaling coefficients are only mildly sensitive to the choice of training data, and a single global parametrization of the model can describe both thermochemistry and noncovalent interactions.The resulting dispersion-corrected, spin-component-scaled MP2 (SCS-MP2D) model predicts conformational energies and intermolecular interactions with accuracy comparable to or better than those of many range-separated and double-hybrid density functionals, as is demonstrated on a variety of benchmark tests. Among the functionals considered here, only the revDSD-PBEP86-D3(BJ) functional gives consistently smaller errors in benchmark tests. The results presented also hint that further improvements of SCS-MP2D may be possible through a more robust fitting procedure for the seven empirical parameters.To demonstrate the performance of SCS-MP2D further, several applications to molecular crystal problems are presented. The three chosen examples all represent cases where density-driven delocalization error causes GGA or hybrid density functionals to artificially stabilize crystals exhibiting more extended π-conjugation. Our pragmatic strategy addresses the delocalization error by combining a periodic density functional theory (DFT) treatment of the infinite lattice with intramolecular/conformational energy corrections computed with SCS-MP2D. For the anticancer drug axitinib, applying the SCS-MP2D conformational energy correction produces crystal polymorph stabilities that are consistent with experiment, in contrast to earlier studies. For the crystal structure prediction of the ROY molecule, so named for its colorful red, orange, and yellow crystals, this approach leads to the first plausible crystal energy landscape, and it reveals that the lowest-energy polymorphs have already been found experimentally. Finally, in the context of photomechanical crystals, which transform light into mechanical work, these techniques are used to predict the structural transformations and extract design principles for maximizing the work performed.

Do Optimally Tuned Range-Separated Hybrid Functionals Require a Reparametrization of the Dispersion Correction? It Depends.

For ground- and excited-state studies of large molecules, it is the state of the art to combine (time-dependent) DFT with dispersion-corrected range-separated hybrid functionals (RSHs), which ensures an asymptotically correct description of exchange effects and London dispersion. Specifically for studying excited states, it is common practice to tune the range-separation parameter ω (optimal tuning), which can further improve the accuracy. However, since optimal tuning essentially changes the functional, it is unclear if and how much the parameters used for the dispersion correction depend on the chosen ω value. To answer this question, we explore this interdependency by refitting the DFT-D4 dispersion model for six established RSHs over a wide range of ω values (0.05-0.45 a0-1) using a set of noncovalently bound molecular complexes. The results reveal some surprising differences among the investigated functionals: While PBE-based RSHs and ωB97M-D4 generally exhibit a weak interdependency and robust performance over a wide range of ω values, B88-based RSHs, specifically LC-BLYP, are strongly affected. For these, even a minor reduction of ω from the default value manifests in strong systematic overbinding and poor performance in the typical range of optimally tuned ω values. Finally, we discuss strategies to mitigate these issues and reflect the results in the context of the employed D4 parameter optimization algorithm and fit set, outlining strategies for future improvements.

Evaluation of impairment mitigations for optical fiber communications using dispersion compensation techniques

This study discusses dispersion as a key issue influencing the performance of optical fiber communication technology. Poor bit rate, pulse broadening, and transmission distance restrictions are the outcomes. The dispersion correction method frequently employs fiber with dispersion compensation (DCF) and fiber Bragg grating (FBG). A model transmission system has been analyzed in light of numerous parameters using the Opti System 20.0 simulator. The DCF and FBG dispersion compensation methods are two that have been proposed. The Q-factor, BER, and eye height measured at the receiver end will be utilized as three metrics to describe the results of the two dispersion-correcting strategies. Additionally, the system's transmission performance will vary as a result of varied compensatory approaches. The DCF approach does not significantly increase system performance compared to the compensation technique. In addition, the properties of the optical fiber used for transmission in the system affect the characteristics of the FBG and the fiber used for the compensating process. It has been demonstrated that post-compensation, followed by symmetrical compensation, produces the best result. In addition, the transmission performance of the system with different compensation methods will be different, and the compensation technique will improve the system performance much better than the DCF technique. Moreover, the characteristics of the FBG and the fiber that is used for the compensation process depend on the characteristics of the optical fiber that is used for the transmission in system. It is shown that the symmetrical compensation performance is the best. Due to the nonlinear nature of propagation, system performance depends on power levels, followed by post-compensation.

Low thermal expansion of layered electrides predicted by density-functional theory.

Layered electrides are a unique class of materials with anionic electrons bound in interstitial regions between thin, positively charged atomic layers. While density-functional theory is the tool of choice for computational study of electrides, there has to date been no systematic comparison of density functionals or dispersion corrections for their accurate simulation. There has also been no research into the thermomechanical properties of layered electrides, with computational predictions considering only static lattices. In this work, we investigate the thermomechanical properties of five layered electrides using density-functional theory to evaluate the magnitude of thermal effects on their lattice constants and cell volumes. We also assess the accuracy of five popular dispersion corrections with both planewave and numerical atomic orbital calculations.