No Arabic abstract
We investigate the molecular geometries of the ground state and the minimal energy conical intersections (MECIs) between the ground and first excited states of the models for the retinal protonated Schiff base in the gas phase using the extended multistate complete active space second-order perturbation theory (XMS-CASPT2). The biggest model in this work is the rhodopsin chromophore truncated between the {epsilon} and {delta} carbon atoms, which consists of 54 atoms and 12-orbital {pi} conjugation. The results are compared with those obtained by the state-averaged complete active space self-consistent field (SA-CASSCF). The XMS-CASPT2 results suggest that the minimum energy conical intersection associated with the so-called 13-14 isomerization is thermally inaccessible, which is in contrast to the SA-CASSCF results. The differences between the geometries of the conical intersections computed by SA-CASSCF and XMS-CASPT2 are ascribed to the fact that the charge transfer states are more stabilized by dynamical electron correlation than the diradicaloid states. The impact of the various choices of active spaces, basis sets, and state averaging schemes is also examined.
The fitting of physical models is often done only using a single target observable. However, when multiple targets are considered, the fitting procedure becomes cumbersome, there being no easy way to quantify the robustness of the model for all different observables. Here, we illustrate that one can jointly search for the best model for each desired observable through multi-objective optimization. To do so we construct the Pareto front to study if there exists a set of parameters of the model that can jointly describe multiple, or all, observables. To alleviate the computational cost, the predicted error for each targeted objective is approximated with a Gaussian process model, as it is commonly done in the Bayesian optimization framework. We applied this methodology to improve three different models used in the simulation of stationary state $cis-trans$ photoisomerization of retinal in rhodopsin. Optimization was done with respect to different experimental measurements, including emission spectra, peak absorption frequencies for the $cis$ and $trans$ conformers, and the energy storage.
A new one-dimensional model is proposed for the low-energy vibrational quantum dynamics of CH5+ based on the motion of an effective particle confined to a 60-vertex graph ${Gamma}_{60}$ with a single edge length parameter. Within this model, the quantum states of CH5+ are obtained in analytic form and are related to combinatorial properties of ${Gamma}_{60}$. The bipartite structure of ${Gamma}_{60}$ gives a simple explanation for curious symmetries observed in numerically exact variational calculations on CH5+.
The SEHRS spectrum of 4,4 Bipyridine is analyzed on the base of the Dipole Quadrupole theory. It is demonstrated that there appear strong lines caused by vibrations transforming after a unit irreducible representation of the D2 symmetry group. which is most probably describes the symmetry properties of the molecule. These lines are nearly forbidden for the molecule, adsorbed on rough metal surface. Appearance of these lines is associated with a strong quadrupole light molecule interaction, which exists in this system. In addition, there are lines, caused by contributions from both the vibrations transforming after the unit irreducible representation A and the representation B1, which describes transformational properties of the dz component of the dipole moment, which is perpendicular to the surface. This result is associated with the specific geometry of the molecule, when the indicated vibrations can be nearly degenerated and cannot be resolved by the SEHRS spectra analysis. Analysis of the SEHRS spectra for the possible geometry of the molecule with the D2h symmetry group leads to similar results. This issue is in a full coincidence with the results of the SEHRS Dipole Quadrupole theory.
We consider the emergent behavior of viral spread when agents in a large population interact with each other over a contact network. When the number of agents is large and the contact network is a complete graph, it is well known that the population behavior -- that is, the fraction of susceptible, infected and recovered agents -- converges to the solution of an ordinary differential equation (ODE) known as the classical SIR model as the population size approaches infinity. In contrast, we study interactions over contact networks with generic topologies and derive conditions under which the population behavior concentrates around either the classic SIR model or other deterministic models. Specifically, we show that when most vertex degrees in the contact network are sufficiently large, the population behavior concentrates around an ODE known as the network SIR model. We then study the short and intermediate-term evolution of the network SIR model and show that if the contact network has an expander-type property or the initial set of infections is well-mixed in the population, the network SIR model reduces to the classical SIR model. To complement these results, we illustrate through simulations that the two models can yield drastically different predictions, hence use of the classical SIR model can be misleading in certain cases.
In nanopore sequencing, electrical signal is measured as DNA molecules pass through the sequencing pores. Translating these signals into DNA bases (base calling) is a highly non-trivial task, and its quality has a large impact on the sequencing accuracy. The most successful nanopore base callers to date use convolutional neural networks (CNN) to accomplish the task. Convolutional layers in CNNs are typically composed of filters with constant window size, performing best in analysis of signals with uniform speed. However, the speed of nanopore sequencing varies greatly both within reads and between sequencing runs. Here, we present dynamic pooling, a novel neural network component, which addresses this problem by adaptively adjusting the pooling ratio. To demonstrate the usefulness of dynamic pooling, we developed two base callers: Heron and Osprey. Heron improves the accuracy beyond the experimental high-accuracy base caller Bonito developed by Oxford Nanopore. Osprey is a fast base caller that can compete in accuracy with Guppy high-accuracy mode, but does not require GPU acceleration and achieves a near real-time speed on common desktop CPUs. Availability: https://github.com/fmfi-compbio/osprey, https://github.com/fmfi-compbio/heron Keywords: nanopore sequencing, base calling, convolutional neural networks, pooling