Accurate Quantum Monte Carlo Forces for Machine-Learned Force Fields: Ethanol as a Benchmark.

Quantum Monte Carlo (QMC) is a powerful method to calculate accurate energies and forces for molecular systems. In this work, we demonstrate how we can obtain accurate QMC forces for the fluxional ethanol molecule at room temperature by using either multideterminant Jastrow-Slater wave functions in variational Monte Carlo or just a single determinant in diffusion Monte Carlo. The excellent performance of our protocols is assessed against high-level coupled cluster calculations on a diverse set of representative configurations of the system. Finally, we train machine-learning force fields on the QMC forces and compare them to models trained on coupled cluster reference data, showing that a force field based on the diffusion Monte Carlo forces with a single determinant can faithfully reproduce coupled cluster power spectra in molecular dynamics simulations.

Effectiveness of Nifurtimox in the Treatment of Chagas Disease: a Long-Term Retrospective Cohort Study in Children and Adults.

Chagas disease (ChD), caused by Trypanosoma cruzi, has a global prevalence due to patient migration. However, despite its worldwide distribution, long-term follow-up efficacy studies with nifurtimox (NF) are scarce and have been conducted with only small numbers of patients. A retrospective study of a large cohort of ChD treated children and adults with NF. Treatment response was evaluated by clinical, parasitological, and serological after-treatment evaluation. A total of 289 patients were enrolled, of which 199 were children and 90 adults. At diagnosis, 89.6% of patients were asymptomatic. Overall, all symptomatic patients showed clinical improvement. At baseline, parasitemia was positive in 130 of 260 (50%) patients. All but one adult patient had cleared their parasitemia by the end of treatment. That patient was considered a treatment failure. Median follow-up time for children was 37.7 months, with an interquartile range of (IQR25-75 12.2 to 85.3), and for adults was 14.2 months (IQR25-75, 1.9 to 33.8). After treatment, a decrease of T. cruzi antibodies and seroconversion were observed in 34.6% of patients. The seroconversion profile showed that, the younger the patient, the higher the rate of seroconversion (log rank test; P value, <0.01). At least 20% seroreduction at 1 year follow-up was observed in 33.2% of patients. Nifurtimox was highly effective for ChD treatment. Patients had excellent treatment responses with fully resolved symptoms related to acute T. cruzi infection. Clearance of parasitemia and a decrease in T. cruzi antibodies were observed as markers of treatment response. This study reinforces the importance of treating patients during childhood since the treatment response was more marked in younger subjects. (This protocol was registered at ClinicalTrials.gov under registration number NCT04274101).

Adverse Events Associated with Nifurtimox Treatment for Chagas Disease in Children and Adults.

Nifurtimox (NF) is one of the only two drugs currently available for Chagas disease (ChD) treatment. However, data on NF safety are scarce, and many physicians defer or refuse NF treatment because of concerns about drug tolerance. In a retrospective study of adverse drug reactions (ADRs) associated with NF treatment of ChD, children received NF doses of 10 to 15 mg/kg/day for 60 to 90 days, and adults received 8 to 10 mg/kg/day for 30 days. A total of 215 children (median age, 2.6 years; range, 0 to 17 years) and 105 adults (median age, 34 years; range, 18 to 57 years) were enrolled. Overall, 127/320 (39.7%) patients developed ADRs, with an incidence of 64/105 adults and 63/215 children (odds ratio [OR] = 3.7; 95% confidence interval [CI], 2.2 to 6.3). We observed 215 ADRs, 131 in adults (median, 2 events/patient; interquartile range for the 25th to 75th percentiles [IQR25-75], 1 to 3) and 84 in children (median, 1 event/patient; IQR25-75 = 1 to 1.5) (Padjusted < 0.001). ADRs were mainly mild and moderate. Severe ADRs were infrequent (1.2% in children and 0.9% in adults). Nutritional, central nervous, and digestive systems were the most frequently affected, without differences between groups. Treatment was discontinued in 31/320 (9.7%) patients without differences between groups. However, ADR-related discontinuations occurred more frequently in adults than in children (OR = 5.5, 95% CI = 1.5 to 24). Our study supports the safety of NF for ChD treatment. Delaying NF treatment due to safety concerns does not seem to be supported by the evidence. (This study has been registered in ClinicalTrials.gov under identifier NCT04274101.).

Optimizing the Energy with Quantum Monte Carlo: A Lower Numerical Scaling for Jastrow-Slater Expansions.

We present an improved formalism for quantum Monte Carlo calculations of energy derivatives and properties (e.g., the interatomic forces), with a multideterminant Jastrow-Slater function. As a function of the number Ne of Slater determinants, the numerical scaling of O(Ne) per derivative we have recently reported is here lowered to O(Ne) for the entire set of derivatives. As a function of the number of electrons N, the scaling to optimize the wave function and the geometry of a molecular system is lowered to O(N3) + O(NNe), the same as computing the energy alone in the sampling process. The scaling is demonstrated on linear polyenes up to C60H62 and the efficiency of the method is illustrated with the structural optimization of butadiene and octatetraene with Jastrow-Slater wave functions comprising as many as 200 000 determinants and 60 000 parameters.

Imaginary time density-density correlations for two-dimensional electron gases at high density.

We evaluate imaginary time density-density correlation functions for two-dimensional homogeneous electron gases of up to 42 particles in the continuum using the phaseless auxiliary field quantum Monte Carlo method. We use periodic boundary conditions and up to 300 plane waves as basis set elements. We show that such methodology, once equipped with suitable numerical stabilization techniques necessary to deal with exponentials, products, and inversions of large matrices, gives access to the calculation of imaginary time correlation functions for medium-sized systems. We discuss the numerical stabilization techniques and the computational complexity of the methodology and we present the limitations related to the size of the systems on a quantitative basis. We perform the inverse Laplace transform of the obtained density-density correlation functions, assessing the ability of the phaseless auxiliary field quantum Monte Carlo method to evaluate dynamical properties of medium-sized homogeneous fermion systems.

Imaginary time correlations and the phaseless auxiliary field quantum Monte Carlo.

The phaseless Auxiliary Field Quantum Monte Carlo (AFQMC) method provides a well established approximation scheme for accurate calculations of ground state energies of many-fermions systems. Here we address the possibility of calculating imaginary time correlation functions with the phaseless AFQMC. We give a detailed description of the technique and test the quality of the results for static properties and imaginary time correlation functions against exact values for small systems.

Practical Schemes for Accurate Forces in Quantum Monte Carlo.

While the computation of interatomic forces has become a well-established practice within variational Monte Carlo (VMC), the use of the more accurate Fixed-Node Diffusion Monte Carlo (DMC) method is still largely limited to the computation of total energies on structures obtained at a lower level of theory. Algorithms to compute exact DMC forces have been proposed in the past, and one such scheme is also put forward in this work, but remain rather impractical due to their high computational cost. As a practical route to DMC forces, we therefore revisit here an approximate method, originally developed in the context of correlated sampling and named here the Variational Drift-Diffusion (VD) approach. We thoroughly investigate its accuracy by checking the consistency between the approximate VD force and the derivative of the DMC potential energy surface for the SiH and C2 molecules and employ a wide range of wave functions optimized in VMC to assess its robustness against the choice of trial function. We find that, for all but the poorest wave function, the discrepancy between force and energy is very small over all interatomic distances, affecting the equilibrium bond length obtained with the VD forces by less than 0.004 au. Furthermore, when the VMC forces are approximate due to the use of a partially optimized wave function, the DMC forces have smaller errors and always lead to an equilibrium distance in better agreement with the experimental value. We also show that the cost of computing the VD forces is only slightly larger than the cost of calculating the DMC energy. Therefore, the VD approximation represents a robust and efficient approach to compute accurate DMC forces, superior to the VMC counterparts.

Minimum energy pathways via quantum Monte Carlo.

We perform quantum Monte Carlo (QMC) calculations to determine minimum energy pathways of simple chemical reactions, and compare the computed geometries and reaction barriers with those obtained with density functional theory (DFT) and quantum chemistry methods. We find that QMC performs in general significantly better than DFT, being also able to treat cases in which DFT is inaccurate or even unable to locate the transition state. Since the wave function form employed here is particularly simple and can be transferred to larger systems, we suggest that a QMC approach is both viable and useful for reactions difficult to address by DFT and system sizes too large for high level quantum chemistry methods.

Excitation spectrum of a supersolid.

Conclusive experimental evidence of a supersolid phase in any known condensed matter system is presently lacking. On the other hand, a supersolid phase has been recently predicted for a system of spinless bosons in continuous space, interacting via a broad class of soft-core, repulsive potentials. Such an interaction can be engineered in assemblies of ultracold atoms, providing a well-defined pathway to the unambiguous observation of this fascinating phase of matter. In this Letter, we study by first principles computer simulations the elementary excitation spectrum of the supersolid, and show that it features two distinct modes, namely, a solidlike phonon and a softer collective excitation, related to broken translation and gauge symmetry, respectively.

Rotational spectrum of cyanoacetylene solvated with helium atoms.

The high resolution microwave spectra of He(N)-HCCCN clusters were studied in the size ranges of 1-18 and 25-31. In the absence of an accompanying infrared study, rotational excitation energies were computed by the reptation quantum Monte Carlo method and used to facilitate the search and assignment of R(0) transitions from N > 6, as well as R(1) transitions with N > 1. The assignments in the range of 25-31 are accurate to +/-2 cluster size units, with an essentially certain relative ordering. The rotational transition frequencies decrease with N = 1-6 and then show oscillatory behavior for larger cluster sizes, which is now recognized to be a manifestation of the onset and microscopic evolution of superfluidity. For cluster sizes beyond completion of the first solvation shell the rotational frequencies increase significantly above the large-droplet limit. This behavior, common to other linear molecules whose interaction with He features a strong nearly equatorial minimum, is analyzed using path integral Monte Carlo simulations. The He density in the incipient second solvation shell is shown to open a new channel for long permutation cycles, thus increasing the decoupling of the quantum solvent from the rotation of the dopant molecule.

Effects of thickness on the spin susceptibility of the two dimensional electron gas.

Using available quantum Monte Carlo predictions for a strictly 2D electron gas, we estimate the spin susceptibility of electrons in actual devices taking into account the effect of the finite transverse thickness and finding very good agreement with experiments. A weak disorder, as found in very clean devices and/or at densities not too low, just brings about a minor enhancement of the susceptibility.

Age, menopausal status and homocysteine levels in women around menopause.

OBJECTIVE: To obtain data on the relation between age, menopause and homocysteine levels in women around menopause, we collected data for a sample of about 500 women attending a menopause clinic in Milan, Italy. STUDY DESIGN: Eligible for the study were all women aged 45-75 years, never HRT users consecutively observed for the first time at the Menopause Center of the 1st Obstetric Gynecological Clinic of the University of Milan. Fasting blood samples for total homocysteine plasma levels were collected during the visit. Of the 490 study subjects, 107 were pre-menopausal and 383 post-menopausal. RESULTS: In the total series, the mean homocysteine level was 8.3 micromol/L (S.D. 3.7, range 3.2-48.8). The values increased from 7.8 micromol/L in women aged <47 years to 9.0 micromol/L in those aged >59. Among pre-menopausal women the mean homocysteine level was 7.7 micromol/L compared to 8.3 micromol/L in post-menopausal women: this difference disappeared on the adjusted values. In post-menopausal women only, no clear relationship emerged between years since menopause and homocysteine levels. CONCLUSION: In our population, age, and not menopausal status, was the main determinant of homocysteine levels in women around menopause.

Small para-hydrogen clusters doped with carbon monoxide: quantum Monte Carlo simulations and observed infrared spectra.

The structures and rotational dynamics of clusters of a single carbon monoxide molecule solvated in para-hydrogen, (paraH(2))(N)-CO, have been simulated for sizes up to N=17 using the reptation Monte Carlo technique. The calculations indicate the presence of two series of R(0) rotational transitions with J=1<--0 for cold clusters, similar to those predicted and observed in the case of He(N)-CO. Infrared spectra of these clusters have been observed in the region of the C-O stretch ( approximately 2143 cm(-1)) in a pulsed supersonic jet expansion using a tunable diode laser probe. With the help of the calculations, the observed R(0) rotational transitions have been assigned up to N=9 for the b-type series and N=14 for the a-type series. Theory and experiment agree rather well, except that theory tends to overestimate the b-type energies. The (paraH(2))(12)-CO cluster is calculated to be particularly stable and (relatively) rigid, corresponding to completion of the first solvation shell, and it is observed to have the strongest a-type transition.

Bridging the gap between small clusters and nanodroplets: spectroscopic study and computer simulation of carbon dioxide solvated with helium atoms.

High resolution infrared spectra of He(N)-CO2 clusters with N up to 17 have been studied in the region of the CO2 nu(3) fundamental band. The B rotational constant initially drops as expected for a normal molecule, reaching a minimum for N=5. Its subsequent rise for N=6 to 11 can be interpreted as the transition from a normal (though floppy) molecule to a quantum solvation regime. For N>13, the B value becomes approximately constant with a value about 17% larger than that measured in much larger helium nanodroplets. Quantum Monte Carlo calculations of pure rotational spectra are in excellent agreement with the measured B in this size range and complement the experimental study with detailed structural information. For a larger cluster size (N=30-50) the simulations show a clear sign of convergence towards the nanodroplet B value.

Risk factors for endometrial hyperplasia: results from a case-control study.

We analyzed epidemiologic characteristics of women at risk for endometrial hyperplasia (EH), using data from a case-control study. One hundred twenty nine women aged 35-73 (median 51 years) with histologically confirmed complex endometrial hyperplasia without atypies identified at the University of Milan during the period 1990-99 were examined. Controls were 258 non hysterectomized women aged 36-74 (median 52 years), admitted to a network of hospitals covering the same area where cases had been identified for conditions other than gynecological, malignant, or hormone-related. Cases with EH were more educated than controls (OR > 12 years of education vs. < 7: 2.8, 95% CI 1.7-4.8), more frequently obese (OR 2.7, 95% CI 1.5-5.0) and diabetic (OR 2.4, 95% CI 0.8-6.9). Parous women (OR 1.8) and women reporting induced abortions (OR 1.6) showed an increased risk of EH, but the associations were not statistically significant. Compared to premenopausal women, the OR of EH was 0.2 (95% 0.1-0.5) for postmenopausal ones. Compared to women reporting menopause at age 50 or less, the OR of endometrial hyperplasia was 1.5 (95% CI 0.6-3.5) and 2.2 (95%CI 0.7-6.7), respectively, in women with menopause at age 50-52 and > or = 53. Considering postmenopausal women only the OR was 3.1 (95% CI 1.1-9.3) for use of hormonal replacement therapy (HRT). We conclude that this study indicates that high education, obesity, diabetes, and HRT use increase the risk of endometrial hyperplasia.

Identification of AFLP and microsatellite markers linked with an aluminium tolerance gene in barley ( Hordeum vulgare L.).

Barley is the most sensitive among the cereals to aluminium (Al) stress and breeding for more tolerant cultivars is a priority. To enhance selection efficiency for Al tolerance in barley, PCR-based AFLP and microsatellite markers linked to a locus conferring tolerance to aluminium were identified. The study used F(2) progeny derived from a single cross between Yambla (moderately tolerant of Al) and WB229 (tolerant of Al) and developed hydroponic pulse-recovery screening methods to assess tolerance of phenotypes based on root growth. The segregation ratios of tolerant and sensitive genotypes and F(3) progeny testing suggest that a single major gene controlled Al tolerance ( Alt). In order to determine the chromosomal location of the Alt gene, we used the AFLP technique coupled with bulk segregant analysis. We evaluated tolerant and sensitive bulks using 30 combinations of EcoRI/ MseI primers, and 12 of these permitted differentiation of the sensitive and tolerant bulks. More than 1,000 amplified fragments were obtained, and 98 polymorphic bands were scored. AFLP analysis of wheat-barley chromosome addition lines indicated that the Alt gene was located on barley chromosome 4H. Four chromosome 4H-specific microsatellite markers (Bmac310, Bmag353, HVM68 and HVMCABG) were tightly linked to Alt. The large allelic variation detected with microsatellite marker Bmag353 allowed us to implement this marker for routine marker-assisted selection for Al tolerance, and 396 plants could be screened on a single gel.

Energy expenditure and nutritional adequacy of rehabilitation paraplegics with asymptomatic bacteriuria and pressure sores.

OBJECTIVE: To measure resting energy expenditure (REE) in a group of people with postacute paraplegia, quantify the impact of asymptomatic bacteriuria and pressure sore(s) on patients' metabolic rate, and estimate the adequacy of patients' nutritional intakes to preserve patients' protein levels. MATERIAL AND METHODS: Ten males with post-acute paraplegia aged 42.1+/-18.7 years. We evaluated: height, body mass index (BMI), resting energy expenditure (REE), total daily calorie requirement (E), 24-h urine creatinine excretion (Cru), creatinine index (CI), and nitrogen balance (NB). RESULTS: Subjects with paraplegia showed high erythrocyte sedimentation rates. As a group, they had normal resting calorie consumption when REE was normalized for unit of urine creatinine (REE/Cru), it was higher in patients than in controls. Six of the 10 patients had a low calorie intake: of these only three had a negative nitrogen balance. CONCLUSION: In conclusion, the resting energy expenditure of the subjects with significant bacteriuria and pressure sore(s) of 23.7 kcal/kg/day suggests that a large portion of patients may have an inadequate calorie protein intake to preserve their nutritional status. The clinical significance of this study is that 28.5 kcal/kg/day may be the lower calorie threshold to meet the metabolic demands of people with apyretic paraplegia with bacteriuria and pressure sore(s).

Condensation of helium in nanotube bundles

Helium atoms are strongly attracted to the interstitial channels within a bundle of carbon nanotubes. The strong corrugation of the axial potential within a channel can produce a lattice gas system wherein the weak mutual attraction between atoms in neighboring channels induces a transition to an anisotropic condensed phase. At low temperatures, the specific heat of the adsorbate phase (with fewer than 2% of the atoms) greatly exceeds that of the host.

Equation of state of solid 3He

We present results of diffusion Monte Carlo calculations for the bcc and hcp phases of solid 3He, using a recent ab initio interatomic potential, including two- and three-body terms. This potential is found to yield an equation of state for condensed 4He in excellent agreement with experiment, in a wide density range. For 3He, we find a systematic discrepancy, worth 0.7 K, between our computed equation of state and a commonly accepted experimental one. We attribute such a discrepancy to an improper choice of reference energy in the determination of the experimental equation of state.

Estroprogestin vs. gonadotrophin agonists plus estroprogestin in the treatment of endometriosis-related pelvic pain: a randomized trial. Gruppo Italiano per lo Studio dell'Endometriosi.

OBJECTIVE: This is a randomized clinical trial comparing estroprogestin (E/P) pill given for 12 months vs. gonadotrophin releasing hormone agonist (GNRHa) given for 4 months followed by E/P pill treatment for 8 months in the relief of endometriosis-related pelvic pain. METHODS: Eligible for the study were women with laparoscopically confirmed endometriosis and pelvic pain lasting 3-12 months after diagnosis. Eligible women were randomly assigned to treatment with E/P pill (gestroden 0.75 mg and ethynlestradiol 0.03 mg) for 12 months (47 patients) vs. tryptorelin 3.75 mg slow release every 28 days for 4 months followed by E/P pill for 8 months (55 patients). RESULTS: At baseline, dysmenorrhea was reported in 46 women allocated to E/P pill only (97.9%), and in all the 55 women allocated to GNRHa+E/P pill. The corresponding value at the 12 months follow-up visit was 14 subjects (35.9%) and 16 subjects (34.8%). The baseline median values of the multidimensional and analog scale were for dysmenorrhea 4 and 6 in the EP only and 3 and 6 in the GNRHa+E/P group. The corresponding value at the 12 months follow-up visit were 2 and 6 and 0 and 5. Non-menstrual pain was reported, respectively, at baseline and 12 month visit by 46 (97.9%) and 15 (38.5%) subjects in the E/P pill group and 49 (89.1%) and 17 (37.0%) of the GNRHa+E/P pill one. The baseline median values of the multidimensional and analog scale were for non-menstrual pain 3 and 5 in the E/P only and 2 and 6 in the GNRHa+E/P group. The corresponding values at the 12 month follow-up visit were 0 and 4 and 0 and 4. These differences between the two groups were not statistically significant. CONCLUSIONS: 1 year after randomization, the two treatment schedules show similar relief of pelvic pain in women with endometriosis.