Hamiltonian geophysical fluid dynamics

Hamiltonian dynamics describes the evolution of conservative physical systems. Originally developed as a generalization of Newtonian mechanics, describing gravitationally driven motion from the simple pendulum to celestial mechanics, it also applies to such diverse areas of physics as quantum mechanics, quantum field theory, statistical mechanics, electromagnetism, and optics – in short, to any physical system for which dissipation is negligible. Dynamical meteorology consists of the fundamental laws of physics, including Newton’s second law. For many purposes, diabatic and viscous processes can be neglected and the equations are then conservative. (For example, in idealized modeling studies, dissipation is often only present for numerical reasons and is kept as small as possible.) In such cases dynamical meteorology obeys Hamiltonian dynamics. Even when nonconservative processes are not negligible, it often turns out that separate analysis of the conservative dynamics, which fully describes the nonlinear interactions, is essential for an understanding of the complete system, and the Hamiltonian description can play a useful role in this respect. Energy budgets and momentum transfer by waves are but two examples.

Modeling the effects of anesthesia on the electroencephalogram

Changes to the electroencephalogram (EEG) observed during general anesthesia are modeled with a physiological mean field theory of electrocortical activity. To this end a parametrization of the postsynaptic impulse response is introduced which takes into account pharmacological effects of anesthetic agents on neuronal ligand-gated ionic channels. Parameter sets for this improved theory are then identified which respect known anatomical constraints and predict mean firing rates and power spectra typically encountered in human subjects. Through parallelized simulations of the eight nonlinear, two-dimensional partial differential equations on a grid representing an entire human cortex, it is demonstrated that linear approximations are sufficient for the prediction of a range of quantitative EEG variables. More than 70 000 plausible parameter sets are finally selected and subjected to a simulated induction with the stereotypical inhaled general anesthetic isoflurane. Thereby 86 parameter sets are identified that exhibit a strong “biphasic” rise in total power, a feature often observed in experiments. A sensitivity study suggests that this “biphasic” behavior is distinguishable even at low agent concentrations. Finally, our results are briefly compared with previous work by other groups and an outlook on future fits to experimental data is provided.

A mathematical model of crystallization in an emulsion

A mathematical model incorporating many of the important processes at work in the crystallization of emulsions is presented. The model describes nucleation within the discontinuous domain of an emulsion, precipitation in the continuous domain, transport of monomers between the two domains, and formation and subsequent growth of crystals in both domains. The model is formulated as an autonomous system of nonlinear, coupled ordinary differential equations. The description of nucleation and precipitation is based upon the Becker–Döring equations of classical nucleation theory. A particular feature of the model is that the number of particles of all species present is explicitly conserved; this differs from work that employs Arrhenius descriptions of nucleation rate. Since the model includes many physical effects, it is analyzed in stages so that the role of each process may be understood. When precipitation occurs in the continuous domain, the concentration of monomers falls below the equilibrium concentration at the surface of the drops of the discontinuous domain. This leads to a transport of monomers from the drops into the continuous domain that are then incorporated into crystals and nuclei. Since the formation of crystals is irreversible and their subsequent growth inevitable, crystals forming in the continuous domain effectively act as a sink for monomers “sucking” monomers from the drops. In this case, numerical calculations are presented which are consistent with experimental observations. In the case in which critical crystal formation does not occur, the stationary solution is found and a linear stability analysis is performed. Bifurcation diagrams describing the loci of stationary solutions, which may be multiple, are numerically calculated.

Step edges in thin films of lamellar-forming diblock copolymer

Self-consistent field theory (SCFT) is used to study the step edges that occur in thin films of lamellar-forming diblock copolymer, when the surfaces each have an affinity for one of the polymer components. We examine film morphologies consisting of a stack of ν continuous monolayers and one semi-infinite bilayer, the edge of which creates the step. The line tension of each step morphology is evaluated and phase diagrams are constructed showing the conditions under which the various morphologies are stable. The predicted behavior is then compared to experiment. Interestingly, our atomic force microscopy (AFM) images of terraced films reveal a distinct change in the character of the steps with increasing ν, which is qualitatively consistent with our SCFT phase diagrams. Direct quantitative comparisons are not possible because the SCFT is not yet able to probe the large polymer/air surface tensions characteristic of experiment.

Morphology induced spinodal decomposition at the surface of symmetric diblock copolymer films

Atomic force microscopy is used to study the ordering dynamics of symmetric diblock copolymer films. The films order to form a lamellar structure which results in a frustration when the film thickness is incommensurate with the lamellae. By probing the morphology of incommensurate films in the early ordering stages, we discover an intermediate phase of lamellae arranged perpendicular to the film surface. This morphology is accompanied by a continuous growth in amplitude of the film surface topography with a characteristic wavelength, indicative of a spinodal process. Using selfconsistent field theory, we show that the observation of perpendicular lamellae suggests an intermediate state with parallel lamellae at the substrate and perpendicular lamellae at the free surface. The calculations confirm that the intermediate state is unstable to thickness fluctuations, thereby driving the spinodal growth of surface structures.

Applying the chain formation model to magnetic properties of aggregated ferrofluids

The magnetization properties of aggregated ferrofluids are calculated by combining the chain formation model developed by Zubarev with the modified mean-field theory. Using moderate assumptions for the inter- and intrachain interactions we obtain expressions for the magnetization and initial susceptibility. When comparing the results of our theory to molecular dynamics simulations of the same model we find that at large dipolar couplings (lambda>3) the chain formation model appears to give better predictions than other analytical approaches. This supports the idea that chain formation is an important structural ingredient of strongly interacting dipolar particles.

Polymer brushes under shear: molecular dynamics simulations compared to experiments

Surfaces coated with polymer brushes in a good solvent are known to exhibit excellent tribological properties. We have performed coarse-grained equilibrium and nonequilibrium molecular dynamics (MD) simulations to investigate dextran polymer brushes in an aqueous environment in molecular detail. In a first step, we determined simulation parameters and units by matching experimental results for a single dextran chain. Analyzing this model when applied to a multichain system, density profiles of end-tethered polymer brushes obtained from equilibrium MD simulations compare very well with expectations based on self-consistent field theory. Simulation results were further validated against and correlated with available experimental results. The simulated compression curves (normal force as a function of surface separation) compare successfully with results obtained with a surface forces apparatus. Shear stress (friction) obtained via nonequilibrium MD is contrasted with nanoscale friction studies employing colloidal-probe lateral force microscopy. We find good agreement in the hydrodynamic regime and explain the observed leveling-off of the friction forces in the boundary regime by means of an effective polymer–wall attraction.

Risk taking of executives under different incentive contracts: experimental evidence

Classic financial agency theory recommends compensation through stock options rather than shares to counteract excessive risk aversion in agents. In a setting where any kind of risk taking is suboptimal for shareholders, we show that excessive risk taking may occur for one of two reasons: risk preferences or incentives. Even when compensated through restricted company stock, experimental CEOs take large amounts of excessive risk. This contradicts classical financial theory, but can be explained through risk preferences that are not uniform over the probability and outcome spaces, and in particular, risk seeking for small probability gains and large probability losses. Compensation through options further increases risk taking as expected. We show that this effect is driven mainly by the personal asset position of the experimental CEO, thus having deleterious effects on company performance.

Transcendental Brauer groups of products of CM elliptic curves

Let L be a number field and let E/L be an elliptic curve with complex multiplication by the ring of integers O_K of an imaginary quadratic field K. We use class field theory and results of Skorobogatov and Zarhin to compute the transcendental part of the Brauer group of the abelian surface ExE. The results for the odd order torsion also apply to the Brauer group of the K3 surface Kum(ExE). We describe explicitly the elliptic curves E/Q with complex multiplication by O_K such that the Brauer group of ExE contains a transcendental element of odd order. We show that such an element gives rise to a Brauer-Manin obstruction to weak approximation on Kum(ExE), while there is no obstruction coming from the algebraic part of the Brauer group.

Equilibrium behavior of symmetric ABA triblock copolymer melts

Melts of ABA triblock copolymer molecules with identical end blocks are examined using self-consistent field theory (SCFT). Phase diagrams are calculated and compared with those of homologous AB diblock copolymers formed by snipping the triblocks in half. This creates additional end segments which decreases the degree of segregation. Consequently, triblock melts remain ordered to higher temperatures than their diblock counterparts. We also find that middle-block domains are easier to stretch than end-block domains. As a result, domain spacings are slightly larger, the complex phase regions are shifted towards smaller A-segment compositions, and the perforated-lamellar phase becomes more metastable in triblock melts as compared to diblock melts. Although triblock and diblock melts exhibit very similar phase behavior, their mechanical properties can differ substantially due to triblock copolymers that bridge between otherwise disconnected A domains. We evaluate the bridging fraction for lamellar, cylindrical, and spherical morphologies to be about 40%–45%, 60%–65%, and 75%–80%, respectively. These fractions only depend weakly on the degree of segregation and the copolymer composition.

Effective interaction between monolayers of block copolymer compatiblizer in a polymer blend

The stability of ternary blends of two immiscible homopolymers and a block copolymer compatiblizer depends crucially on the effective interaction between the copolymermonolayers that form between the unlike homopolymer domains. Here, the interaction is calculated for blends involving A and B homopolymers of equal size with ABABdiblock copolymers of symmetric composition using both self-consistent field theory (SCFT) and strong-segregation theory (SST). If the homopolymers are larger than the copolymer molecules, an attractive interaction is predicted which would destroy the blend. This conclusion coupled with considerations regarding the elastic properties of the monolayer suggests that the optimum size of the homopolymer molecules is about 80% that of the copolymer molecule. A detailed examination of the theory demonstrates that the attraction results from the configurational entropy loss of the homopolymer molecules trapped between the copolymermonolayers. We conclude by suggesting how the monolayers can be altered in order to suppress this attraction and thus improve compatiblization.

Instabilities in thermal gravity with a cosmological constant

It is shown that in quantum gravity at finite temperature, the effective potential evaluated in the tadpole approximation can have a local minimum below a certain critical temperature. However, when the leading higher order thermal loop corrections are included, one finds that no static solution exists at high temperature. (C) 2008 Elsevier B.V. All rights reserved.

Critical behavior of a contact process with aperiodic transition rates

We performed Monte Carlo simulations to investigate the steady-state critical behavior of a one-dimensional contact process with an aperiodic distribution of rates of transition. As in the presence of randomness, spatial fluctuations can lead to changes of critical behavior. For sufficiently weak fluctuations, we give numerical evidence to show that there is no departure from the universal critical behavior of the underlying uniform model. For strong spatial fluctuations, the analysis of the data indicates a change of critical universality class.

Finite temperature effective actions

We present, from first principles, a direct method for evaluating the exact fermion propagator in the presence of a general background held at finite temperature, which can be used to determine the finite temperature effective action for the system. As applications, we determine the complete one loop finite temperature effective actions for (0 + 1)-dimensional QED as well as the Schwinger model. These effective actions, which are derived in the real time (closed time path) formalism, generate systematically all the Feynman amplitudes calculated in thermal perturbation theory and also show that the retarded (advanced) amplitudes vanish in these theories. (C) 2009 Elsevier B.V. All rights reserved.

A Precise Formulation of the Third Law of Thermodynamics

The third law of thermodynamics is formulated precisely: all points of the state space of zero temperature I""(0) are physically adiabatically inaccessible from the state space of a simple system. In addition to implying the unattainability of absolute zero in finite time (or ""by a finite number of operations""), it admits as corollary, under a continuity assumption, that all points of I""(0) are adiabatically equivalent. We argue that the third law is universally valid for all macroscopic systems which obey the laws of quantum mechanics and/or quantum field theory. We also briefly discuss why a precise formulation of the third law for black holes remains an open problem.