Search

In this paper we discuss the performance of the non-iterative state-specific multireference coupled cluster (SS-MRCC) methods accounting for the effect of triply excited cluster amplitudes. The corrections to the Brillouin-Wigner and Mukherjee's MRCC models based on the manifold of singly and doubly excited cluster amplitudes (BW-MRCCSD and Mk-MRCCSD, respectively) are tested and compared with exact full configuration interaction results for small systems (H2O, N2, and Be3). For the larger systems (naphthyne isomers) the BW-MRCC and Mk-MRCC methods with iterative singles, doubles, and non-iterative triples (BW-MRCCSD(T) and Mk-MRCCSD(T)) are compared against the results obtained with single reference coupled cluster methods. We also report on the parallel performance of the non-iterative implementations based on the use of processor groups. [ABSTRACT FROM AUTHOR]

Coarse-graining (CG) techniques have recently attracted great interest for providing descriptions at a mesoscopic level of resolution that preserve fluid thermodynamic and transport behaviors with a reduced number of degrees of freedom and hence less computational effort. One fundamental question arises: how well and to what extent can a 'bottom-up' developed mesoscale model recover the physical properties of a molecular scale system? To answer this question, we explore systematically the properties of a CG model that is developed to represent an intermediate mesoscale model between the atomistic and continuum scales. This CG model aims to reduce the computational cost relative to a full atomistic simulation, and we assess to what extent it is possible to preserve both the thermodynamic and transport properties of an underlying reference all-atom Lennard-Jones (LJ) system. In this paper, only the thermodynamic properties are considered in detail. The transport properties will be examined in subsequent work. To coarse-grain, we first use the iterative Boltzmann inversion (IBI) to determine a CG potential for a (1-[lowercase_phi_synonym])N mesoscale particle system, where [lowercase_phi_synonym] is the degree of coarse-graining, so as to reproduce the radial distribution function (RDF) of an N atomic particle system. Even though the uniqueness theorem guarantees a one to one relationship between the RDF and an effective pairwise potential, we find that RDFs are insensitive to the long-range part of the IBI-determined potentials, which provides some significant flexibility in further matching other properties. We then propose a reformulation of IBI as a robust minimization procedure that enables simultaneous matching of the RDF and the fluid pressure. We find that this new method mainly changes the attractive tail region of the CG potentials, and it improves the isothermal compressibility relative to pure IBI. We also find that there are optimal interaction cutoff lengths for the CG system, as a function of [lowercase_phi_synonym], that are required to attain an adequate potential while maintaining computational speedup. To demonstrate the universality of the method, we test a range of state points for the LJ liquid as well as several LJ chain fluids. [ABSTRACT FROM AUTHOR]

As a consequence of the rugged landscape of RNA molecules their folding is described by the kinetic partitioning mechanism according to which only a small fraction (φF) reaches the folded state while the remaining fraction of molecules is kinetically trapped in misfolded intermediates. The transition from the misfolded states to the native state can far exceed biologically relevant time. Thus, RNA folding in vivo is often aided by protein cofactors, called RNA chaperones, that can rescue RNAs from a multitude of misfolded structures. We consider two models, based on chemical kinetics and chemical master equation, for describing assisted folding. In the passive model, applicable for class I substrates, transient interactions of misfolded structures with RNA chaperones alone are sufficient to destabilize the misfolded structures, thus entropically lowering the barrier to folding. For this mechanism to be efficient the intermediate ribonucleoprotein complex between collapsed RNA and protein cofactor should have optimal stability. We also introduce an active model (suitable for stringent substrates with small φF), which accounts for the recent experimental findings on the action of CYT-19 on the group I intron ribozyme, showing that RNA chaperones do not discriminate between the misfolded and the native states. In the active model, the RNA chaperone system utilizes chemical energy of adenosine triphosphate hydrolysis to repeatedly bind and release misfolded and folded RNAs, resulting in substantial increase of yield of the native state. The theory outlined here shows, in accord with experiments, that in the steady state the native state does not form with unit probability. [ABSTRACT FROM AUTHOR]