Observation of Rabi dynamics with a short-wavelength free-electron laser.

Rabi oscillations are periodic modulations of populations in two-level systems interacting with a time-varying field1. They are ubiquitous in physics with applications in different areas such as photonics2, nano-electronics3, electron microscopy4 and quantum information5. While the theory developed by Rabi was intended for fermions in gyrating magnetic fields, Autler and Townes realized that it could also be used to describe coherent light-matter interactions within the rotating-wave approximation6. Although intense nanometre-wavelength light sources have been available for more than a decade7-9, Rabi dynamics at such short wavelengths has not been directly observed. Here we show that femtosecond extreme-ultraviolet pulses from a seeded free-electron laser10 can drive Rabi dynamics between the ground state and an excited state in helium atoms. The measured photoelectron signal reveals an Autler-Townes doublet and an avoided crossing, phenomena that are both fundamental to coherent atom-field interactions11. Using an analytical model derived from perturbation theory on top of the Rabi model, we find that the ultrafast build-up of the doublet structure carries the signature of a quantum interference effect between resonant and non-resonant photoionization pathways. Given the recent availability of intense attosecond12 and few-femtosecond13 extreme-ultraviolet pulses, our results unfold opportunities to carry out ultrafast manipulation of coherent processes at short wavelengths using free-electron lasers.

Phylogenetic inference of where species spread or split across barriers.

SignificanceGeography molds how species evolve in space. Strong geographical barriers to movement, for instance, both inhibit dispersal between regions and allow isolated populations to diverge as new species. Weak barriers, by contrast, permit species range expansion and persistence. These factors present a conundrum: How strong must a barrier be before between-region speciation outpaces dispersal? We designed a phylogenetic model of dispersal, extinction, and speciation that allows regional features to influence rates of biogeographic change and applied it to the neotropical radiation of Anolis lizards. Separation by water induces a threefold steeper barrier to movement than equivalent distances over land. Our model will help biologists detect relationships between evolutionary processes and the spatial contexts in which they operate.

Relict inland mangrove ecosystem reveals Last Interglacial sea levels.

Climatic oscillations during the Pleistocene played a major role in shaping the spatial distribution and demographic dynamics of Earth's biota, including our own species. The Last Interglacial (LIG) or Eemian Period (ca. 130 to 115 thousand years B.P.) was particularly influential because this period of peak warmth led to the retreat of all ice sheets with concomitant changes in global sea level. The impact of these strong environmental changes on the spatial distribution of marine and terrestrial ecosystems was severe as revealed by fossil data and paleogeographic modeling. Here, we report the occurrence of an extant, inland mangrove ecosystem and demonstrate that it is a relict of the LIG. This ecosystem is currently confined to the banks of the freshwater San Pedro Mártir River in the interior of the Mexico-Guatemala El Petén rainforests, 170 km away from the nearest ocean coast but showing the plant composition and physiognomy typical of a coastal lagoon ecosystem. Integrating genomic, geologic, and floristic data with sea level modeling, we present evidence that this inland ecosystem reached its current location during the LIG and has persisted there in isolation ever since the oceans receded during the Wisconsin glaciation. Our study provides a snapshot of the Pleistocene peak warmth and reveals biotic evidence that sea levels substantially influenced landscapes and species ranges in the tropics during this period.

Evolution of Gene Expression across Species and Specialized Zooids in Siphonophora.

Siphonophores are complex colonial animals, consisting of asexually produced bodies (zooids) that are functionally specialized for specific tasks, including feeding, swimming, and sexual reproduction. Though this extreme functional specialization has captivated biologists for generations, its genomic underpinnings remain unknown. We use RNA-seq to investigate gene expression patterns in five zooids and one specialized tissue across seven siphonophore species. Analyses of gene expression across species present several challenges, including identification of comparable expression changes on gene trees with complex histories of speciation, duplication, and loss. We examine gene expression within species, conduct classical analyses examining expression patterns between species, and introduce species branch filtering, which allows us to examine the evolution of expression across species in a phylogenetic framework. Within and across species, we identified hundreds of zooid-specific and species-specific genes, as well as a number of putative transcription factors showing differential expression in particular zooids and developmental stages. We found that gene expression patterns tended to be largely consistent in zooids with the same function across species, but also some large lineage-specific shifts in gene expression. Our findings show that patterns of gene expression have the potential to define zooids in colonial organisms. Traditional analyses of the evolution of gene expression focus on the tips of gene phylogenies, identifying large-scale expression patterns that are zooid or species variable. The new explicit phylogenetic approach we propose here focuses on branches (not tips) offering a deeper evolutionary perspective into specific changes in gene expression within zooids along all branches of the gene (and species) trees.

Fast Likelihood Calculations for Automatic Identification of Macroevolutionary Rate Heterogeneity in Continuous and Discrete Traits.

Understanding phenotypic disparity across the tree of life requires identifying where and when evolutionary rates change on phylogeny. A primary methodological challenge in macroevolution is therefore to develop methods for accurate inference of among-lineage variation in rates of phenotypic evolution. Here, we describe a method for inferring among-lineage evolutionary rate heterogeneity in both continuous and discrete traits. The method assumes that the present-day distribution of a trait is shaped by a variable-rate process arising from a mixture of constant-rate processes and uses a single-pass tree traversal algorithm to estimate branch-specific evolutionary rates. By employing dynamic programming optimization techniques and approximate maximum likelihood estimators where appropriate, our method permits rapid exploration of the tempo and mode of phenotypic evolution. Simulations indicate that the method reconstructs rates of trait evolution with high accuracy. Application of the method to data sets on squamate reptile reproduction and turtle body size recovers patterns of rate heterogeneity identified by previous studies but with computational costs reduced by many orders of magnitude. Our results expand the set of tools available for detecting macroevolutionary rate heterogeneity and point to the utility of fast, approximate methods for studying large-scale biodiversity dynamics. [Brownian motion; continuous characters; discrete characters; macroevolution; Markov process; rate heterogeneity.].

Natural selection maintains species despite frequent hybridization in the desert shrub Encelia.

Natural selection is an important driver of genetic and phenotypic differentiation between species. For species in which potential gene flow is high but realized gene flow is low, adaptation via natural selection may be a particularly important force maintaining species. For a recent radiation of New World desert shrubs (Encelia: Asteraceae), we use fine-scale geographic sampling and population genomics to determine patterns of gene flow across two hybrid zones formed between two independent pairs of species with parapatric distributions. After finding evidence for extremely strong selection at both hybrid zones, we use a combination of field experiments, high-resolution imaging, and physiological measurements to determine the ecological basis for selection at one of the hybrid zones. Our results identify multiple ecological mechanisms of selection (drought, salinity, herbivory, and burial) that together are sufficient to maintain species boundaries despite high rates of hybridization. Given that multiple pairs of Encelia species hybridize at ecologically divergent parapatric boundaries, such mechanisms may maintain species boundaries throughout Encelia.

Reference genome of the color polymorphic desert annual plant sandblossoms, Linanthus parryae.

Sandblossoms, Linanthus parryae is a widespread annual plant species found in washes and sandy open habitats across the Mojave Desert and Eastern Sierra Nevada of California. Studies in this species have played a central role in evolutionary biology, serving as the first test cases of the shifting balance theory of evolution, models of isolation by distance, and metrics to describe the genetic structure of natural populations. Despite the importance of L. parryae in the development of landscape genetics and phylogeography, there are no genomic resources available for the species. Through the California Conservation Genomics Project, we assembled the first genome in the genus Linanthus. Using PacBio HiFi long reads and Hi-C chromatin conformation capture, we assembled 123 scaffolds spanning 1.51 Gb of the 1.96 Gb estimated genome, with a contig N50 of 18.7 Mb and a scaffold N50 of 124.8 Mb. This assembly, with a BUSCO completeness score of 88.7%, will allow us to revisit foundational ideas central to our understanding of how evolutionary forces operate in a geographic landscape. In addition, it will be a new resource to uncover adaptations to arid environments in the fragile desert habitat threatened by urban and solar farm development, climate change, and off-road vehicles.

Photoionization and core resonances from range-separated density-functional theory: General formalism and example of the beryllium atom.

We explore the merits of linear-response range-separated time-dependent density-functional theory (TDDFT) for the calculation of photoionization spectra. We consider two variants of range-separated TDDFT, namely, the time-dependent range-separated hybrid (TDRSH) scheme, which uses a global range-separation parameter, and the time-dependent locally range-separated hybrid (TDLRSH), which uses a local range-separation parameter, and compare with standard time-dependent local-density approximation (TDLDA) and time-dependent Hartree-Fock (TDHF). We show how to calculate photoionization spectra with these methods using the Sternheimer approach formulated in a non-orthogonal B-spline basis set with appropriate frequency-dependent boundary conditions. We illustrate these methods on the photoionization spectrum of the Be atom, focusing, in particular, on the core resonances. Both the TDRSH and TDLRSH photoionization spectra are found to constitute a large improvement over the TDLDA photoionization spectrum and a more modest improvement over the TDHF photoionization spectrum.

Photoionization and core resonances from range-separated time-dependent density-functional theory for open-shell states: Example of the lithium atom.

We consider the calculations of photoionization spectra and core resonances of open-shell systems using range-separated time-dependent density-functional theory. Specifically, we use the time-dependent range-separated hybrid (TDRSH) scheme, combining a long-range Hartree-Fock exchange potential and kernel with a short-range potential and kernel from a local density-functional approximation, and the time-dependent locally range-separated hybrid (TDLRSH) scheme, which uses a local range-separation parameter. To efficiently perform the calculations, we formulate a spin-unrestricted linear-response Sternheimer approach in a non-orthogonal B-spline basis set using appropriate frequency-dependent boundary conditions. We illustrate this approach on the Li atom, which suggests that TDRSH and TDLRSH are adequate simple methods for estimating the single-electron photoionization spectra of open-shell systems.

Diversification, disparification and hybridization in the desert shrubs Encelia.

There are multiple hypotheses for the spectacular plant diversity found in deserts. We explore how different factors, including the roles of ecological opportunity and selection, promote diversification and disparification in Encelia, a lineage of woody plants in the deserts of the Americas. Using a nearly complete species-level phylogeny based on double-digest restriction-aided sequencing along with a broad set of phenotypic traits, we estimate divergence times and diversification rates, identify instances of hybridization, quantify trait disparity and assess phenotypic divergence across environmental gradients. We show that Encelia originated and diversified recently (mid-Pleistocene) and rapidly, with rates comparable to notable adaptive radiations in plants. Encelia probably originated in the hot deserts of North America, with subsequent diversification across steep environmental gradients. We uncover multiple instances of gene flow between species. The radiation of Encelia is characterized by fast rates of phenotypic evolution, trait lability and extreme disparity across environments and between species pairs with overlapping geographic ranges. Encelia exemplifies how interspecific gene flow in combination with high trait lability can enable exceptionally fast diversification and disparification across steep environmental gradients.

Pairwise comparisons across species are problematic when analyzing functional genomic data.

There is considerable interest in comparing functional genomic data across species. One goal of such work is to provide an integrated understanding of genome and phenotype evolution. Most comparative functional genomic studies have relied on multiple pairwise comparisons between species, an approach that does not incorporate information about the evolutionary relationships among species. The statistical problems that arise from not considering these relationships can lead pairwise approaches to the wrong conclusions and are a missed opportunity to learn about biology that can only be understood in an explicit phylogenetic context. Here, we examine two recently published studies that compare gene expression across species with pairwise methods, and find reason to question the original conclusions of both. One study interpreted pairwise comparisons of gene expression as support for the ortholog conjecture, the hypothesis that orthologs tend to have more similar attributes (expression in this case) than paralogs. The other study interpreted pairwise comparisons of embryonic gene expression across distantly related animals as evidence for a distinct evolutionary process that gave rise to phyla. In each study, distinct patterns of pairwise similarity among species were originally interpreted as evidence of particular evolutionary processes, but instead, we find that they reflect species relationships. These reanalyses concretely show the inadequacy of pairwise comparisons for analyzing functional genomic data across species. It will be critical to adopt phylogenetic comparative methods in future functional genomic work. Fortunately, phylogenetic comparative biology is also a rapidly advancing field with many methods that can be directly applied to functional genomic data.

Phonon-Mediated and Weakly Size-Dependent Electron and Hole Cooling in CsPbBr3 Nanocrystals Revealed by Atomistic Simulations and Ultrafast Spectroscopy.

We combine state-of-the-art ultrafast photoluminescence and absorption spectroscopy and nonadiabatic molecular dynamics simulations to investigate charge-carrier cooling in CsPbBr3 nanocrystals over a very broad size regime, from 0.8 to 12 nm. Contrary to the prevailing notion that polaron formation slows down charge-carrier cooling in lead-halide perovskites, no suppression of carrier cooling is observed in CsPbBr3 nanocrystals except for a slow cooling (over ∼10 ps) of "warm" electrons in the vicinity (within ∼0.1 eV) of the conduction band edge. At higher excess energies, electrons and holes cool with similar rates, on the order of 1 eV ps-1 carrier-1, increasing weakly with size. Our ab initio simulations suggest that cooling proceeds via fast phonon-mediated intraband transitions driven by strong and size-dependent electron-phonon coupling. The presented experimental and computational methods yield the spectrum of involved phonons and may guide the development of devices utilizing hot charge carriers.

Continuous-wave infrared optical gain and amplified spontaneous emission at ultralow threshold by colloidal HgTe quantum dots.

Colloidal quantum dots (QDs) raise more and more interest as solution-processable and tunable optical gain materials. However, especially for infrared active QDs, optical gain remains inefficient. Since stimulated emission involves multifold degenerate band-edge states, population inversion can be attained only at high pump power and must compete with efficient multi-exciton recombination. Here, we show that mercury telluride (HgTe) QDs exhibit size-tunable stimulated emission throughout the near-infrared telecom window at thresholds unmatched by any QD studied before. We attribute this unique behaviour to surface-localized states in the bandgap that turn HgTe QDs into 4-level systems. The resulting long-lived population inversion induces amplified spontaneous emission under continuous-wave optical pumping at power levels compatible with solar irradiation and direct current electrical pumping. These results introduce an alternative approach for low-threshold QD-based gain media based on intentional trap states that paves the way for solution-processed infrared QD lasers and amplifiers.

Issues and Perspectives in Species Delimitation using Phenotypic Data: Atlantean Evolution in Darwin's Finches.

Progress in the development and use of methods for species delimitation employing phenotypic data lags behind conceptual and practical advances in molecular genetic approaches. The basic evolutionary model underlying the use of phenotypic data to delimit species assumes random mating and quantitative polygenic traits, so that phenotypic distributions within a species should be approximately normal for individuals of the same sex and age. Accordingly, two or more distinct normal distributions of phenotypic traits suggest the existence of multiple species. In light of this model, we show that analytical approaches employed in taxonomic studies using phenotypic data are often compromised by three issues: 1) reliance on graphical analyses that convey little information on phenotype frequencies; 2) exclusion of characters potentially important for species delimitation following reduction of data dimensionality; and 3) use of measures of central tendency to evaluate phenotypic distinctiveness. We outline approaches to overcome these issues based on statistical developments related to normal mixture models (NMMs) and illustrate them empirically with a reanalysis of morphological data recently used to claim that there are no morphologically distinct species of Darwin's ground-finches (Geospiza). We found negligible support for this claim relative to taxonomic hypotheses recognizing multiple species. Although species limits among ground-finches merit further assessments using additional sources of information, our results bear implications for other areas of inquiry including speciation research: because ground-finches have likely speciated and are not trapped in a process of "Sisyphean" evolution as recently argued, they remain useful models to understand the evolutionary forces involved in speciation. Our work underscores the importance of statistical approaches grounded on appropriate evolutionary models for species delimitation. We discuss how NMMs offer new perspectives in the kind of inferences available to systematists, with significant repercussions on ideas about the phenotypic structure of biodiversity.

QMflows: A Tool Kit for Interoperable Parallel Workflows in Quantum Chemistry.

We present the QMflows Python package for quantum chemistry workflow automatization. QMflows allows users to write complex workflows in terms of simple Python scripts. It supports the development of interoperable workflows involving multiple quantum chemistry codes and executes them efficiently on large scale parallel computers. This open source library provides standardized interfaces to a number of quantum chemistry packages and can be easily extended to accommodate additional codes. QMflows features are described and illustrated with a number of representative applications.

Linear-response range-separated density-functional theory for atomic photoexcitation and photoionization spectra.

We investigate the performance of the range-separated hybrid (RSH) scheme, which combines long-range Hartree-Fock (HF) and a short-range density-functional approximation (DFA), for calculating the photoexcitation/photoionization spectra of the H and He atoms, using a B-spline basis set in order to correctly describe the continuum part of the spectra. The study of these simple systems allows us to quantify the influence on the spectra of the errors coming from the short-range exchange-correlation DFA and from the missing long-range correlation in the RSH scheme. We study the differences using the long-range HF exchange (nonlocal) potential and the long-range exact exchange (local) potential. Contrary to the former, the latter supports a series of Rydberg states and gives reasonable photoexcitation/photoionization spectra, even without applying linear-response theory. The most accurate spectra are obtained with the linear-response time-dependent RSH (TDRSH) scheme. In particular, for the He atom at the optimal value of the range-separation parameter, TDRSH gives slightly more accurate photoexcitation and photoionization spectra than the standard linear-response time-dependent HF. More generally, this work shows the potential of range-separated density-functional theory for calculating linear and nonlinear optical properties involving continuum states.

Characterization of Multipath Effects in Indoor Positioning Systems by AoA and PoA Based on Optical Signals.

In this paper, we characterize and measure the effects of the errors introduced by the multipath when obtaining the position of an agent by means of Indoor Positioning Systems (IPS) based on optical signal. These effects are characterized in Local Positioning Systems (LPSs) based on two different techniques: the first one by determining the Angle of Arrival (AoA) of the infrared signal (IR) to the detector; and the second one by working with the measurement of the Phase shift of signal Arrival from the transmitter to a receiver (PoA). We present the obtained results and conclusions, which indicate that using Position Sensitive Devices (PSD) the multipath effects for AoA have little impact on the measurement, while for PoA the positioning errors are very significant, making the system useless in many cases.

Improved phylogenetic resolution within Siphonophora (Cnidaria) with implications for trait evolution.

Siphonophores are a diverse group of hydrozoans (Cnidaria) that are found at most depths of the ocean - from the surface, like the familiar Portuguese man of war, to the deep sea. They play important roles in ocean ecosystems, and are among the most abundant gelatinous predators. A previous phylogenetic study based on two ribosomal RNA genes provided insight into the internal relationships between major siphonophore groups. There was, however, little support for many deep relationships within the clade Codonophora. Here, we present a new siphonophore phylogeny based on new transcriptome data from 29 siphonophore species analyzed in combination with 14 publicly available genomic and transcriptomic datasets. We use this new phylogeny to reconstruct several traits that are central to siphonophore biology, including sexual system (monoecy vs. dioecy), gain and loss of zooid types, life history traits, and habitat. The phylogenetic relationships in this study are largely consistent with the previous phylogeny, but we find strong support for new clades within Codonophora that were previously unresolved. These results have important implications for trait evolution within Siphonophora, including favoring the hypothesis that monoecy arose at least twice.

Optomechanical Control of Quantum Yield in Trans-Cis Ultrafast Photoisomerization of a Retinal Chromophore Model.

The quantum yield of a photochemical reaction is one of the most fundamental quantities in photochemistry, as it measures the efficiency of the transduction of light energy into chemical energy. Nature has evolved photoreceptors in which the reactivity of a chromophore is enhanced by its molecular environment to achieve high quantum yields. The retinal chromophore sterically constrained inside rhodopsin proteins represents an outstanding example of such a control. In a more general framework, mechanical forces acting on a molecular system can strongly modify its reactivity. Herein, we show that the exertion of tensile forces on a simplified retinal chromophore model provokes a substantial and regular increase in the trans-to-cis photoisomerization quantum yield in a counterintuitive way, as these extension forces facilitate the formation of the more compressed cis photoisomer. A rationale for the mechanochemical effect on this photoisomerization mechanism is also proposed.

Molcas 8: New capabilities for multiconfigurational quantum chemical calculations across the periodic table.

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas-Kroll-Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC-PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large-scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.