Simple finite field method for calculation of static and dynamic vibrational hyperpolarizabilities: curvature contributions

In the static field limit, the vibrational hyperpolarizability consists of two contributions due to: (1) the shift in the equilibrium geometry (known as nuclear relaxation), and (2) the change in the shape of the potential energy surface (known as curvature). Simple finite field methods have previously been developed for evaluating these static field contributions and also for determining the effect of nuclear relaxation on dynamic vibrational hyperpolarizabilities in the infinite frequency approximation. In this paper the finite field approach is extended to include, within the infinite frequency approximation, the effect of curvature on the major dynamic nonlinear optical processes

Automatic generation of active coordinates for quantum dynamics calculations: application to the dynamics of benzene photochemistry

A new practical method to generate a subspace of active coordinates for quantum dynamics calculations is presented. These reduced coordinates are obtained as the normal modes of an analytical quadratic representation of the energy difference between excited and ground states within the complete active space self-consistent field method. At the Franck-Condon point, the largest negative eigenvalues of this Hessian correspond to the photoactive modes: those that reduce the energy difference and lead to the conical intersection; eigenvalues close to 0 correspond to bath modes, while modes with large positive eigenvalues are photoinactive vibrations, which increase the energy difference. The efficacy of quantum dynamics run in the subspace of the photoactive modes is illustrated with the photochemistry of benzene, where theoretical simulations are designed to assist optimal control experiments

Determination of vibrational polarizabilities and hyperpolarizabilities using field-induced coordinates

An analytical set of field-induced coordinates is defined and is used to show that the vibrational degrees of freedom required to completely describe nuclear relaxation polarizabilities and hyperpolarizabilities is reduced from 3N-6 to a relatively small number. As this number does not depend upon the size of the molecule, the process provides computational advantages. A method is provided to separate anharmonic contributions from harmonic contributions as well as effective mechanical from electrical anharmonicity. The procedures are illustrated by Hartree-Fock calculations, indicating that anharmonicity can be very important

Exploring unfrequent events with accelerated molecular dynamics simulations: from photochemistry to drug design

La meva incorporació al grup de recerca del Prof. McCammon (University of California San Diego) en qualitat d’investigador post doctoral amb una beca Beatriu de Pinós, va tenir lloc el passat 1 de desembre de 2010; on vaig dur a terme les meves tasques de recerca fins al darrer 1 d’abril de 2012. El Prof. McCammon és un referent mundial en l’aplicació de simulacions de dinàmica molecular (MD) en sistemes biològics d’interès humà. La contribució més important del Prof. McCammon en la simulació de sistemes biològics és el desenvolupament del mètode de dinàmiques moleculars accelerades (AMD). Les simulacions MD convencionals, les quals estan limitades a l’escala de temps del nanosegon (~10-9s), no son adients per l’estudi de sistemes biològics rellevants a escales de temps mes llargues (μs, ms...). AMD permet explorar fenòmens moleculars poc freqüents però que son clau per l’enteniment de molts sistemes biològics; fenòmens que no podrien ser observats d’un altre manera. Durant la meva estada a la “University of California San Diego”, vaig treballar en diferent aplicacions de les simulacions AMD, incloent fotoquímica i disseny de fàrmacs per ordinador. Concretament, primer vaig desenvolupar amb èxit una combinació dels mètodes AMD i simulacions Car-Parrinello per millorar l’exploració de camins de desactivació (interseccions còniques) en reaccions químiques fotoactivades. En segon lloc, vaig aplicar tècniques estadístiques (Replica Exchange) amb AMD en la descripció d’interaccions proteïna-lligand. Finalment, vaig dur a terme un estudi de disseny de fàrmacs per ordinador en la proteïna-G Rho (involucrada en el desenvolupament de càncer humà) combinant anàlisis estructurals i simulacions AMD. Els projectes en els quals he participat han estat publicats (o estan encara en procés de revisió) en diferents revistes científiques, i han estat presentats en diferents congressos internacionals. La memòria inclosa a continuació conté més detalls de cada projecte esmentat.

Finite horizon bargaining and the consistent field

This paper explores the relationships between noncooperative bargaining games and the consistent value for non-transferable utility (NTU) cooperative games. A dynamic approach to the consistent value for NTU games is introduced: the consistent vector field. The main contribution of the paper is to show that the consistent field is intimately related to the concept of subgame perfection for finite horizon noncooperative bargaining games, as the horizon goes to infinity and the cost of delay goes to zero. The solutions of the dynamic system associated to the consistent field characterize the subgame perfect equilibrium payoffs of the noncooperative bargaining games. We show that for transferable utility, hyperplane and pure bargaining games, the dynamics of the consistent fields converge globally to the unique consistent value. However, in the general NTU case, the dynamics of the consistent field can be complex. An example is constructed where the consistent field has cyclic solutions; moreover, the finite horizon subgame perfect equilibria do not approach the consistent value.

Thermal activation and ATP dependence of the cytoskeleton remodeling dynamics

The cytoskeleton (CSK) is a nonequilibrium polymer network that uses hydrolyzable sources of free energy such as adenosine triphosphate (ATP) to remodel its internal structure. As in inert nonequilibrium soft materials, CSK remodeling has been associated with structural rearrangements driven by energy-activated processes. We carry out particle tracking and traction microscopy measurements of alveolar epithelial cells at various temperatures and ATP concentrations. We provide the first experimental evidence that the remodeling dynamics of the CSK is driven by structural rearrangements over free-energy barriers induced by thermally activated forces mediated by ATP. The measured activation energy of these forces is ~40kBTr (kB being the Boltzmann constant and Tr being the room temperature). Our experiments provide clues to understand the analogy between the dynamics of the living CSK and that of inert nonequilibrium soft materials.

Stability of a nonequilibrium steady-state interface

We study the interfacial modes of a driven diffusive model under suitable nonequilibrium conditions leading to possible instability. The external field parallel to the interface, which sets up a steady-state parallel flux, enhances the growth or decay rates of the interfacial modes. More dramatically, asymmetry in the model can introduce an oscillatory component into the interfacial dispersion relation. In certain circumstances, the applied field behaves as a singular perturbation.

Comparison between molecular dynamics abd Monte Carlo simulations of an ordering process in a binary alloy

Ordering in a binary alloy is studied by means of a molecular-dynamics (MD) algorithm which allows to reach the domain growth regime. Results are compared with Monte Carlo simulations using a realistic vacancy-atom (MC-VA) mechanism. At low temperatures fast growth with a dynamical exponent x>1/2 is found for MD and MC-VA. The study of a nonequilibrium ordering process with the two methods shows the importance of the nonhomogeneity of the excitations in the system for determining its macroscopic kinetics.

Random-field Potts model with dipolarlike interactions: Hysteresis avalanches and microstructure

A model for the study of hysteresis and avalanches in a first-order phase transition from a single variant phase to a multivariant phase is presented. The model is based on a modification of the random-field Potts model with metastable dynamics by adding a dipolar interaction term truncated at nearest neighbors. We focus our study on hysteresis loop properties, on the three-dimensional microstructure formation, and on avalanche statistics.

Gauge-invariant actions from constraint Hamiltonian dynamics

Dirac's constraint Hamiltonian formalism is used to construct a gauge-invariant action for the massive spin-one and -two fields.

Effective field theories for heavy quarkonium

This article reviews recent theoretical developments in heavy-quarkonium physics from the point of view of effective-field theories of QCD. We discuss nonrelativistic QCD and concentrate on potential nonrelativistic QCD. The main goal will be to derive Schrödinger equations based on QCD that govern heavy-quarkonium physics in the weak- and strong-coupling regimes. Finally, the review discusses a selected set of applications, which include spectroscopy, inclusive decays, and electromagnetic threshold production.

Dynamics of domain walls in pattern formation under traveling-wave forcing

We study dynamics of domain walls in pattern forming systems that are externally forced by a moving space-periodic modulation close to 2:1 spatial resonance. The motion of the forcing induces nongradient dynamics, while the wave number mismatch breaks explicitly the chiral symmetry of the domain walls. The combination of both effects yields an imperfect nonequilibrium Ising-Bloch bifurcation, where all kinks (including the Ising-like one) drift. Kink velocities and interactions are studied within the generic amplitude equation. For nonzero mismatch, a transition to traveling bound kink-antikink pairs and chaotic wave trains occurs.

Defect dynamics in viscous fingering

We study the dynamics of Staffman-Taylor fingering in terms of topological defects of the flow field. The defects are created and/or annihilated at the interface. The route towards the single-finger steady state is characterized by a detailed mechanism for defect annihilation. For small viscosity contrast this mechanism is impeded, and creation of new defects leads the system away from a single-finger solution. Strong evidence for a drastic reduction of the basin of attraction of the Saffman-Taylor finger is presented.

Noise-induced phase separation: Mean-field results

We present a study of a phase-separation process induced by the presence of spatially correlated multiplicative noise. We develop a mean-field approach suitable for conserved-order-parameter systems and use it to obtain the phase diagram of the model. Mean-field results are compared with numerical simulations of the complete model in two dimensions. Additionally, a comparison between the noise-driven dynamics of conserved and nonconserved systems is made at the level of the mean-field approximation.

Phase-field model of Hele-Shaw flows in the high-viscosity contrast regime

A one-sided phase-field model is proposed to study the dynamics of unstable interfaces of Hele-Shaw flows in the high viscosity contrast regime. The corresponding macroscopic equations are obtained by means of an asymptotic expansion from the phase-field model. Numerical integrations of the phase-field model in a rectangular Hele-Shaw cell reproduce finger competition with the final evolution to a steady-state finger.