Noncovalent Interactions by QMC: Speedup by One-Particle Basis-Set Size Reduction

While it is empirically accepted that the fixed-node diffusion Monte-Carlo (FN-DMC) depends only weakly on the size of the one-particle basis sets used to expand its guiding functions, limits of this observation are not settled yet. Our recent work indicates that under the FN error cancellation conditions, augmented triple zeta basis sets are sufficient to achieve a benchmark level of 0.1 kcal/mol in a number of small noncovalent complexes. Here we report on a possibility of truncation of the one-particle basis sets used in FN-DMC guiding functions that has no visible effect on the accuracy of the production FN-DMC energy differences. The proposed scheme leads to no significant increase in the local energy variance, indicating that the total CPU cost of large-scale benchmark noncovalent interaction energy FN-DMC calculations may be reduced.Comment: ACS book chapter, accepte

Pink‐ and orange‐pigmented Planctomycetes produce saproxanthin‐type carotenoids including a rare C45 carotenoid

Planctomycetes, are ubiquitous and environmentally important Gram-negative aquatic bacteria with key roles in global carbon and nitrogen cycles. Many planctomycetal species have a pink or orange colour and have been suggested to produce carotenoids. Potential applications as food colorants or anti-oxidants have been proposed. Hitherto, the planctomycetal metabolism is largely unexplored and the strain pigmentation has not been identified. For a holistic view on the complex planctomycetal physiology we analyzed carotenoid profiles of the pink-pigmented strain Rhodopirellula rubra LF2T and of the orange strain Rubinisphaera brasiliensis Gr7. During LC-MS/MS analysis of culture extracts we were able to identify three saproxanthin-type carotenoids including a rare C45 carotenoid. These compounds, saproxanthin, dehydroflexixanthin and 2’-isopentenyldehydrosaproxanthin, derive from the common carotenoid precursor lycopene and are characterized by related end groups, namely a 3-hydroxylated β-carotene-like cyclohexene ring as one end group and simple hydration on the other end of the molecule. Based on the observed molecule structure we present putative pathways for their biosynthesis. Results support Planctomycetes as a promising, yet mostly untapped source of carotenoids

S66: A Well-balanced Database of Benchmark Interaction Energies Relevant to Biomolecular Structures

With numerous new quantum chemistry methods being developed in recent years and the promise of even more new methods to be developed in the near future, it is clearly critical that highly accurate, well-balanced, reference data for many different atomic and molecular properties be available for the parametrization and validation of these methods. One area of research that is of particular importance in many areas of chemistry, biology, and material science is the study of noncovalent interactions. Because these interactions are often strongly influenced by correlation effects, it is necessary to use computationally expensive high-order wave function methods to describe them accurately. Here, we present a large new database of interaction energies calculated using an accurate CCSD(T)/CBS scheme. Data are presented for 66 molecular complexes, at their reference equilibrium geometries and at 8 points systematically exploring their dissociation curves; in total, the database contains 594 points: 66 at equilibrium geometries, and 528 in dissociation curves. The data set is designed to cover the most common types of noncovalent interactions in biomolecules, while keeping a balanced representation of dispersion and electrostatic contributions. The data set is therefore well suited for testing and development of methods applicable to bioorganic systems. In addition to the benchmark CCSD(T) results, we also provide decompositions of the interaction energies by means of DFT-SAPT calculations. The data set was used to test several correlated QM methods, including those parametrized specifically for noncovalent interactions. Among these, the SCS-MI-CCSD method outperforms all other tested methods, with a root-mean-square error of 0.08 kcal/mol for the S66 data set

Ab initio van der Waals interactions in simulations of water alter structure from mainly tetrahedral to high-density-like

The structure of liquid water at ambient conditions is studied in ab initio molecular dynamics simulations using van der Waals (vdW) density-functional theory, i.e. using the new exchange-correlation functionals optPBE-vdW and vdW-DF2. Inclusion of the more isotropic vdW interactions counteracts highly directional hydrogen-bonds, which are enhanced by standard functionals. This brings about a softening of the microscopic structure of water, as seen from the broadening of angular distribution functions and, in particular, from the much lower and broader first peak in the oxygen-oxygen pair-correlation function (PCF), indicating loss of structure in the outer solvation shells. In combination with softer non-local correlation terms, as in the new parameterization of vdW-DF, inclusion of vdW interactions is shown to shift the balance of resulting structures from open tetrahedral to more close-packed. The resulting O-O PCF shows some resemblance with experiment for high-density water (A. K. Soper and M. A. Ricci, Phys. Rev. Lett., 84:2881, 2000), but not directly with experiment for ambient water. However, an O-O PCF consisting of a linear combination of 70% from vdW-DF2 and 30% from experiment on low-density liquid water reproduces near-quantitatively the experimental O-O PCF for ambient water, indicating consistency with a two-liquid model with fluctuations between high- and low-density regions

Random-phase approximation and its applications in computational chemistry and materials science

The random-phase approximation (RPA) as an approach for computing the electronic correlation energy is reviewed. After a brief account of its basic concept and historical development, the paper is devoted to the theoretical formulations of RPA, and its applications to realistic systems. With several illustrating applications, we discuss the implications of RPA for computational chemistry and materials science. The computational cost of RPA is also addressed which is critical for its widespread use in future applications. In addition, current correction schemes going beyond RPA and directions of further development will be discussed.Comment: 25 pages, 11 figures, published online in J. Mater. Sci. (2012

Rationale, study design, and analysis plan of the Alveolar Recruitment for ARDS Trial (ART): study protocol for a randomized controlled trial

BACKGROUND: Acute respiratory distress syndrome (ARDS) is associated with high in-hospital mortality. Alveolar recruitment followed by ventilation at optimal titrated PEEP may reduce ventilator-induced lung injury and improve oxygenation in patients with ARDS, but the effects on mortality and other clinical outcomes remain unknown. This article reports the rationale, study design, and analysis plan of the Alveolar Recruitment for ARDS Trial (ART). METHODS/DESIGN: ART is a pragmatic, multicenter, randomized (concealed), controlled trial, which aims to determine if maximum stepwise alveolar recruitment associated with PEEP titration is able to increase 28-day survival in patients with ARDS compared to conventional treatment (ARDSNet strategy). We will enroll adult patients with ARDS of less than 72 h duration. The intervention group will receive an alveolar recruitment maneuver, with stepwise increases of PEEP achieving 45 cmH2O and peak pressure of 60 cmH2O, followed by ventilation with optimal PEEP titrated according to the static compliance of the respiratory system. In the control group, mechanical ventilation will follow a conventional protocol (ARDSNet). In both groups, we will use controlled volume mode with low tidal volumes (4 to 6 mL/kg of predicted body weight) and targeting plateau pressure 6430 cmH2O. The primary outcome is 28-day survival, and the secondary outcomes are: length of ICU stay; length of hospital stay; pneumothorax requiring chest tube during first 7 days; barotrauma during first 7 days; mechanical ventilation-free days from days 1 to 28; ICU, in-hospital, and 6-month survival. ART is an event-guided trial planned to last until 520 events (deaths within 28 days) are observed. These events allow detection of a hazard ratio of 0.75, with 90% power and two-tailed type I error of 5%. All analysis will follow the intention-to-treat principle. DISCUSSION: If the ART strategy with maximum recruitment and PEEP titration improves 28-day survival, this will represent a notable advance to the care of ARDS patients. Conversely, if the ART strategy is similar or inferior to the current evidence-based strategy (ARDSNet), this should also change current practice as many institutions routinely employ recruitment maneuvers and set PEEP levels according to some titration metho

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level.Comment: 24 pages, 5 figure