Interaction between water and carbon nanostructures: How good are current density functional approximations?

Due to their current and future technological applications, including realisation of water filters and desalination membranes, water adsorption on graphitic sp$^{2}$-bonded carbon is of overwhelming interest. However, these systems are notoriously challenging to model, even for electronic structure methods such as density functional theory (DFT), because of the crucial role played by London dispersion forces and non-covalent interactions in general. Recent efforts have established reference quality interactions of several carbon nanostructures interacting with water. Here, we compile a new benchmark set (dubbed \textbf{WaC18}), which includes a single water molecule interacting with a broad range of carbon structures, and various bulk (3D) and two dimensional (2D) ice polymorphs. The performance of 28 approaches, including semi-local exchange-correlation functionals, non-local (Fock) exchange contributions, and long-range van der Waals (vdW) treatments, are tested by computing the deviations from the reference interaction energies. The calculated mean absolute deviations on the WaC18 set depends crucially on the DFT approach, ranging from 135 meV for LDA to 12 meV for PBE0-D4. We find that modern vdW corrections to DFT significantly improve over their precursors. Within the 28 tested approaches, we identify the best performing within the functional classes of: generalized gradient approximated (GGA), meta-GGA, vdW-DF, and hybrid DF, which are BLYP-D4, TPSS-D4, rev-vdW-DF2, and PBE0-D4, respectively

Evidence for stable square ice from quantum Monte Carlo

Recent experiments on ice formed by water under nanoconfinement provide evidence for a two-dimensional (2D) “square ice” phase. However, the interpretation of the experiments has been questioned and the stability of square ice has become a matter of debate. Partially this is because the simulation approaches employed so far (force fields and density functional theory) struggle to accurately describe the very small energy differences between the relevant phases. Here we report a study of 2D ice using an accurate wave-function based electronic structure approach, namely diffusion Monte Carlo (DMC). We find that at relatively high pressure, square ice is indeed the lowest enthalpy phase examined, supporting the initial experimental claim. Moreover, at lower pressures, a “pentagonal ice” phase (not yet observed experimentally) has the lowest enthalpy, and at ambient pressure, the “pentagonal ice” phase is degenerate with a “hexagonal ice” phase. Our DMC results also allow us to evaluate the accuracy of various density functional theory exchange-correlation functionals and force field models, and in doing so we extend the understanding of how such methodologies perform to challenging 2D structures presenting dangling hydrogen bonds

Correction to Thermochemistry of Racemic and Enantiopure Organic Crystals for Predicting Enantiomer Separation

Correction to Thermochemistry of Racemic and Enantiopure Organic Crystals for Predicting Enantiomer Separatio

Interactions between large molecules pose a puzzle for reference quantum mechanical methods

Quantum-mechanical methods are used for understanding molecular interactions throughout the natural sciences. Quantum diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] are state-of-the-art trusted wavefunction methods that have been shown to yield accurate interaction energies for small organic molecules. These methods provide valuable reference information for widely-used semi-empirical and machine learning potentials, especially where experimental information is scarce. However, agreement for systems beyond small molecules is a crucial remaining milestone for cementing the benchmark accuracy of these methods. We show that CCSD(T) and DMC interaction energies are not consistent for a set of polarizable supramolecules. Whilst there is agreement for some of the complexes, in a few key systems disagreements of up to 8 kcal mol−1 remain. These findings thus indicate that more caution is required when aiming at reproducible non-covalent interactions between extended molecules

Crystal structure prediction of flexible pharmaceutical-like molecules: Density functional tight-binding as an intermediate optimization method and for free energy estimation

Successful methodologies for theoretical crystal structure prediction (CSP) on flexible pharmaceutical-like organic molecules explore the lattice energy surface to find a set of plausible crystal structures. The initial search stage of CSP studies uses a relatively simple lattice energy approximation as hundreds of thousands of minima have to be considered. These generated crystal structures often have poor molecular geometries, as well as inaccurate lattice-energy rankings, and performing reasonably accurate but computationally affordable optimisations of the crystal structures generated in a search would be highly desirable. Here, we seek to explore whether semi-empirical quantum-mechanical methods can perform this task. We employed the dispersion-corrected tight-binding Hamiltonian (DFTB3-D3) to relax all inter and intra-molecular degrees of freedom of several thousands of generated crystal structures of five pharmaceutical-like molecules, saving a large amount of computational effort compared to earlier studies. The computational cost scales better with molecular size and flexibility than other CSP methods, suggesting it could be extended to even larger and more flexible molecules. On average, this optimisation improved the average reproduction of the eight experimental crystal structures (RMSD15 ) and experimental conformers (RMSD1) by 4% and 23%, respectively. The intermolecular interactions were then further optimised using distributed multipoles, derived from the molecular wave-function, to accurately describe the electrostatic component of the intermolecular energy. In all cases, the experimental crystal structures are close to the top of the lattice energy ranking. Phonon calculations on some of the lowest energy structures were also performed with DFTB3-D3 methods to calculate the vibrational component of the Helmholtz free energy, providing further insights into the solid-state behaviour of the target molecules. We conclude that DFTB3-D3 is a cost-effective method for optimising flexible molecules, bridging the gap between the approximate methods used in CSP searches for generating crystal structures and more accurate methods required in the final energy ranking

Measurement of the Relative Branching Fraction of Υ(4S)\Upsilon(4S) to Charged and Neutral B-Meson Pairs

We analyze 9.7 x 10^6 B\bar{B}$ pairs recorded with the CLEO detector to determine the production ratio of charged to neutral B-meson pairs produced at the Y(4S) resonance. We measure the rates for B^0 -> J/psi K^{(*)0} and B^+ -> J/psi K^{(*)+} decays and use the world-average B-meson lifetime ratio to extract the relative widths f+-/f00 = Gamma(Y(4S) -> B+B-)/Gamma(Y(4S) -> B0\bar{B0}) = = 1.04 +/- 0.07(stat) +/- 0.04(syst). With the assumption that f+- + f00 = 1, we obtain f00 = 0.49 +/- 0.02(stat) +/- 0.01(syst) and f+- = 0.51 +/- 0.02(stat) +/- 0.01(syst). This production ratio and its uncertainty apply to all exclusive B-meson branching fractions measured at the Y(4S) resonance.Comment: 11 pages postscript, also available through http://w4.lns.cornell.edu/public/CLN

First Observation of the Decays B0→D∗−ppˉπ+B^{0}\to D^{*-}p\bar{p}\pi^{+} and B^{0}\to D^{*-}p\bar{n}$

We report the first observation of exclusive decays of the type B to D^* N anti-N X, where N is a nucleon. Using a sample of 9.7 times 10^{6} B-Bbar pairs collected with the CLEO detector operating at the Cornell Electron Storage Ring, we measure the branching fractions B(B^0 \to D^{*-} proton antiproton \pi^+) = ({6.5}^{+1.3}_{-1.2} +- 1.0) \times 10^{-4} and B(B^0 \to D^{*-} proton antineutron) = ({14.5}^{+3.4}_{-3.0} +- 2.7) times 10^{-4}. Antineutrons are identified by their annihilation in the CsI electromagnetic calorimeter.Comment: 9 pages postscript, also available through http://w4.lns.cornell.edu/public/CLN

Study of the Decays B0 --> D(*)+D(*)-

The decays B0 --> D*+D*-, B0 --> D*+D- and B0 --> D+D- are studied in 9.7 million Y(4S) --> BBbar decays accumulated with the CLEO detector. We determine Br(B0 --> D*+D*-) = (9.9+4.2-3.3+-1.2)e-4 and limit Br(B0 --> D*+D-) < 6.3e-4 and Br(B0 --> D+D-) < 9.4e-4 at 90% confidence level (CL). We also perform the first angular analysis of the B0 --> D*+D*- decay and determine that the CP-even fraction of the final state is greater than 0.11 at 90% CL. Future measurements of the time dependence of these decays may be useful for the investigation of CP violation in neutral B meson decays.Comment: 21 pages, 5 figures, submitted to Phys. Rev.

A Search for B→τνB\to \tau\nu

We report results of a search for $B\to\tau\nu$ in a sample of 9.7 million charged $B$ meson decays. The search uses both $\pi\nu$ and $\ell\nu\bar\nu$ decay modes of the $\tau$, and demands exclusive reconstruction of the companion $\bar B$ decay to suppress background. We set an upper limit on the branching fraction ${\cal B}(B\to \tau\nu) < 8.4\times 10^{-4}$ at 90% confidence level. With slight modification to the analysis we also establish ${\cal B}(B^\pm\to K^\pm\nu\bar\nu) < 2.4\times 10^{-4}$ at 90% confidence level.Comment: 10 ages postscript, also available through http://w4.lns.cornell.edu/public/CLN

Measurements of B --> D_s^{(*)+} D^{*(*)} Branching Fractions

This article describes improved measurements by CLEO of the $B^0 \to D_s^+ D^{*-}$ and $B^0 \to D_s^{*+} D^{*-}$ branching fractions, and first evidence for the decay $B^+ \to D_s^{(*)+} \bar{D}^{**0}$, where $\bar{D}^{**0}$ represents the sum of the $\bar{D}_1(2420)^0$, $\bar{D}_2^*(2460)^0$, and $\bar{D}_1(j=1/2)^0$ L=1 charm meson states. Also reported is the first measurement of the $D_s^{*+}$ polarization in the decay $B^0 \to D_s^{*+} D^{*-}$. A partial reconstruction technique, employing only the fully reconstructed $D_s^+$ and slow pion $\pi_s^-$ from the $D^{*-} \to \bar{D}^0 \pi^-_s$ decay, enhances sensitivity. The observed branching fractions are ${\mathcal B} (B^0 \to D_s^+ D^{*-}) = (1.10 \pm 0.18 \pm 0.10 \pm 0.28)$, ${\mathcal B} (B^0 \to D_s^{*+} D^{*-}) = (1.82 \pm 0.37 \pm 0.24 \pm 0.46)$, and ${\mathcal B} (B^+ \to D_s^{(*)+} \bar{D}^{**0}) = (2.73 \pm 0.78 \pm 0.48 \pm 0.68)$, where the first error is statistical, the second systematic, and the third is due to the uncertainty in the $D_s^+ \to \phi \pi^+$ branching fraction. The measured $D_s^{*+}$ longitudinal polarization, $\Gamma_L/\Gamma = (50.6 \pm 13.9 \pm 3.6)$, is consistent with the factorization prediction of 54%.Comment: 26 pages (LaTeX), 15 figures. To be submitted to PR