Efficient many-body non-Markovian dynamics of organic polaritons

Funding: P.F.-W. acknowledges support from EPSRC (EP/T518062/1). B.W.L. and J.K. acknowledge support from EPSRC (EP/T014032/1).We show how to simulate a model of many molecules with both strong coupling to many vibrational modes and collective coupling to a single photon mode. We do this by combining process tensor matrix product operator methods with a mean-field approximation which reduces the dimension of the problem. We analyze the steady-state of the model under incoherent pumping to determine the dependence of the polariton lasing threshold on cavity detuning, light-matter coupling strength, and environmental temperature. Moreover, by measuring two-time correlations, we study quadratic fluctuations about the mean-field to calculate the photoluminescence spectrum. Our method enables one to simulate many-body systems with strong coupling to multiple environments, and to extract both static and dynamical properties.Publisher PDFPeer reviewe

Sublinear scaling in non-Markovian open quantum systems simulations

Funder: M.C. and E.M.G. acknowledge funding from EPSRC grant no. EP/T01377X/1. B.W.L. and J.K. were supported by EPSRC grant no. EP/T014032/1.While several numerical techniques are available for predicting the dynamics of non-Markovian open quantum systems, most struggle with simulations for very long memory and propagation times, e.g., due to superlinear scaling with the number of time steps n. Here, we introduce a numerically exact algorithm to calculate process tensors—compact representations of environmental influences—which provides a scaling advantage over previous algorithms by leveraging self-similarity of the tensor networks that represent the environment. It is applicable to environments with Gaussian statistics, such as for spin-boson-type open quantum systems. Based on a divide-and-conquer strategy, our approach requires only (n log n) singular value decompositions for environments with infinite memory. Where the memory can be truncated after nc time steps, a nominal scaling (nc log nc) is found, which is independent of n. This improved scaling is enabled by identifying process tensors with repeatable blocks. To demonstrate the power and utility of our approach, we provide three examples. (1) We calculate the fluorescence spectra of a quantum dot under both strong driving and strong dot-phonon couplings, a task requiring simulations over millions of time steps, which we are able to perform in minutes. (2) We efficiently find process tensors describing superradiance of multiple emitters. (3) We explore the limits of our algorithm by considering coherence decay with a very strongly coupled environment. The observed computation time is not necessarily proportional to the number of singular value decompositions because the matrix dimensions also depend on the number of time steps. Nevertheless, quasilinear and sublinear scaling of computation time is found in practice for a wide range of parameters. While there are instances where existing methods can achieve comparable nominal scaling by precalculating effective propagators for time-independent or periodic system Hamiltonians, process tensors contain all the information needed to extract arbitrary multitime correlation functions of the system when driven by arbitrary time-dependent system Hamiltonians. The algorithm we present here not only significantly extends the scope of numerically exact techniques to open quantum systems with long memory times, but it also has fundamental implications for the simulation complexity of tensor network approaches.Publisher PDFPeer reviewe

Sublinear scaling in non-Markovian open quantum systems simulations

While several numerical techniques are available for predicting the dynamics of non-Markovian open quantum systems, most struggle with simulations for very long memory and propagation times, e.g., due to superlinear scaling with the number of time steps $n$. Here, we introduce a numerically exact algorithm to calculate process tensors -- compact representations of environmental influences -- which provides a scaling advantage over previous algorithms by leveraging self-similarity of the tensor networks that represent Gaussian environments. Based on a divide-and-conquer strategy, our approach requires only $\mathcal{O}(n\log n)$ singular value decompositions for environments with infinite memory. Where the memory can be truncated after $n_c$ time steps, a scaling $\mathcal{O}(n_c\log n_c)$ is found, which is independent of $n$. This improved scaling is enabled by identifying process tensors with repeatable blocks. To demonstrate the power and utility of our approach we provide three examples. (1) We calculate the fluorescence spectra of a quantum dot under both strong driving and strong dot-phonon couplings, a task requiring simulations over millions of time steps, which we are able to perform in minutes. (2) We efficiently find process tensors describing superradiance of multiple emitters. (3) We explore the limits of our algorithm by considering coherence decay with a very strongly coupled environment. The algorithm we present here not only significantly extends the scope of numerically exact techniques to open quantum systems with long memory times, but also has fundamental implications for simulation complexity

Determining the validity of cumulant expansions for central spin models

For a model with many-to-one connectivity it is widely expected that mean-field theory captures the exact many-particle $N\to\infty$ limit, and that higher-order cumulant expansions of the Heisenberg equations converge to this same limit whilst providing improved approximations at finite $N$. Here we show that this is in fact not always the case. Instead, whether mean-field theory correctly describes the large-$N$ limit depends on how the model parameters scale with $N$, and we show that convergence of cumulant expansions may be non-uniform across even and odd orders. Further, even when a higher-order cumulant expansion does recover the correct limit, the error is not monotonic with $N$ and may exceed that of mean-field theory.Comment: 7 pages, 3 figures plus supplementary materia

Few-emitter lasing in single ultra-small nanocavities

Funding: We acknowledge support from EPSRC grants EP/G060649/1, EP/L027151/1, EP/G037221/1, EP/T014032/1, EPSRC NanoDTC, and from the European Research Council (ERC) under Horizon 2020 research and innovation programme PICOFORCE (Grant Agreement No. 883703), THOR (Grant Agreement No. 829067) and POSEIDON (Grant Agreement No. 861950). O.S.O acknowledges the support of a Rubicon fellowship from the Netherlands Organisation for Scientific Research.Lasers are ubiquitous for information storage, processing, communications, sensing, biological research, and medical applications. To decrease their energy and materials usage, a key quest is to miniaturize lasers down to nanocavities. Obtaining the smallest mode volumes demands plasmonic nanocavities, but for these, gain comes from only single or few emitters. Until now, lasing in such devices was unobtainable due to low gain and high cavity losses. Here, we demonstrate a form of “few emitter lasing” in a plasmonic nanocavity approaching the single-molecule emitter regime. The few-emitter lasing transition significantly broadens, and depends on the number of molecules and their individual locations. We show this non-standard few-emitter lasing can be understood by developing a theoretical approach extending previous weak-coupling theories. Our work paves the way for developing nanolaser applications as well as fundamental studies at the limit of few emitters.Publisher PDFPeer reviewe

Preliminary Findings from an Interventional Study using Network Analysis to Support Management in Local Health Departments in Florida

Management is the core service that integrates and coordinates essential public health services. Managers of local health departments (LHDs) are experts in practice but may not have expertise in organizational management. We conducted an evidence-based training intervention in 10 LHDs in Florida to support managers’ decision-making on organizational integration and coordination. We deployed a standard survey to collect organizational network measurements pre and post intervention. We presented results as evidence-based performance feedback and interviewed managers to document how they used the results in the context of each organization. Post intervention we found unexpected, significantly higher network centralization in daily work. We attributed this increase in hierarchical communication to preparations for a statewide accreditation initiative. When QI initiatives are undertaken globally within a state, managers and leaders need to be alert for possible impact on autonomous decision-making of professionals at the point of service which could affect service delivery

Optimizing performance of quantum operations with non-Markovian decoherence: the tortoise or the hare?

The interaction between a quantum system and its environment limits our ability to control it and perform quantum operations on it. We present an efficient method to find optimal controls for quantum systems coupled to non-Markovian environments, by using the process tensor to compute the gradient of an objective function. We consider state transfer for a driven two-level system coupled to a bosonic environment, and characterize performance in terms of speed and fidelity. We thus determine the best achievable fidelity as a function of process duration. We show there is a trade-off between speed and fidelity, and that slower processes can have higher fidelity by exploiting non-Markovian effects.Comment: 6 pages, 3 figure

Modelling the impact of local reactive school closures on critical care provision during an influenza pandemic

Despite the fact that the 2009 H1N1 pandemic influenza strain was less severe than had been feared, both seasonal epidemics of influenza-like-illness and future influenza pandemics have the potential to place a serious burden on health services. The closure of schools has been postulated as a means of reducing transmission between children and hence reducing the number of cases at the peak of an epidemic; this is supported by the marked reduction in cases during school holidays observed across the world during the 2009 pandemic. However, a national policy of long-duration school closures could have severe economic costs. Reactive short-duration closure of schools in regions where health services are close to capacity offers a potential compromise, but it is unclear over what spatial scale and time frame closures would need to be made to be effective. Here, using detailed geographical information for England, we assess how localized school closures could alleviate the burden on hospital intensive care units (ICUs) that are reaching capacity. We show that, for a range of epidemiologically plausible assumptions, considerable local coordination of school closures is needed to achieve a substantial reduction in the number of hospitals where capacity is exceeded at the peak of the epidemic. The heterogeneity in demand per hospital ICU bed means that even widespread school closures are unlikely to have an impact on whether demand will exceed capacity for many hospitals. These results support the UK decision not to use localized school closures as a control mechanism, but have far wider international public-health implications. The spatial heterogeneities in both population density and hospital capacity that give rise to our results exist in many developed countries, while our model assumptions are sufficiently general to cover a wide range of pathogens. This leads us to believe that when a pandemic has severe implications for ICU capacity, only widespread school closures (with their associated costs and organizational challenges) are sufficient to mitigate the burden on the worst-affected hospitals

Does egg deposition by herbivorous pine sawflies affect transcription of sesquiterpene synthases in pine?

Scots pine (Pinus sylvestris; Pinaceae, Pinales) is known to defend against egg deposition by herbivorous sawflies by changing its terpenoid volatile blend. The oviposition-induced pine odor attracts egg parasitoids that kill the sawfly eggs. Here, we investigated whether sawfly egg deposition activates genes encoding pine terpene synthases by extracting mRNA from oviposition-induced P. sylvestris. Three new sesquiterpene synthases, PsTPS 1, PsTPS 2, and PsTPS 3, were isolated that were shown on heterologous expression in Escherichia coli to produce (E)-β-caryophyllene and α-humulene (PsTPS 1), 1(10),5-germacradiene-4-ol (PsTPS 2), and longifolene and α-longipinene (PsTPS 3) as their principal products. Quantitative RT-PCR analyses revealed that transcript levels of PsTPS 1 and PsTPS 2 were significantly higher in oviposition-induced twigs that were attractive to the parasitoids than in non-attractive, artificially damaged twigs. Thus, our results demonstrate a specific transcription response to egg deposition, distinct from that caused by artificial wounding. Transcripts of PsTPS 3 did not change in response to egg deposition. The transcript levels of PsTPS 1, PsTPS 2, and PsTPS 3 were also determined in relation to time after egg deposition, since pine odor is attractive to the parasitoid only 72 h after egg deposition. Transcription rates of PsTPS 1 and PsTPS 2 were significantly enhanced only 72 h after egg deposition, thus matching the timing of odor attractiveness, while for PsTPS 3, enhanced transcription was not detected at any time period studied after egg deposition. The ecological significance of the oviposition-induced increase of sesquiterpene synthase transcripts is discussed

Increased Terpenoid Accumulation in Cotton (Gossypium hirsutum) Foliage is a General Wound Response

The subepidermal pigment glands of cotton accumulate a variety of terpenoid products, including monoterpenes, sesquiterpenes, and terpenoid aldehydes that can act as feeding deterrents against a number of insect herbivore species. We compared the effect of herbivory by Spodoptera littoralis caterpillars, mechanical damage by a fabric pattern wheel, and the application of jasmonic acid on levels of the major representatives of the three structural classes of terpenoids in the leaf foliage of 4-week-old Gossypium hirsutum plants. Terpenoid levels increased successively from control to mechanical damage, herbivory, and jasmonic acid treatments, with E-β-ocimene and heliocide H1 and H4 showing the highest increases, up to 15-fold. Herbivory or mechanical damage to older leaves led to terpenoid increases in younger leaves. Leaf-by-leaf analysis of terpenes and gland density revealed that higher levels of terpenoids were achieved by two mechanisms: (1) increased filling of existing glands with terpenoids and (2) the production of additional glands, which were found to be dependent on damage intensity. As the relative response of individual terpenoids did not differ substantially among herbivore, mechanical damage, and jasmonic acid treatments, the induction of terpenoids in cotton foliage appears to represent a non-specific wound response mediated by jasmonic acid