Metal-Insulator Transition in a Semiconductor Heterobilayer Model.

Transition metal dichalcogenide superlattices provide an exciting new platform for exploring and understanding a variety of phases of matter. The moiré continuum Hamiltonian, of two-dimensional jellium in a modulating potential, provides a fundamental model for such systems. Accurate computations with this model are essential for interpreting experimental observations and making predictions for future explorations. In this work, we combine two complementary quantum Monte Carlo (QMC) methods, phaseless auxiliary field quantum Monte Carlo and fixed-phase diffusion Monte Carlo, to study the ground state of this Hamiltonian. We observe a metal-insulator transition between a paramagnet and a 120° Néel ordered state as the moiré potential depth and the interaction strength are varied. We find significant differences from existing results by Hartree-Fock and exact diagonalization studies. In addition, we benchmark density-functional theory, and suggest an optimal hybrid functional which best approximates our QMC results.

Association between diabetes and cancer. Current mechanistic insights into the association and future challenges.

Compelling pieces of epidemiological, clinical, and experimental research have demonstrated that Diabetes mellitus (DM) is a major risk factor associated with increased cancer incidence and mortality in many human neoplasms. In the pathophysiology context of DM, many of the main classical actors are relevant elements that can fuel the different steps of the carcinogenesis process. Hyperglycemia, hyperinsulinemia, metabolic inflammation, and dyslipidemia are among the classic contributors to this association. Furthermore, new emerging actors have received particular attention in the last few years, and compelling data support that the microbiome, the epigenetic changes, the reticulum endoplasmic stress, and the increased glycolytic influx also play important roles in promoting the development of many cancer types. The arsenal of glucose-lowering therapeutic agents used for treating diabetes is wide and diverse, and a growing body of data raised during the last two decades has tried to clarify the contribution of therapeutic agents to this association. However, this research area remains controversial, because some anti-diabetic drugs are now considered as either promotors or protecting elements. In the present review, we intend to highlight the compelling epidemiological shreds of evidence that support this association, as well as the mechanistic contributions of many of these potential pathological mechanisms, some controversial points as well as future challenges.

Long-term potentiation and its neurotrophin-dependent modulation in the superior cervical ganglion of the rat are influenced by KCNQ channel function.

Ganglionic long-term potentiation (gLTP) in the rat superior cervical ganglion (SCG) is differentially modulated by neurotrophic factors (Nts): brain-derived neurotrophic factor (BDNF) and nerve growth factor (NGF). KCNQ/M channels, key regulators of neuronal excitability, and firing pattern are modulated by Nts; therefore, they might contribute to gLTP expression and to the Nts-dependent modulation of gLTP. In the SCG of rats, we characterized the presence of the KCNQ2 isoform and the effects of opposite KCNQ/M channel modulators on gLTP in control condition and under Nts modulation. Immunohistochemical and reverse transcriptase polymerase chain reaction analyses showed the expression of the KCNQ2 isoform. We found that 1 µmol/L XE991, a channel inhibitor, significantly reduced gLTP (∼50%), whereas 5 µmol/L flupirtine, a channel activator, significantly increased gLTP (1.3- to 1.7-fold). Both modulators counterbalanced the effects of the Nts on gLTP. Data suggest that KCNQ/M channels are likely involved in gLTP expression and in the modulation exerted by BDNF and NGF.

Optical properties of high-pressure fluid hydrogen across molecular dissociation.

Optical properties of compressed fluid hydrogen in the region where dissociation and metallization is observed are computed by ab initio methods and compared with recent experimental results. We confirm that at T > 3,000 K, both processes are continuous, while at T < 1,500 K, the first-order phase transition is accompanied by a discontinuity of the dc conductivity and the thermal conductivity, while both the reflectivity and absorption coefficient vary rapidly but continuously. Our results support the recent analysis of National Ignition Facility (NIF) experiments [Celliers PM, et al. (2018) Science 361:677-682], which assigned the inception of metallization to pressures where the reflectivity is ∼0.3. Our results also support the conclusion that the temperature plateau seen in laser-heated diamond-anvil cell (DAC) experiments at temperatures higher than 1,500 K corresponds to the onset of optical absorption, not to the phase transition.

Reduction in the potential distribution of bumble bees (Apidae: Bombus) in Mesoamerica under different climate change scenarios: Conservation implications.

Bumble bees are an ecologically and economically important group of pollinating insects worldwide. Global climate change is predicted to affect bumble bee ecology including habitat suitability and geographic distribution. Our study aims to estimate the impact of projected climate change on 18 Mesoamerican bumble bee species. We used ecological niche modeling (ENM) using current and future climate emissions scenarios (representative concentration pathway 4.5, 6.0, and 8.5) and models (CCSM4, HadGEM2-AO, and MIROC-ESM-CHEM). Regardless of the scenario and model applied, our results suggest that all bumble bee species are predicted to undergo a reduction in their potential distribution and habitat suitability due to projected climate change. ENMs based on low emission scenarios predict a distribution loss ranging from 7% to 67% depending on the species for the year 2050. Furthermore, we discovered that the reduction of bumble bee geographic range shape will be more evident at the margins of their distribution. The reduction of suitable habitat is predicted to be accompanied by a 100-500 m upslope change in altitude and 1-581 km shift away from the current geographic centroid of a species' distribution. On average, protected natural areas in Mesoamerica cover ~14% of each species' current potential distribution, and this proportion is predicted to increase to ~23% in the high emission climate change scenarios. Our models predict that climate change will reduce Mesoamerican bumble bee habitat suitability, especially for rare species, by reducing their potential distribution ranges and suitability. The small proportion of current and future potential distribution falling in protected natural areas suggests that such areas will likely have marginal contribution to bumble bee habitat conservation. Our results have the capacity to inform stakeholders in designing effective landscape management for bumble bees, which may include developing restoration plans for montane pine oak forests habitats and native flowering plants.

A phaseless auxiliary-field quantum Monte Carlo perspective on the uniform electron gas at finite temperatures: Issues, observations, and benchmark study.

We investigate the viability of the phaseless finite-temperature auxiliary-field quantum Monte Carlo (ph-FT-AFQMC) method for ab initio systems using the uniform electron gas as a model. Through comparisons with exact results and FT coupled cluster theory, we find that ph-FT-AFQMC is sufficiently accurate at high to intermediate electronic densities. We show, both analytically and numerically, that the phaseless constraint at FT is fundamentally different from its zero-temperature counterpart (i.e., ph-ZT-AFQMC), and generally, one should not expect ph-FT-AFQMC to agree with ph-ZT-AFQMC in the low-temperature limit. With an efficient implementation, we are able to compare exchange-correlation energies to the existing results in the thermodynamic limit and find that the existing parameterizations are highly accurate. In particular, we found that ph-FT-AFQMC exchange-correlation energies are in better agreement with a known parameterization than is restricted path-integral MC in the regime of Θ ≤ 0.5 and rs ≤ 2, which highlights the strength of ph-FT-AFQMC.

Accelerating the convergence of auxiliary-field quantum Monte Carlo in solids with optimized Gaussian basis sets.

We investigate the use of optimized correlation-consistent Gaussian basis sets for the study of insulating solids with auxiliary-field quantum Monte Carlo (AFQMC). The exponents of the basis set are optimized through the minimization of the second-order Møller-Plesset perturbation theory (MP2) energy in a small unit cell of the solid. We compare against other alternative basis sets proposed in the literature, namely, calculations in the Kohn-Sham basis and in the natural orbitals of an MP2 calculation. We find that our optimized basis sets accelerate the convergence of the AFQMC correlation energy compared to a Kohn-Sham basis and offer similar convergence to MP2 natural orbitals at a fraction of the cost needed to generate them. We also suggest the use of an improved, method independent, MP2-based basis set correction that significantly reduces the required basis set sizes needed to converge the correlation energy. With these developments, we study the relative performance of these basis sets in LiH, Si, and MgO and determine that our optimized basis sets yield the most consistent results as a function of volume. Using these optimized basis sets, we systematically converge the AFQMC calculations to the complete basis set and thermodynamic limit and find excellent agreement with experiment for the systems studied. Although we focus on AFQMC, our basis set generation procedure is independent of the subsequent correlated wavefunction method used.

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion quantum Monte Carlo.

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing the accuracy. Advances in real space methods include techniques for accurate computation of bandgaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods, including GW and density functional based techniques. To provide an improved foundation for these calculations, we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

Non-orthogonal multi-Slater determinant expansions in auxiliary field quantum Monte Carlo.

We investigate the use of non-orthogonal multi-Slater determinant (NOMSD) expansions as trial wavefunctions in auxiliary field quantum Monte Carlo simulations of molecular systems. We show that NOMSD trial wavefunctions with as few as twenty determinants are sufficient in order to achieve chemical accuracy across most of the G1 molecular test set. We also show that NOMSD trial wavefunctions are useful for more challenging strongly correlated systems by computing relative energies along the isomerization path of the [Cu2O2]2+ molecule. Our results for [Cu2O2]2+ compare favorably with other accurate quantum many-body methods, including density matrix renormalization group and completely renormalized coupled cluster methods.

Liquid-liquid phase transition in hydrogen by coupled electron-ion Monte Carlo simulations.

The phase diagram of high-pressure hydrogen is of great interest for fundamental research, planetary physics, and energy applications. A first-order phase transition in the fluid phase between a molecular insulating fluid and a monoatomic metallic fluid has been predicted. The existence and precise location of the transition line is relevant for planetary models. Recent experiments reported contrasting results about the location of the transition. Theoretical results based on density functional theory are also very scattered. We report highly accurate coupled electron-ion Monte Carlo calculations of this transition, finding results that lie between the two experimental predictions, close to that measured in diamond anvil cell experiments but at 25-30 GPa higher pressure. The transition along an isotherm is signaled by a discontinuity in the specific volume, a sudden dissociation of the molecules, a jump in electrical conductivity, and loss of electron localization.

Ganglionic Long-Term Potentiation in Prehypertensive and Hypertensive Stages of Spontaneously Hypertensive Rats Depends on GABA Modulation.

The sympathetic nervous system (SNS) regulates body functions in normal and pathological conditions and is characterized by the presence of a neuroplastic phenomenon, termed ganglionic long-term potentiation (gLTP). In hypertension, either in spontaneously hypertensive rats (SHR) or in humans, sympathetic hyperfunction, such as elevated SNS outflow and changes in synaptic plasticity have been described. Because enhanced SNS outflow is detected in the hypertensive stage and, more importantly, in the prehypertensive phase of SHR, here we explored whether synaptic plasticity, particularly gLTP, was modified in the superior cervical ganglia (SCG) of prehypertensive SHR. Furthermore, considering that GABA modulates sympathetic synaptic transmission and gLTP in Wistar rats, we studied whether GABA might modulate gLTP expression in SHR. We characterized gLTP in the SCG of young prehypertensive 6-week-old (wo) and adult hypertensive (12 wo) SHR and in the SCG of Wistar Kyoto (WKy) normotensive control rats of the same ages. We found that gLTP was expressed in 6 wo SHR, but not in 12 wo rats. By contrast, in WKy, gLTP was expressed in 12 wo, but not in 6 wo rats. We also found that gLTP depends on GABA modulation, as blockade of GABA-A subtype receptors with its antagonist bicuculline unmasked gLTP expression in adult SHR and young WKy. We propose that (1) activity-dependent changes in synaptic efficacy are altered not only during hypertension but also before its onset and (2) GABA may play a modulatory role in the changes in synaptic plasticity in SHR, because the blockade of GABA-A receptors unmasked the expression of gLTP. These early changes in neuroplasticity and GABA modulation of gLTP could be part of the sympathetic hyperfunction observed in hypertension.

Nuclear Quantum Effects in Water and Aqueous Systems: Experiment, Theory, and Current Challenges.

Nuclear quantum effects influence the structure and dynamics of hydrogen-bonded systems, such as water, which impacts their observed properties with widely varying magnitudes. This review highlights the recent significant developments in the experiment, theory, and simulation of nuclear quantum effects in water. Novel experimental techniques, such as deep inelastic neutron scattering, now provide a detailed view of the role of nuclear quantum effects in water's properties. These have been combined with theoretical developments such as the introduction of the principle of competing quantum effects that allows the subtle interplay of water's quantum effects and their manifestation in experimental observables to be explained. We discuss how this principle has recently been used to explain the apparent dichotomy in water's isotope effects, which can range from very large to almost nonexistent depending on the property and conditions. We then review the latest major developments in simulation algorithms and theory that have enabled the efficient inclusion of nuclear quantum effects in molecular simulations, permitting their combination with on-the-fly evaluation of the potential energy surface using electronic structure theory. Finally, we identify current challenges and future opportunities in this area of research.

Auxiliary-field quantum Monte Carlo calculations of the structural properties of nickel oxide.

Auxiliary-field quantum Monte Carlo (AFQMC) has repeatedly demonstrated itself as one of the most accurate quantum many-body methods, capable of simulating both real and model systems. In this article, we investigate the application of AFQMC to realistic strongly correlated materials in periodic Gaussian basis sets. Using nickel oxide (NiO) as an example, we investigate the importance of finite size effects and basis set errors on the structural properties of the correlated solid. We provide benchmark calculations for NiO and compare our results to both experimental measurements and existing theoretical methods.

Coupled electron-ion Monte Carlo simulation of hydrogen molecular crystals.

We performed simulations for solid molecular hydrogen at high pressures (250 GPa ≤ P ≤ 500 GPa) along two isotherms at T = 200 K (phase III) and at T = 414 K (phase IV). At T = 200 K, we considered likely candidates for phase III, the C2c and Cmca12 structures, while at T = 414 K in phase IV, we studied the Pc48 structure. We employed both Coupled Electron-Ion Monte Carlo (CEIMC) and Path Integral Molecular Dynamics (PIMD). The latter is based on Density Functional Theory (DFT) with the van der Waals approximation (vdW-DF). The comparison between the two methods allows us to address the question of the accuracy of the exchange-correlation approximation of DFT for thermal and quantum protons without recurring to perturbation theories. In general, we find that atomic and molecular fluctuations in PIMD are larger than in CEIMC which suggests that the potential energy surface from vdW-DF is less structured than the one from quantum Monte Carlo. We find qualitatively different behaviors for systems prepared in the C2c structure for increasing pressure. Within PIMD, the C2c structure is dynamically partially stable for P ≤ 250 GPa only: it retains the symmetry of the molecular centers but not the molecular orientation; at intermediate pressures, it develops layered structures like Pbcn or Ibam and transforms to the metallic Cmca-4 structure at P ≥ 450 GPa. Instead, within CEIMC, the C2c structure is found to be dynamically stable at least up to 450 GPa; at increasing pressure, the molecular bond length increases and the nuclear correlation decreases. For the other two structures, the two methods are in qualitative agreement although quantitative differences remain. We discuss various structural properties and the electrical conductivity. We find that these structures become conducting around 350 GPa but the metallic Drude-like behavior is reached only at around 500 GPa, consistent with recent experimental claims.

Leishmania phosphatase PP5 is a regulator of HSP83 phosphorylation and essential for parasite pathogenicity.

Leishmania parasites are responsible for important neglected diseases in humans and animals, ranging from self-healing cutaneous lesions to fatal visceral manifestations. During the infectious cycle, Leishmania differentiates from the extracellular flagellated promastigote to the intracellular pathogenic amastigote. Parasite differentiation is triggered by changes in environmental cues, mainly pH and temperature. In general, extracellular signals are translated into stage-specific gene expression by a cascade of reversible protein phosphorylation regulated by protein kinases and phosphatases. Though protein kinases have been actively studied as potential anti-parasitic drug targets, our understanding of the biology of protein phosphatases in Leishmania is poor. We have previously reported the principal analysis of a novel protein phosphatase 5 (PP5) in Leishmania species. Here, we assessed the role of PP5 in parasite pathogenicity, where we uncovered, using transgenic PP5 over-expressing and PP5 null-mutant parasites, its importance in metacyclogeneisis, maintaining HSP83 phosphorylation homeostasis and virulence. All together, our results indicate the importance of PP5 in regulating parasite stress and adaptation during differentiation, making this protein an attractive potential target for therapeutic intervention.

Cellulases immobilization on chitosan-coated magnetic nanoparticles: application for Agave Atrovirens lignocellulosic biomass hydrolysis.

In the present study, Trichoderma reesei cellulase was covalently immobilized on chitosan-coated magnetic nanoparticles using glutaraldehyde as a coupling agent. The average diameter of magnetic nanoparticles before and after enzyme immobilization was about 8 and 10 nm, respectively. The immobilized enzyme retained about 37 % of its initial activity, and also showed better thermal and storage stability than free enzyme. Immobilized cellulase retained about 80 % of its activity after 15 cycles of carboxymethylcellulose hydrolysis and was easily separated with the application of an external magnetic field. However, in this reaction, K m was increased eight times. The immobilized enzyme was able to hydrolyze lignocellulosic material from Agave atrovirens leaves with yield close to the amount detected with free enzyme and it was re-used in vegetal material conversion up to four cycles with 50 % of activity decrease. This provides an opportunity to reduce the enzyme consumption during lignocellulosic material saccharification for bioethanol production.

Zinc(II)-Dipicolylamine Coordination Complexes as Targeting and Chemotherapeutic Agents for Leishmania major.

Cutaneous leishmaniasis is a neglected tropical disease that causes painful lesions and severe disfigurement. Modern treatment relies on a few chemotherapeutics with serious limitations, and there is a need for more effective alternatives. This study describes the selective targeting of zinc(II)-dipicolylamine (ZnDPA) coordination complexes toward Leishmania major, one of the species responsible for cutaneous leishmaniasis. Fluorescence microscopy of L. major promastigotes treated with a fluorescently labeled ZnDPA probe indicated rapid accumulation of the probe within the axenic promastigote cytosol. The antileishmanial activities of eight ZnDPA complexes were measured using an in vitro assay. All tested complexes exhibited selective toxicity against L. major axenic promastigotes, with 50% effective concentration values in the range of 12.7 to 0.3 µM. Similar toxicity was observed against intracellular amastigotes, but there was almost no effect on the viability of mammalian cells, including mouse peritoneal macrophages. In vivo treatment efficacy studies used fluorescence imaging to noninvasively monitor changes in the red fluorescence produced by an infection of mCherry-L. major in a mouse model. A ZnDPA treatment regimen reduced the parasite burden nearly as well as the reference care agent, potassium antimony(III) tartrate, and with less necrosis in the local host tissue. The results demonstrate that ZnDPA coordination complexes are a promising new class of antileishmanial agents with potential for clinical translation.

Exploring the free energy surface using ab initio molecular dynamics.

Efficient exploration of configuration space and identification of metastable structures in condensed phase systems are challenging from both computational and algorithmic perspectives. In this regard, schemes that utilize a set of pre-defined order parameters to sample the relevant parts of the configuration space [L. Maragliano and E. Vanden-Eijnden, Chem. Phys. Lett. 426, 168 (2006); J. B. Abrams and M. E. Tuckerman, J. Phys. Chem. B 112, 15742 (2008)] have proved useful. Here, we demonstrate how these order-parameter aided temperature accelerated sampling schemes can be used within the Born-Oppenheimer and the Car-Parrinello frameworks of ab initio molecular dynamics to efficiently and systematically explore free energy surfaces, and search for metastable states and reaction pathways. We have used these methods to identify the metastable structures and reaction pathways in SiO2 and Ti. In addition, we have used the string method [W. E, W. Ren, and E. Vanden-Eijnden, Phys. Rev. B 66, 052301 (2002); L. Maragliano et al., J. Chem. Phys. 125, 024106 (2006)] within the density functional theory to study the melting pathways in the high pressure cotunnite phase of SiO2 and the hexagonal closed packed to face centered cubic phase transition in Ti.

On the accuracy of the MB-pol many-body potential for water: Interaction energies, vibrational frequencies, and classical thermodynamic and dynamical properties from clusters to liquid water and ice.

The MB-pol many-body potential has recently emerged as an accurate molecular model for water simulations from the gas to the condensed phase. In this study, the accuracy of MB-pol is systematically assessed across the three phases of water through extensive comparisons with experimental data and high-level ab initio calculations. Individual many-body contributions to the interaction energies as well as vibrational spectra of water clusters calculated with MB-pol are in excellent agreement with reference data obtained at the coupled cluster level. Several structural, thermodynamic, and dynamical properties of the liquid phase at atmospheric pressure are investigated through classical molecular dynamics simulations as a function of temperature. The structural properties of the liquid phase are in nearly quantitative agreement with X-ray diffraction data available over the temperature range from 268 to 368 K. The analysis of other thermodynamic and dynamical quantities emphasizes the importance of explicitly including nuclear quantum effects in the simulations, especially at low temperature, for a physically correct description of the properties of liquid water. Furthermore, both densities and lattice energies of several ice phases are also correctly reproduced by MB-pol. Following a recent study of DFT models for water, a score is assigned to each computed property, which demonstrates the high and, in many respects, unprecedented accuracy of MB-pol in representing all three phases of water.

Probing druggability and biological function of essential proteins in Leishmania combining facilitated null mutant and plasmid shuffle analyses.

Leishmania parasites cause important human morbidity and mortality. Essential Leishmania genes escape genetic assessment by loss-of-function analyses due to lethal null mutant phenotypes, even though these genes and their products are biologically most significant and represent validated drug targets. Here we overcome this limitation using a facilitated null mutant approach applied for the functional genetic analysis of the MAP kinase LmaMPK4. This system relies on the episomal expression of the target gene from vector pXNG that expresses the Herpes simplex virus thymidine kinase gene thus rendering transgenic parasites susceptible for negative selection using the antiviral drug ganciclovir. Using this system we establish the genetic proof of LmaMPK4 as essential kinase in promastigotes. LmaMPK4 structure/function analysis by plasmid shuffle allowed us to identify regulatory kinase sequence elements relevant for chemotherapeutic intervention. A partial null mutant, expressing an MPK4 derivative with altered ATP-binding properties, showed defects in metacyclogenesis, establishing a first link of MPK4 function to parasite differentiation. The approaches presented here are broadly applicable to any essential gene in Leishmania thus overcoming major bottlenecks for their functional genetic analysis and their exploitation for structure-informed drug development.