Dissipative dynamics of a harmonic oscillator: A nonperturbative approach

Starting from a microscopic theory, we derive a master equation for a harmonic oscillator coupled to a bath of noninteracting oscillators. We follow a nonperturbative approach, proposed earlier by us for the free Brownian particle. The diffusion constants are calculated analytically and the positivity of the master equation is shown to hold above a critical temperature. We compare the long time behavior of the average kinetic and potential energies with known thermodynamic results. In the limit of vanishing oscillator frequency of the system, we recover the results of the free Brownian particle.

Photophysics of conjugated polymers: interplay between Forster energy migration and defect concentration in shaping a photochemical funnel in PPV

Recent single molecule experiments have suggested the existence of a photochemical funnel in the photophysics of conjugated polymers, like poly[2-methoxy-5-(2'-ethylhexyl)oxy-1,4-phenylenevinylene] (MEH-PPV). The funnel is believed to be a consequence of the presence of conformational or chemical defects along the polymer chain and efficient non-radiative energy transfer among different chromophore segments. Here we address the effect of the excitation energy dynamics on the photophysics of PPV. The PPV chain is modeled as a polymer with the length distribution of chromophores given either by a Gaussian or by a Poisson distribution. We observe that the Poisson distribution of the segment lengths explains the photophysics of PPV better than the Gaussian distribution. A recently proposed version of an extended particle-in-a-box' model is used to calculate the exciton energies and the transition dipole moments of the chromophores, and a master equation to describe the excitation energy transfer among different chromophores. The rate of energy transfer is assumed to be given here, as a first approximation, by the well-known Forster expression. The observed excitation population dynamics confirms the photochemical funneling of excitation energy from shorter to longer chromophores of the polymer chain. The time scale of spectral shift and energy transfer for our model polymer, with realistic values of optical parameters, is in the range of 200-300 ps. We find that the excitation energy may not always migrate towards the longest chromophore segments in the polymer chain as the efficiency of energy transfer between chromophores depends on the separation distance between the two and their relative orientation.

Monte carlo simulation for molecular gas dynamics

The dynamics of low-density flows is governed by the Boltzmann equation of the kinetic theory of gases. This is a nonlinear integro-differential equation and, in general, numerical methods must be used to obtain its solution. The present paper, after a brief review of Direct Simulation Monte Carlo (DSMC) methods due to Bird, and Belotserkovskii and Yanitskii, studies the details of theDSMC method of Deshpande for mono as well as multicomponent gases. The present method is a statistical particle-in-cell method and is based upon the Kac-Prigogine master equation which reduces to the Boltzmann equation under the hypothesis of molecular chaos. The proposed Markoff model simulating the collisions uses a Poisson distribution for the number of collisions allowed in cells into which the physical space is divided. The model is then extended to a binary mixture of gases and it is shown that it is necessary to perform the collisions in a certain sequence to obtain unbiased simulation.

Characterization of the dynamics of glass-forming liquids from the properties of the potential energy landscape

We develop a framework for understanding the difference between strong and fragile behavior in the dynamics of glass-forming liquids from the properties of the potential energy landscape. Our approach is based on a master equation description of the activated jump dynamics among the local minima of the potential energy (the so-called inherent structures) that characterize the potential energy landscape of the system. We study the dynamics of a small atomic cluster using this description as well as molecular dynamics simulations and demonstrate the usefulness of our approach for this system. Many of the remarkable features of the complex dynamics of glassy systems emerge from the activated dynamics in the potential energy landscape of the atomic cluster. The dynamics of the system exhibits typical characteristics of a strong supercooled liquid when the system is allowed to explore the full configuration space. This behavior arises because the dynamics is dominated by a few lowest-lying minima of the potential energy and the potential energy barriers between these minima. When the system is constrained to explore only a limited region of the potential energy landscape that excludes the basins of attraction of a few lowest-lying minima, the dynamics is found to exhibit the characteristics of a fragile liquid.

Adiabatic eigenfunction-based approach for coherent excitation transfer: An almost analytical treatment of the Fenna-Matthews-Olson complex

We suggest a method of studying coherence in finite-level systems coupled to the environment and use it for the Hamiltonian that has been used to describe the light-harvesting pigment-protein complex. The method works with the adiabatic states and transforms the Hamiltonian to a form in which the terms responsible for decoherence and population relaxation are separated out. Decoherence is then accounted for nonperturbatively and population relaxation using a Markovian master equation. Almost analytical results can be obtained for the seven-level system, and the calculations are very simple for systems with more levels. We apply the treatment to the seven-level system, and the results are in excellent agreement with the exact numerical results of Nalbach et al. Nalbach, Braun, and Thorwart, Phys. Rev. E 84, 041926 (2011)]. Our approach is able to account for decoherence and population relaxation separately. It is found that decoherence causes only damping of oscillations and does not lead to transfer to the reaction center. Population relaxation is necessary for efficient transfer to the reaction center, in agreement with earlier findings. Our results show that the transformation to the adiabatic basis followed by a Redfield type of approach leads to results in good agreement with exact simulation.

Adiabatic Eigenfunction Based Approach to Coherent Transfer: Application to the Fenna-Matthews-Olson (FMO) Complex and the Role of Correlations in the Efficiency of Energy Transfer

We have recently suggested a method (Pallavi Bhattacharyya and K. L. Sebastian, Physical Review E 2013, 87, 062712) for the analysis of coherence in finite-level systems that are coupled to the surroundings and used it to study the process of energy transfer in the Fenna-Matthews-Olson (FMO) complex. The method makes use of adiabatic eigenstates of the Hamiltonian, with a subsequent transformation of the Hamiltonian into a form where the terms responsible for decoherence and population relaxation could be separated out at the lowest order. Thus one can account for decoherence nonperturbatively, and a Markovian type of master equation could be used for evaluating the population relaxation. In this paper, we apply this method to a two-level system as well as to a seven-level system. Comparisons with exact numerical results show that the method works quite well and is in good agreement with numerical calculations. The technique can be applied with ease to systems with larger numbers of levels as well. We also investigate how the presence of correlations among the bath degrees of freedom of the different bacteriochlorophyll a molecules of the FMO Complex affect the rate of energy transfer. Surprisingly, in the cases that we studied, our calculations suggest that the presence of anticorrelations, in contrast to correlations, make the excitation transfer more facile.

Globally coupled stochastic two-state oscillators: Fluctuations due to finite numbers

Infinite arrays of coupled two-state stochastic oscillators exhibit well-defined steady states. We study the fluctuations that occur when the number N of oscillators in the array is finite. We choose a particular form of global coupling that in the infinite array leads to a pitchfork bifurcation from a monostable to a bistable steady state, the latter with two equally probable stationary states. The control parameter for this bifurcation is the coupling strength. In finite arrays these states become metastable: The fluctuations lead to distributions around the most probable states, with one maximum in the monostable regime and two maxima in the bistable regime. In the latter regime, the fluctuations lead to transitions between the two peak regions of the distribution. Also, we find that the fluctuations break the symmetry in the bimodal regime, that is, one metastable state becomes more probable than the other, increasingly so with increasing array size. To arrive at these results, we start from microscopic dynamical evolution equations from which we derive a Langevin equation that exhibits an interesting multiplicative noise structure. We also present a master equation description of the dynamics. Both of these equations lead to the same Fokker-Planck equation, the master equation via a 1/N expansion and the Langevin equation via standard methods of Ito calculus for multiplicative noise. From the Fokker-Planck equation we obtain an effective potential that reflects the transition from the monomodal to the bimodal distribution as a function of a control parameter. We present a variety of numerical and analytic results that illustrate the strong effects of the fluctuations. We also show that the limits N -> infinity and t -> infinity(t is the time) do not commute. In fact, the two orders of implementation lead to drastically different results.

Frequency-domain stimulated and spontaneous light emission signals at molecular junctions

Using a diagrammatic superoperator formalism we calculate optical signals at molecular junctions where a single molecule is coupled to two metal leads which are held at different chemical potentials. The molecule starts in a nonequilibrium steady state whereby it continuously exchanges electrons with the leads with a constant electron flux. Expressions for frequency domain optical signals measured in response to continuous laser fields are derived by expanding the molecular correlation functions in terms of its many-body states. The nonunitary evolution of molecular states is described by the quantum master equation. (C) 2014 AIP Publishing LLC.

Statistics of an adiabatic charge pump

We investigate the effect of time-dependent cyclic-adiabatic driving on the charge transport in a quantum junction. We propose a nonequilibrium Green's function formalism to study the statistics of the charge pumped (at zero bias) through the junction. The formulation is used to demonstrate charge pumping in a single electronic level coupled to two (electronic) reservoirs with time-dependent couplings. An analytical expression for the average pumped current for a general cyclic driving is derived. It is found that for zero bias, for a certain class of driving, the Berry phase contributes only to the odd cumulants. In contrast, a quantum master equation formulation does not show a Berry-phase effect at all.

A comment on the consistency of truncated nonlinear integral equation based theories of freezing

We report the results of two studies of aspects of the consistency of truncated nonlinear integral equation based theories of freezing: (i) We show that the self-consistent solutions to these nonlinear equations are unfortunately sensitive to the level of truncation. For the hard sphere system, if the Wertheim–Thiele representation of the pair direct correlation function is used, the inclusion of part but not all of the triplet direct correlation function contribution, as has been common, worsens the predictions considerably. We also show that the convergence of the solutions found, with respect to number of reciprocal lattice vectors kept in the Fourier expansion of the crystal singlet density, is slow. These conclusions imply great sensitivity to the quality of the pair direct correlation function employed in the theory. (ii) We show the direct correlation function based and the pair correlation function based theories of freezing can be cast into a form which requires solution of isomorphous nonlinear integral equations. However, in the pair correlation function theory the usual neglect of the influence of inhomogeneity of the density distribution on the pair correlation function is shown to be inconsistent to the lowest order in the change of density on freezing, and to lead to erroneous predictions. The Journal of Chemical Physics is copyrighted by The American Institute of Physics.

A comment on the consistency of truncated nonlinear integral equation based theories of freezing

We report the results of two studies of aspects of the consistency of truncated nonlinear integral equation based theories of freezing: (i) We show that the self-consistent solutions to these nonlinear equations are unfortunately sensitive to the level of truncation. For the hard sphere system, if the Wertheim–Thiele representation of the pair direct correlation function is used, the inclusion of part but not all of the triplet direct correlation function contribution, as has been common, worsens the predictions considerably. We also show that the convergence of the solutions found, with respect to number of reciprocal lattice vectors kept in the Fourier expansion of the crystal singlet density, is slow. These conclusions imply great sensitivity to the quality of the pair direct correlation function employed in the theory. (ii) We show the direct correlation function based and the pair correlation function based theories of freezing can be cast into a form which requires solution of isomorphous nonlinear integral equations. However, in the pair correlation function theory the usual neglect of the influence of inhomogeneity of the density distribution on the pair correlation function is shown to be inconsistent to the lowest order in the change of density on freezing, and to lead to erroneous predictions. The Journal of Chemical Physics is copyrighted by The American Institute of Physics.

Analysis of Paschen Curves for air, N2 and SF6 Using the Townsend Breakdown Equation

It has been shown that it is possible to extend the validity of the Townsend breakdown criterion for evaluating the breakdown voltages in the complete pd range in which Paschen curves are available. Evaluation of the breakdown voltages for air (pd=0.0133 to 1400 kPa · cm), N2(pd=0.0313 to 1400 kPa · cm) and SF6 (pd=0.3000 to 1200 kPa · cm) has been done and in most cases the computed values are accurate to ±3% of the measured values. The computations show that it is also possible to estimate the secondary ionization coefficient ¿ in the pd ranges mentioned above.

Multinomial solutions of the beach equation

Exact multinomial solutions of the beach equation for shallow water waves on a uniformly sloping beach are found and related to solution of the same equation found earlier by other investigators, using integral transform techniques. The use of these solutions for a general initialvalue problem for the equation under investigation is briefly discussed.

Convenient construction of tetrahedral finite elements for the quasi-harmonic equation

It is shown that in the finite-element formulation of the general quasi-harmonic equation using tetrahedral elements, for every member of the element family there exists just one numerical universal matrix indpendent of the size, shape and material properties of the element. Thus the element matrix is conveniently constructed by manipulating this single matrix along with a set of reverse sequence codes at the same time accounting for the size, shape and material properties in a simple manner.

Multinomial solutions of the beach equation

Exact multinomial solutions of the beach equation for shallow water waves on a uniformly sloping beach are found and related to solution of the same equation found earlier by other investigators, using integral transform techniques. The use of these solutions for a general initialvalue problem for the equation under investigation is briefly discussed.