Determinant- and Derivative-Free Quantum Monte Carlo Within the Stochastic Representation of Wavefunctions.

Describing the ground states of continuous, real-space quantum many-body systems, like atoms and molecules, is a significant computational challenge with applications throughout the physical sciences. Recent progress was made by variational methods based on machine learning (ML) ansatzes. However, since these approaches are based on energy minimization, ansatzes must be twice differentiable. This (a) precludes the use of many powerful classes of ML models; and (b) makes the enforcement of bosonic, fermionic, and other symmetries costly. Furthermore, (c) the optimization procedure is often unstable unless it is done by imaginary time propagation, which is often impractically expensive in modern ML models with many parameters. The stochastic representation of wavefunctions (SRW), introduced in Nat Commun 14, 3601 (2023), is a recent approach to overcoming (c). SRW enables imaginary time propagation at scale, and makes some headway towards the solution of problem (b), but remains limited by problem (a). Here, we argue that combining SRW with path integral techniques leads to a new formulation that overcomes all three problems simultaneously. As a demonstration, we apply the approach to generalized ``Hooke's atoms'': interacting particles in harmonic wells. We benchmark our results against state-of-the-art data where possible, and use it to investigate the crossover between the Fermi liquid and the Wigner molecule within closed-shell systems. Our results shed new light on the competition between interaction-driven symmetry breaking and kinetic-energy-driven delocalization.

Stark Many-Body Localization in Interacting Infinite Dimensional Systems.

We study bulk particle transport in a Fermi-Hubbard model on an infinite-dimensional Bethe lattice, driven by a constant electric field. Previous numerical studies showed that one dimensional analogs of this system exhibit a breakdown of diffusion due to Stark many-body localization at least up to time that scales exponentially with the system size. Here, we consider systems initially in a spin density wave state using a combination of numerically exact and approximate techniques. We show that for sufficiently weak electric fields, the wave's momentum component decays exponentially with time in a way consistent with normal diffusion. By studying different wavelengths, we extract the dynamical exponent and the generalized diffusion coefficient at each field strength. Interestingly, we find a nonmonotonic dependence of the dynamical exponent on the electric field. As the field increases toward a critical value proportional to the Hubbard interaction strength, transport slows down, becoming subdiffusive. At large interaction strengths, however, transport speeds up again with increasing field, exhibiting superdiffusive characteristics when the electric field is comparable to the interaction strength. Eventually, at the large field limit, localization occurs and the current through the system is suppressed.

Numerically Exact Simulation of Photodoped Mott Insulators.

A description of long-lived photodoped states in Mott insulators is challenging, as it needs to address exponentially separated timescales. We demonstrate how properties of such states can be computed using numerically exact steady state techniques, in particular, the quantum Monte Carlo algorithm, by using a time-local ansatz for the distribution function with separate Fermi functions for the electron and hole quasiparticles. The simulations show that the Mott gap remains robust to large photodoping, and the photodoped state has hole and electron quasiparticles with strongly renormalized properties.

Nonequilibrium steady state full counting statistics in the noncrossing approximation.

Quantum transport is often characterized not just by mean observables like the particle or energy current but by their fluctuations and higher moments, which can act as detailed probes of the physical mechanisms at play. However, relatively few theoretical methods are able to access the full counting statistics (FCS) of transport processes through electronic junctions in strongly correlated regimes. While most experiments are concerned with steady state properties, most accurate theoretical methods rely on computationally expensive propagation from a tractable initial state. Here, we propose a simple approach for computing the FCS through a junction directly at the steady state, utilizing the propagator noncrossing approximation. Compared to time propagation, our method offers reduced computational cost at the same level of approximation, but the idea can also be used within other approximations or as a basis for numerically exact techniques. We demonstrate the method's capabilities by investigating the impact of lead dimensionality on electronic transport in the nonequilibrium Anderson impurity model at the onset of Kondo physics. Our results reveal a distinct signature of one dimensional leads in the noise and Fano factor not present for other dimensionalities, showing the potential of FCS measurements as a probe of the environment surrounding a quantum dot.

Shaping Electronic Flows with Strongly Correlated Physics.

Nonequilibrium quantum transport is of central importance in nanotechnology. Its description requires the understanding of strong electronic correlations that couple atomic-scale phenomena to the nanoscale. So far, research in correlated transport has focused predominantly on few-channel transport, precluding the investigation of cross-scale effects. Recent theoretical advances enable the solution of models that capture the interplay between quantum correlations and confinement beyond a few channels. This problem is the focus of this study. We consider an atomic impurity embedded in a metallic nanosheet spanning two leads, showing that transport is significantly altered by tuning only the phase of a single local hopping parameter. Furthermore─depending on this phase─correlations reshape the electronic flow throughout the sheet, either funneling it through the impurity or scattering it away from a much larger region. This demonstrates the potential for quantum correlations to bridge length scales in the design of nanoelectronic devices and sensors.

Benzyl-para-di-[5-methyl-4-(n-octylamino) pyrimidin-2(1H)one] as an interferon beta (IFN-ß) modulator.

IFN-ß is a cytokine that plays a significant role in the immune system. Inhibition of IFN-ß might be used as a therapeutic approach to treat septic shock. A peptidomimetic previously developed by our research team, 1-benzyl-5-methyl-4-(n-octylamino)pyrimidin-2(1H)-one (LT87), was used as an cardioprotective agent in a myocardial ischemia (MI) mouse model. We have developed new LT87 derivatives by synthetizing its dimers in an attempt to extend its structural variety and enhance its biological activity. A dimeric derivative, LT127, exhibited a dose-dependent inhibition of LPS-mediated IFN-ß and subsequent CXCL10 mRNA transcription. The effect was selective and transduced through TLR4- and TRAM/TRIF-mediated signaling, with no significant effect on MyD88-dependent signaling. However, this effect was not specific to TLR4, since a similar effect was observed both on TLR8- and MDA5/RIG-I-stimulated IFN-ß expression. Nevertheless, LT127 might serve as a drug candidate, specifically as an inhibitor for IFN-ß production in order to develop a novel therapeutic approach to prevent septic shock.

Dynamic graph and polynomial chaos based models for contact tracing data analysis and optimal testing prescription.

In this study, we address three important challenges related to disease transmissions such as the COVID-19 pandemic, namely, (a) providing an early warning to likely exposed individuals, (b) identifying individuals who are asymptomatic, and (c) prescription of optimal testing when testing capacity is limited. First, we present a dynamic-graph based SEIR epidemiological model in order to describe the dynamics of the disease propagation. Our model considers a dynamic graph/network that accounts for the interactions between individuals over time, such as the ones obtained by manual or automated contact tracing, and uses a diffusion-reaction mechanism to describe the state dynamics. This dynamic graph model helps identify likely exposed/infected individuals to whom we can provide early warnings, even before they display any symptoms and/or are asymptomatic. Moreover, when the testing capacity is limited compared to the population size, reliable estimation of individual's health state and disease transmissibility using epidemiological models is extremely challenging. Thus, estimation of state uncertainty is paramount for both eminent risk assessment, as well as for closing the tracing-testing loop by optimal testing prescription. Therefore, we propose the use of arbitrary Polynomial Chaos Expansion, a popular technique used for uncertainty quantification, to represent the states, and quantify the uncertainties in the dynamic model. This design enables us to assign uncertainty of the state of each individual, and consequently optimize the testing as to reduce the overall uncertainty given a constrained testing budget. These tools can also be used to optimize vaccine distribution to curb the disease spread when limited vaccines are available. We present a few simulation results that illustrate the performance of the proposed framework, and estimate the impact of incomplete contact tracing data.

Anti-Herpes Simplex 1 Activity of Simmondsia chinensis (Jojoba) Wax.

Jojoba (Simmondsia chinensis (Link) Schneider) wax is used for various dermatological and pharmaceutical applications. Several reports have previously shown beneficial properties of Jojoba wax and extracts, including antimicrobial activity. The current research aimed to elucidate the impact of Jojoba wax on skin residential bacterial (Staphylococcus aureus and Staphylococcus epidermidis), fungal (Malassezia furfur), and virus infection (herpes simplex 1; HSV-1). First, the capacity of four commercial wax preparations to attenuate their growth was evaluated. The results suggest that the growth of Staphylococcus aureus, Staphylococcus epidermidis, and Malassezia furfur was unaffected by Jojoba in pharmacologically relevant concentrations. However, the wax significantly attenuated HSV-1 plaque formation. Next, a complete dose-response analysis of four different Jojoba varieties (Benzioni, Shiloah, Hatzerim, and Sheva) revealed a similar anti-viral effect with high potency (EC50 of 0.96 ± 0.4 µg/mL) that blocked HSV-1 plaque formation. The antiviral activity of the wax was also confirmed by real-time PCR, as well as viral protein expression by immunohistochemical staining. Chemical characterization of the fatty acid and fatty alcohol composition was performed, showing high similarity between the wax of the investigated varieties. Lastly, our results demonstrate that the observed effects are independent of simmondsin, repeatedly associated with the medicinal impact of Jojoba wax, and that Jojoba wax presence is required to gain protection against HSV-1 infection. Collectively, our results support the use of Jojoba wax against HSV-1 skin infections.

Multiorbital Quantum Impurity Solver for General Interactions and Hybridizations.

We present a numerically exact inchworm Monte Carlo method for equilibrium multiorbital quantum impurity problems with general interactions and hybridizations. We show that the method, originally developed to overcome the dynamical sign problem in certain real-time propagation problems, can also overcome the sign problem as a function of temperature for equilibrium quantum impurity models. This is shown in several cases where the current method of choice, the continuous-time hybridization expansion, fails due to the sign problem. Our method therefore enables simulations of impurity problems as they appear in embedding theories without further approximations, such as the truncation of the hybridization or interaction structure or a discretization of the impurity bath with a set of discrete energy levels, and eliminates a crucial bottleneck in the simulation of ab initio embedding problems.

Green's function methods for single molecule junctions.

We present a brief pedagogical review of theoretical Green's function methods applicable to open quantum systems out of equilibrium, in general, and single molecule junctions, in particular. We briefly describe experimental advances in molecular electronics and then discuss different theoretical approaches. We then focus on Green's function methods. Two characteristic energy scales governing the physics are many-body interactions within the junctions and molecule-contact coupling. We, therefore, discuss weak interactions and weak coupling as two limits that can be conveniently treated within, respectively, the standard nonequilibrium Green's function (NEGF) method and its many-body flavors (pseudoparticle and Hubbard NEGF). We argue that the intermediate regime, where the two energy scales are comparable, can in many cases be efficiently treated within the recently introduced superperturbation dual fermion approach. Finally, we review approaches for going beyond these analytically accessible limits, as embodied by recent developments in numerically exact methods based on Green's functions.

Advances in Understanding the Initial Steps of Pruritoceptive Itch: How the Itch Hits the Switch.

Pruritoceptive (dermal) itch was long considered an accompanying symptom of diseases, a side effect of drug applications, or a temporary sensation induced by invading pruritogens, as produced by the stinging nettle. Due to extensive research in recent years, it was possible to provide detailed insights into the mechanism of itch mediation and modulation. Hence, it became apparent that pruritus is a complex symptom or disease in itself, which requires particular attention to improve patients' health. Here, we summarize recent findings in pruritoceptive itch, including how this sensation is triggered and modulated by diverse endogenous and exogenous pruritogens and their receptors. A differentiation between mediating pruritogen and modulating pruritogen seems to be of great advantage to understand and decipher the molecular mechanism of itch perception. Only a comprehensive view on itch sensation will provide a solid basis for targeting this long-neglected adverse sensation accompanying numerous diseases and many drug side effects. Finally, we identify critical aspects of itch perception that require future investigation.

Auxiliary Master Equation for Nonequilibrium Dual-Fermion Approach.

We introduce an auxiliary quantum master equation dual fermion method and argue that it presents a convenient way to describe steady states of correlated impurity models. The scheme yields an expansion around a reference that is much closer to the true nonequilibrium state than that in the original dual fermion formulation. In steady-state situations, the scheme is numerically inexpensive and avoids time propagation. The Anderson impurity model is used to test the approach against numerically exact benchmarks.

Peri-pharyngeal muscle response to inspiratory loading: comparison of patients with OSA and healthy subjects.

Upper airway patency to airflow and the occurrence of obstructive sleep apnea involve a complex interplay between pharyngeal anatomy and synergic co-activation of peri-pharyngeal muscles. In previous studies we observed large differences in the response to sleep-associated flow limitation between the genioglossus and other (non-GG) peri-pharyngeal muscles. We hypothesized that similar differences are present also during wakefulness. In the present study we compared the response to inspiratory loading of the genioglossus electromyogram and four other peri-pharyngeal muscles. Studies were performed in eight obstructive sleep apnea patients, seven age-matched healthy subjects and five additional younger subjects. Electromyogram activity was evaluated over a range of negative oesophageal pressures and expressed as % of maximal electromyograms. In healthy subjects, the slope response to inspiratory loading (electromyogram/pressures) was similar for the genioglossus and non-GG muscles studied. However, the electromyogram responses were significantly higher in the young subjects compared with older subjects. In contrast, in the obstructive sleep apnea patients, the electromyogram/pressure response of the non-GG muscles was similar to that of the age-matched healthy subjects, whereas the slope response of the genioglossus electromyogram was significantly higher than non-GG muscles. We conclude that both age and the presence of obstructive sleep apnea affect the response of peri-pharyngeal muscles to inspiratory loading. In patients with obstructive sleep apnea the genioglossus seems to compensate for mechanical disadvantages, but non-GG muscles apparently are not included in this neuromuscular compensatory mechanism. Our current and previous findings suggest that attempts to improve obstructive sleep apnea with myofunctional therapy should put added emphasis on the training of non-GG muscles.

Lead geometry and transport statistics in molecular junctions.

We present a numerically exact study of charge transport and its fluctuations through a molecular junction driven out of equilibrium by a bias voltage, using the inchworm quantum Monte Carlo method. After showing how the technique can be used to address any lead geometry, we concentrate on one dimensional chains as an example. The finite bandwidth of the leads is shown to affect transport properties in ways that cannot be fully captured by quantum master equations: in particular, we reveal an interaction-induced broadening of transport channels that is visible at all voltages and show how fluctuations of the current are a more sensitive probe of this effect than the mean current.

Nrf2 Activation by SK-119 Attenuates Oxidative Stress, UVB, and LPS-Induced Damage.

BACKGROUND/AIMS: The Nrf2 signaling pathway plays a pivotal role in neutralizing excess reactive oxygen species formation and therefore enhancing the endogenous cellular protection mechanism. Thus, activating this pathway may provide therapeutic options against oxidative stress-related disorders. We have recently applied a computer-aided drug design approach to the design and synthesis of novel Nrf2 enhancers. The current study was aimed at investigating the potential beneficial impact of (E)-5-oxo-1-(4-((2,4,6-trihydroxybenzylidene)amino)phenyl)pyrrolidine-3-carboxylic acid (SK-119) in skin oxidative damage models. METHODS: SK-119, tested initially in PC-12 cells, attenuated oxidative stress-induced cytotoxicity concomitantly with Nrf2 activation. The potential impact of this compound was evaluated in skin-based disease models both in vitro (HaCaT cells) and ex vivo (human skin organ culture). RESULTS: The data clearly showed the marked anti-inflammatory and photoprotection properties of the compound; SK-119-treated cells or tissues displayed a reduction in cytokine secretion induced by lipopolysaccharides (LPS) in a manner comparable with dexamethasone. In addition, topical application of SK-119 was able to block UVB-induced oxidative stress and attenuated caspase-mediated apoptosis, DNA adduct formation, and the concomitant cellular damage. CONCLUSION: These results indicate that SK-119 is an Nrf2 activator that can be used as a prototype molecule for the development of novel treatments of dermatological disorders related to oxidative stress.

Regulation of GLUT4 activity in myotubes by 3-O-methyl-d-glucose.

The rate of glucose influx to skeletal muscles is determined primarily by the number of functional units of glucose transporter-4 (GLUT4) in the myotube plasma membrane. The abundance of GLUT4 in the plasma membrane is tightly regulated by insulin or contractile activity, which employ distinct pathways to translocate GLUT4-rich vesicles from intracellular compartments. Various studies have indicated that GLUT4 intrinsic activity is also regulated by conformational changes and/or interactions with membrane components and intracellular proteins in the vicinity of the plasma membrane. Here we show that the non-metabolizable glucose analog 3-O-methyl-d-glucose (MeGlc) augmented the rate of hexose transport into myotubes by increasing GLUT4 intrinsic activity without altering the content of the transporter in the plasma membrane. This effect was not a consequence of ATP depletion or hyperosmolar stress and did not involve Akt/PKB or AMPK signal transduction pathways. MeGlc reduced the inhibitory potency (increased Ki) of indinavir, a selective inhibitor of GLUT4, in a dose-dependent manner. Kinetic analyses indicate that MeGlc induced changes in GLUT4 or GLUT4 complexes within the plasma membrane, which enhanced the hexose transport activity and reduced the potency of indinavir inhibition. Finally, we present a simple kinetic analysis for screening and discovering low molecular weight compounds that augment GLUT4 activity.

A Novel Phenylchromane Derivative Increases the Rate of Glucose Uptake in L6 Myotubes and Augments Insulin Secretion from Pancreatic Beta-Cells by Activating AMPK.

PURPOSE: A series of novel polycyclic aromatic compounds that augment the rate of glucose uptake in L6 myotubes and increase glucose-stimulated insulin secretion from beta-cells were synthesized. Designing these molecules, we have aimed at the two main pathogenic mechanisms of T2D, deficient insulin secretion and diminished glucose clearance. The ultimate purpose of this work was to create a novel antidiabetic drug candidate with bi-functional mode of action. METHODS: All presented compounds were synthesized, and characterized in house. INS-1E cells and L6 myoblasts were used for the experiments. The rate of glucose uptake, mechanism of action, level of insulin secretion and the druggability of the lead compound were studied. RESULTS: The lead compound (6-(1,3-dithiepan-2-yl)-2-phenylchromane), dose- and time-dependently at the low µM range increased the rate of glucose uptake in L6 myotubes and insulin secretion in INS-1E cells. The compound exerted its effects through the activation of the LKB1 (Liver Kinase B1)-AMPK pathway. In vitro metabolic parameters of this lead compound exhibited good druggability. CONCLUSIONS: We anticipate that bi-functionality (increased rate of glucose uptake and augmented insulin secretion) will allow the lead compound to be a starting point for the development of a novel class of antidiabetic drugs.

Inchworm Monte Carlo for exact non-adiabatic dynamics. II. Benchmarks and comparison with established methods.

In this second paper of a two part series, we present extensive benchmark results for two different inchworm Monte Carlo expansions for the spin-boson model. Our results are compared to previously developed numerically exact approaches for this problem. A detailed discussion of convergence and error propagation is presented. Our results and analysis allow for an understanding of the benefits and drawbacks of inchworm Monte Carlo compared to other approaches for exact real-time non-adiabatic quantum dynamics.

Inchworm Monte Carlo for exact non-adiabatic dynamics. I. Theory and algorithms.

In this paper, we provide a detailed description of the inchworm Monte Carlo formalism for the exact study of real-time non-adiabatic dynamics. This method optimally recycles Monte Carlo information from earlier times to greatly suppress the dynamical sign problem. Using the example of the spin-boson model, we formulate the inchworm expansion in two distinct ways: The first with respect to an expansion in the system-bath coupling and the second as an expansion in the diabatic coupling. The latter approach motivates the development of a cumulant version of the inchworm Monte Carlo method, which has the benefit of improved scaling. This paper deals completely with methodology, while Paper II provides a comprehensive comparison of the performance of the inchworm Monte Carlo algorithms to other exact methodologies as well as a discussion of the relative advantages and disadvantages of each.

Beta cell response to nutrient overload involves phospholipid remodelling and lipid peroxidation.

AIMS/HYPOTHESIS: Membrane phospholipids are the major intracellular source for fatty acid-derived mediators, which regulate myriad cell functions. We showed previously that high glucose levels triggered the hydrolysis of polyunsaturated fatty acids from beta cell phospholipids. These fatty acids were subjected to free radical-catalysed peroxidation to generate the bioactive aldehyde 4-hydroxy-2E-nonenal (4-HNE). The latter activated the nuclear peroxisome proliferator-activated receptor-δ (PPARδ), which in turn augmented glucose-stimulated insulin secretion. The present study aimed at investigating the combined effects of glucose and fatty acid overload on phospholipid turnover and the subsequent generation of lipid mediators, which affect insulin secretion and beta cell viability. METHODS: INS-1E cells were incubated with increasing glucose concentrations (5-25 mmol/l) without or with palmitic acid (PA; 50-500 µmol/l) and taken for fatty acid-based lipidomic analysis and functional assays. Rat isolated islets of Langerhans were used similarly. RESULTS: PA was incorporated into membrane phospholipids in a concentration- and time-dependent manner; incorporation was highest at 25 mmol/l glucose. This was coupled to a rapid exchange with saturated, mono-unsaturated and polyunsaturated fatty acids. Importantly, released arachidonic acid and linoleic acid were subjected to peroxidation, resulting in the generation of 4-HNE, which further augmented insulin secretion by activating PPARδ in beta cells. However, this adaptive increase in insulin secretion was abolished at high glucose and PA levels, which induced endoplasmic reticulum stress, apoptosis and cell death. CONCLUSIONS/INTERPRETATION: These findings highlight a key role for phospholipid remodelling and fatty acid peroxidation in mediating adaptive and cytotoxic interactions induced by nutrient overload in beta cells.