Nonadiabatic quantum dynamics explores non-monotonic photodissociation branching of N2 into the N(4S) + N(2D) and N(4S) + N(2P) product channels.

Vacuum ultraviolet (VUV) photodissociation of N2 molecules is a source of reactive N atoms in the interstellar medium. In the energy range of VUV optical excitation of N2, the N-N triple bond cleavage leads to three types of atoms: ground-state N(4S) and excited-state N(2P) and N(2D). The latter is the highest reactive and it is believed to be the primary participant in reactions with hydrocarbons in Titan's atmosphere. Experimental studies have observed a non-monotonic energy dependence and non-statistical character of the photodissociation of N2. This implies different dissociation pathways and final atomic products for different wavelength regions in the sunlight spectrum. We here apply ab initio quantum chemical and nonadiabatic quantum dynamical techniques to follow the path of an electronic state from the excitation of a particular singlet 1Σ+u and 1Πu vibronic level of N2 to its dissociation into different atomic products. We simulate dynamics for two isotopomers of the nitrogen molecule, 14N2 and 14N15N for which experimental data on the branching are available. Our computations capture the non-monotonic energy dependence of the photodissociation branching ratios in the energy range 108 000-116 000 cm-1. Tracing the quantum dynamics in a bunch of electronic states enables us to identify the key components that determine the efficacy of singlet to triplet population transfer and therefore predissociation lifetimes and branching ratios for different energy regions.

Electronic Coherences Excited by an Ultra Short Pulse Are Robust with Respect to Averaging over Randomly Oriented Molecules as Shown by Singular Value Decomposition.

We report a methodology for averaging quantum photoexcitation vibronic dynamics over the initial orientations of the molecules with respect to an ultrashort light pulse. We use singular value decomposition of the ensemble density matrix of the excited molecules, which allows the identification of the few dominant principal molecular orientations with respect to the polarization direction of the electric field. The principal orientations provide insights into the specific stereodynamics of the corresponding principal molecular vibronic states. The massive compaction of the vibronic density matrix of the ensemble of randomly oriented pumped molecules enables a most efficient fully quantum mechanical time propagation scheme. Two examples are discussed for the quantum dynamics of the LiH molecule in the manifolds of its electronically excited Σ and Π states. Our results show that electronic and vibrational coherences between excited states of the same symmetry are resilient to averaging over an ensemble of molecular orientations and can be selectively excited at the ensemble level by tuning the pulse parameters.

Nonadiabatic dynamics in a forest of coupled states: Electronic state branching in the VUV photodissociation of N2.

Multi-state electronic dynamics at higher excitation energies is needed for the understanding of a variety of energy rich situations, including chemistry under extreme conditions, vacuum ultraviolet (VUV) induced astrochemistry, and attochemistry. It calls for an understanding of three stages, energy acquisition, dynamical propagation, and disposal. It is typically not possible to identify a basis of uncoupled quantum states that is sufficient for the three stages. The handicap is the large number of coupled quantum states that is needed to describe the system. Progress in quantum chemistry provides the necessary background to the energetics and the coupling. Progress in quantum dynamics takes this as input for the propagation in time. Right now, it seems that we have come of age with potential detailed applications. We here report a demonstration to a coupled electron-nuclear quantum dynamics through a maze of 47 electronic states and with attention to the order in perturbation theory that is indicated using propensity rules for couplings. Close agreement with experimental results for the VUV photodissociation of 14N2 and its isotopomer 14N15N is achieved. We pay special attention to the coupling between two dissociative continua and an optically accessible bound domain. The computations reproduce and interpret the non-monotonic branching between the two exit channels producing N(2D) and N(2P) atoms as a function of excitation energy and its variation with the mass.

Thermodynamic energetics underlying genomic instability and whole-genome doubling in cancer.

Genomic instability contributes to tumorigenesis through the amplification and deletion of cancer driver genes. DNA copy number (CN) profiling of ensembles of tumors allows a thermodynamic analysis of the profile for each tumor. The free energy of the distribution of CNs is found to be a monotonically increasing function of the average chromosomal ploidy. The dependence is universal across several cancer types. Surprisal analysis distinguishes two main known subgroups: tumors with cells that have or have not undergone whole-genome duplication (WGD). The analysis uncovers that CN states having a narrower distribution are energetically more favorable toward the WGD transition. Surprisal analysis also determines the deviations from a fully stable-state distribution. These deviations reflect constraints imposed by tumor fitness selection pressures. The results point to CN changes that are more common in high-ploidy tumors and thus support altered selection pressures upon WGD.

Intracellular redox potential is correlated with miRNA expression in MCF7 cells under hypoxic conditions.

Hypoxia is a ubiquitous feature of cancers, encouraging glycolytic metabolism, proliferation, and resistance to therapy. Nonetheless, hypoxia is a poorly defined term with confounding features described in the literature. Redox biology provides an important link between the external cellular microenvironment and the cell's response to changing oxygen pressures. In this paper, we demonstrate a correlation between intracellular redox potential (measured using optical nanosensors) and the concentrations of microRNAs (miRNAs) involved in the cell's response to changes in oxygen pressure. The correlations were established using surprisal analysis (an approach derived from thermodynamics and information theory). We found that measured redox potential changes reflect changes in the free energy computed by surprisal analysis of miRNAs. Furthermore, surprisal analysis identified groups of miRNAs, functionally related to changes in proliferation and metastatic potential that played the most significant role in the cell's response to changing oxygen pressure.

Consecutive feedback-driven constitutional dynamic networks.

Cellular transformations are driven by environmentally triggered complex dynamic networks, which include signal-triggered feedback processes, cascaded reactions, and switchable transformations. We apply the structural and functional information encoded in the sequences of nucleic acids to construct signal-triggered constitutional dynamic networks (CDNs) that mimic the functions of natural networks. Using predesigned hairpin structures as triggers, the network generates functional strands, which stabilize one or the other of the constituents of the network, leading to feedback-driven reconfiguration and time-dependent equilibration of the networks. Using structurally designed hairpins, positive-feedback or negative-feedback mechanisms operated by the CDNs are demonstrated. With two predesigned hairpins, the coupled consecutive operations of negative/positive- or positive/positive- feedback cascades are accomplished. The time-dependent composition changes of the networks are well reproduced by chemical kinetics simulations that provide predictive behaviors of the network, under variable auxiliary conditions. Beyond mimicking natural network properties and functions by means of the synthetic nucleic-acid-based CDNs, the systems introduce versatile perspectives for the design of amplified sensors (sensing of miRNA-376a) and the development of logic gate circuits.

Ellipticity controlled dissociative double ionization of ethane by strong fields.

The yields of all dissociation channels of ethane dications produced by strong field double ionization were measured. It was found that the branching ratios can be controlled by varying the ellipticity of laser pulses. The CH3+ formation and H+ formation channels show a clear competition, producing the highest and lowest branching ratios at ellipticity of ∼0.6, respectively. With the help of theoretical calculations, such a control was attributed to the ellipticity dependent yields of different sequential ionization pathways.

Identification of an ultrafast internal conversion pathway of pyrazine by time-resolved vacuum ultraviolet photoelectron spectrum simulations.

The internal conversion from the optically bright S2 (1B2u, ππ*) state to the dark S1 (1B3u, nπ*) state in pyrazine is a standard benchmark for experimental and theoretical studies on ultrafast radiationless decay. Since 2008, a few theoretical groups have suggested significant contributions of other dark states S3 (1Au, nπ*) and S4 (1B2g, nπ*) to the decay of S2. We have previously reported the results of nuclear wave packet simulations [Kanno et al., Phys. Chem. Chem. Phys. 17, 2012 (2015)] and photoelectron spectrum calculations [Mignolet et al., Chem. Phys. 515, 704 (2018)] that support the conventional two-state picture. In this article, the two different approaches, i.e., wave packet simulation and photoelectron spectrum calculation, are combined: We computed the time-resolved vacuum ultraviolet photoelectron spectrum and photoelectron angular distribution for the ionization of the wave packet transferred from S2 to S1. The present results reproduce almost all the characteristic features of the corresponding experimental time-resolved spectrum [Horio et al., J. Chem. Phys. 145, 044306 (2016)], such as a rapid change from a three-band to two-band structure. This further supports the existence and character of the widely accepted pathway (S2 → S1) of ultrafast internal conversion in pyrazine.

Compacting the density matrix in quantum dynamics: Singular value decomposition of the surprisal and the dominant constraints for anharmonic systems.

We introduce a practical method for compacting the time evolution of the quantum state of a closed physical system. The density matrix is specified as a function of a few time-independent observables where their coefficients are time-dependent. The key mathematical step is the vectorization of the surprisal, the logarithm of the density matrix, at each time point of interest. The time span used depends on the required spectral resolution. The entire course of the system evolution is represented as a matrix where each column is the vectorized surprisal at the given time point. Using the singular value decomposition (SVD) of this matrix, we generate realistic approximations for the time-independent observables and their respective time-dependent coefficients. This allows for a simplification of the algebraic procedure for determining the dominant constraints (the time-independent observables) in the sense of the maximal entropy approach. A non-stationary coherent initial state of a Morse oscillator is used to introduce the approach. We derive the analytical exact expression for the surprisal as a function of time, and this offers a benchmark for comparison with the accurate but approximate SVD results. We discuss two examples of a Morse potential of different anharmonicities, H2 and I2 molecules. We further demonstrate the approach for a two-coupled electronic state problem, the well-studied non-radiative decay of pyrazine from its bright state. Five constraints are found to be enough to capture the ultrafast electronic population exchange and to recover the dynamics of the wave packet in both electronic states.

Time-dependent view of an isotope effect in electron-nuclear nonequilibrium dynamics with applications to N2.

Isotopic fractionation in the photodissociation of N2 could explain the considerable variation in the 14N/15N ratio in different regions of our galaxy. We previously proposed that such an isotope effect is due to coupling of photoexcited bound valence and Rydberg electronic states in the frequency range where there is strong state mixing. We here identify features of the role of the mass in the dynamics through a time-dependent quantum-mechanical simulation. The photoexcitation of N2 is by an ultrashort pulse so that the process has a sharply defined origin in time and so that we can monitor the isolated molecule dynamics in time. An ultrafast pulse is necessarily broad in frequency and spans several excited electronic states. Each excited molecule is therefore not in a given electronic state but in a superposition state. A short time after excitation, there is a fairly sharp onset of a mass-dependent large population transfer when wave packets on two different electronic states in the same molecule overlap. This coherent overlap of the wave packets on different electronic states in the region of strong coupling allows an effective transfer of population that is very mass dependent. The extent of the transfer depends on the product of the populations on the two different electronic states and on their relative phase. It is as if two molecules collide but the process occurs within one molecule, a molecule that is simultaneously in both states. An analytical toy model recovers the (strong) mass and energy dependence.

Ultrafast 25-fs relaxation in highly excited states of methyl azide mediated by strong nonadiabatic coupling.

Highly excited electronic states are challenging to explore experimentally and theoretically-due to the large density of states and the fact that small structural changes lead to large changes in electronic character with associated strong nonadiabatic dynamics. They can play a key role in astrophysical and ionospheric chemistry, as well as the detonation chemistry of high-energy density materials. Here, we implement ultrafast vacuum-UV (VUV)-driven electron-ion coincidence imaging spectroscopy to directly probe the reaction pathways of highly excited states of energetic molecules-in this case, methyl azide. Our data, combined with advanced theoretical simulations, show that photoexcitation of methyl azide by a 10-fs UV pulse at 8 eV drives fast structural changes and strong nonadiabatic coupling that leads to relaxation to other excited states on a surprisingly fast timescale of 25 fs. This ultrafast relaxation differs from dynamics occurring on lower excited states, where the timescale required for the wavepacket to reach a region of strong nonadiabatic coupling is typically much longer. Moreover, our theoretical calculations show that ultrafast relaxation of the wavepacket to a lower excited state occurs along one of the conical intersection seams before reaching the minimum energy conical intersection. These findings are important for understanding the unique strongly coupled non-Born-Oppenheimer molecular dynamics of VUV-excited energetic molecules. Although such observations have been predicted for many years, this study represents one of the few where such strongly coupled non-Born-Oppenheimer molecular dynamics of VUV-excited energetic molecules have been conclusively observed directly, making it possible to identify the ultrafast reaction pathways.

Understanding chemical interaction between phosphonate-derivative molecules and a silver surface cluster in SERS: a combined experimental and computational approach.

The interaction between phosphonate functions and a silver surface cluster is investigated using Surface-Enhanced Raman Spectroscopy (SERS). Changing the functional group (methylphosphonic acid based molecule) by studying the effect of protonation, methylation and substitution of the side chain with amine and carboxylate functions enabled us to modulate the chemical interactions between the different functions and the metal cluster. We find that the adsorption energy of the methylphosphonic acid decreases with the protonation, the methylation processes and the substitution of the side chain. In all cases, only the deprotonated phosphonate forms are SERS active. To understand how the molecules interact with the nanoparticle, the electronic structure, adsorption energies and Raman spectra were computed for molecules adsorbed on a 20 atom silver cluster representing a nanoparticle surface. The qualitative agreement between computed static Raman spectra and experimental SERS spectra makes it possible to determine stable geometries of the analyte-silver cluster complexes and to characterize the adsorption modes. The findings presented here provide a framework for designing analytical developments based on SERS for simultaneous detection of phosphonated molecules, including pesticides such as glyphosate, creating practical opportunities in key areas such as environmental and water resource in situ monitoring.

Time resolved mechanism of the isotope selectivity in the ultrafast light induced dissociation in N2.

The time evolution of a vacuum ultraviolet excited N2 molecule is followed all the way from an ultrafast excitation to dissociation by a quantum mechanical simulation. The primary aim is to discern the role of the excitation by a pulse short compared to the vibrational period, to discern the different coupling mechanisms between different electronic states, nonadiabatic, spin orbit, and to analyze the origin of any isotopic effect. We compare the picture in the time and energy domains. The initial ultrafast excitation pumps the molecule to a coherent electronic wave packet to which several singlet bound electronic states contribute. The total nonstationary wave function is given as a coherent sum of nuclear wave packets on each electronic state times the stationary electronic wave function. When the wave packets on different electronic states overlap, they are coupled in a mass-dependent manner whether one uses an adiabatic or a diabatic electronic basis. A weak spin-orbit coupling acts as a bottleneck between the bound singlet part of phase space and the triplet manifold of states in which dissociation takes place. To describe the spin-orbit perturbation that is ongoing in time, an energy-resolved eigenstate representation appears to be more intuitive. In the eigenstate basis, the singlet-to-triplet population transfer is large only between those vibronic eigenstates that are quasiresonant in energy. The states in resonance are different for different excitation energy ranges. The resonances are mass dependent, which explains the control of the isotope effect through the profile of the pulse.

Intercommunication of DNA-Based Constitutional Dynamic Networks.

Intercommunication between dynamic chemical networks plays a major role in cellular transformations. Inspired by nature, we introduce the intercommunication between two constitutional dynamic networks, CDNs, "S" and "T" composed, each, of four equilibrated supramolecular constituents AA', AB', BA', and BB', and of CC', CD', DC', and DD', respectively. Each of the constituents is conjugated to a Mg2+-ion-dependent DNAzyme unit that acts as a reporter element for the concentration of the respective constituent via the catalyzed cleavage of the fluorophore/quencher-functionalized substrate associated with the respective DNAzyme reporter. Also, constituents BB' (in CDN "S") and CC' (in CDN "T") include Mg2+-ion-dependent DNAzymes acting as activator units for generating triggering signals between the networks. Subjecting CDNs "S" and "T" to the catalytically cleavable hairpin trigger Hdd' or Haa', respectively, yields input strands that intercommunicate the CDNs by affecting the time-dependent re-equilibration of the constituents of the counter CDN without affecting the dynamic equilibrium of the constituents of the CDN that generates the triggering strands. Treatment of CDNs "S" and "T" with hairpins Hdd' and Haa' (or Hba'), respectively, stimulates autonomous positive/positive or positive/negative feedback to the programmed time-dependent up-regulation or down-regulation of the equilibrated constituents in the two CDNs.

Gas-Phase Photoluminescence and Photodissociation of Silver-Capped Hexagold Clusters.

We report on the radiative and nonradiative deactivation pathways of selected charge states of the stoichiometric hexagold phosphine-stabilized ionic clusters, [(C)(AuDppy)6Ag2·(BF4) x](4- x)+ with x = 2 and 3 (Dppy = diphenylphosphino-2-pyridine), combining gas-phase photoluminescence and photodissociation with quantum chemical computations. These clusters possess an identical isostructural core made of a hyper-coordinated carbon at their center octahedrally surrounded by six gold ions, and two silver ions at their apexes. Their luminescence and fragmentation behavior upon photoexcitation was investigated under mass and charge control in an ion trap. The experimental and computational results shed light on the electronic states involved in the optical transitions as well as on their core, ligand, or charge transfer character. Gas-phase results are discussed in relation with condensed phase measurements, as well as previous observations in solution and on metal-organic frameworks. The monocationic species ( x = 3) is found to be less stable than the dicationic one ( x = 2). In the luminescence spectrum of the monocationic species, a shoulder at short wavelength can be observed and is assigned to fragment emission. This fragment formation appears to be favored for the monocation by the existence of a low lying singlet state energetically overlapping with the triplet state manifold, which is populated quickly after photoexcitation.

Photoinduced Ultrafast Charge Transfer and Charge Migration in Small Gold Clusters Passivated by a Chromophoric Ligand.

Because the development of attopulses, charge migration induced by short optical pulses has been extensively investigated. We report a computational purely electronic dynamical study of ultrafast few femtoseconds (fs) charge transfer and charge migration in realistic passivated stoichiometric Au11 and Au20 gold nanoclusters functionalized by a bipyridine ligand. We show that a net significant amount of electronic charge (0.1 to 0.4 |e| where |e| is the electron charge) is permanently transferred from the bipyridine chromophore to the gold cluster during the short 5-6 fs UV-vis strong pulse. This electron transfer to the metallic core is induced by the optical excitation of electronic states with a partial charge transfer character involving the chromophore before the onset of nuclei motion. In addition, the photoexcitation by the strong fs pulse builds a nonequilibrium electronic density that beats between the chromophore and the metallic core around the average of the transferred value. Modular systems made of a donor chromophore that can be photoexcited in the UV-vis range coupled to an efficient acceptor that could trap the charge are of interest for applications to nanodevices. Our study provides understanding on the very early, purely electronic dynamics built by the fs optical excitation and the initial charge separation step.

Implementation of Probabilistic Algorithms by Multi-chromophoric Molecular Networks with Application to Multiple Travelling Pathways.

The implementation of probabilistic algorithms by deterministic hardware is demanding and requires hundreds of instructions to generate a pseudo-random sequence of numbers. On the contrary, the dynamics at the molecular scale is physically governed by probabilistic laws because of the stochastic nature of thermally activated and quantum processes. By simulating the exciton transfer dynamics in a multi-chromophoric system, we demonstrate the implementation of a random walk that samples the possible pathways of a traveler through a network and can be probed by time-resolved fluorescence spectroscopy. The ability of controlling the spatial arrangement of the chromophores allows us to design the "landscape" in which the traveler is moving and therefore to program the molecular device.

Implementation of Multivariable Logic Functions in Parallel by Electrically Addressing a Molecule of Three Dopants in Silicon.

To realize low-power, compact logic circuits, one can explore parallel operation on single nanoscale devices. An added incentive is to use multivalued (as distinct from Boolean) logic. Here, we theoretically demonstrate that the computation of all the possible outputs of a multivariate, multivalued logic function can be implemented in parallel by electrical addressing of a molecule made up of three interacting dopant atoms embedded in Si. The electronic states of the dopant molecule are addressed by pulsing a gate voltage. By simulating the time evolution of the non stationary electronic density built by the gate voltage, we show that one can implement a molecular decision tree that provides in parallel all the outputs for all the inputs of the multivariate, multivalued logic function. The outputs are encoded in the populations and in the bond orders of the dopant molecule, which can be measured using an STM tip. We show that the implementation of the molecular logic tree is equivalent to a spectral function decomposition. The function that is evaluated can be field-programmed by changing the time profile of the pulsed gate voltage.

Nuclear Motion Driven Ultrafast Photodissociative Charge Transfer of the PENNA Cation: An Experimental and Computational Study.

Ultrafast nuclear driven charge transfer prior to dissociation is an important process in modular systems as was demonstrated experimentally in the bifunctional molecule 2-phenylethyl-N,N-dimethylamine (PENNA) in work by Lehr et al. ( J. Phys. Chem. A 2005 , 109 , 8074 ). The ultrafast dynamics of PENNA photoexcited to the three lowest electronic states of the cation (D0, D1, and D2) was studied using quantum chemistry and surface hoping. We show that a conical intersection, localized in the Franck-Condon region, between the D0 and the D1 states, leads to an ultrafast charge transfer, computed here to be on a time scale of 65 fs, between the phenyl and the amine charged subunits. On the D0 ground state, the dissociation proceeds on the 60 ps time scale through a 19 kcal/mol late barrier. The computed kinetic energy release is in good agreement with a new experimental measurement of PENNA ionization by an 800 nm 30 fs intense laser pulse.

A new electron-ion coincidence 3D momentum-imaging method and its application in probing strong field dynamics of 2-phenylethyl-N, N-dimethylamine.

We report the development of a new three-dimensional (3D) momentum-imaging setup based on conventional velocity map imaging to achieve the coincidence measurement of photoelectrons and photo-ions. This setup uses only one imaging detector (microchannel plates (MCP)/phosphor screen) but the voltages on electrodes are pulsed to push both electrons and ions toward the same detector. The ion-electron coincidence is achieved using two cameras to capture images of ions and electrons separately. The time-of-flight of ions and electrons are read out from MCP using a digitizer. We demonstrate this new system by studying the dissociative single and double ionization of PENNA (2-phenylethyl-N,N-dimethylamine). We further show that the camera-based 3D imaging system can operate at 10 kHz repetition rate.