Radiative forcing and climate metrics for ozone precursor emissions: the impact of multi-model averaging

Multi-model ensembles are frequently used to assess understanding of the response of ozone and methane lifetime to changes in emissions of ozone precursors such as NOx, VOCs (volatile organic compounds) and CO. When these ozone changes are used to calculate radiative forcing (RF) (and climate metrics such as the global warming potential (GWP) and global temperature-change potential (GTP)) there is a methodological choice, determined partly by the available computing resources, as to whether the mean ozone (and methane) concentration changes are input to the radiation code, or whether each model's ozone and methane changes are used as input, with the average RF computed from the individual model RFs. We use data from the Task Force on Hemispheric Transport of Air Pollution source–receptor global chemical transport model ensemble to assess the impact of this choice for emission changes in four regions (East Asia, Europe, North America and South Asia). We conclude that using the multi-model mean ozone and methane responses is accurate for calculating the mean RF, with differences up to 0.6% for CO, 0.7% for VOCs and 2% for NOx. Differences of up to 60% for NOx 7% for VOCs and 3% for CO are introduced into the 20 year GWP. The differences for the 20 year GTP are smaller than for the GWP for NOx, and similar for the other species. However, estimates of the standard deviation calculated from the ensemble-mean input fields (where the standard deviation at each point on the model grid is added to or subtracted from the mean field) are almost always substantially larger in RF, GWP and GTP metrics than the true standard deviation, and can be larger than the model range for short-lived ozone RF, and for the 20 and 100 year GWP and 100 year GTP. The order of averaging has most impact on the metrics for NOx, as the net values for these quantities is the residual of the sum of terms of opposing signs. For example, the standard deviation for the 20 year GWP is 2–3 times larger using the ensemble-mean fields than using the individual models to calculate the RF. The source of this effect is largely due to the construction of the input ozone fields, which overestimate the true ensemble spread. Hence, while the average of multi-model fields are normally appropriate for calculating mean RF, GWP and GTP, they are not a reliable method for calculating the uncertainty in these fields, and in general overestimate the uncertainty.

VELOCITY FIELD OF COMPRESSIBLE MAGNETOHYDRODYNAMIC TURBULENCE: WAVELET DECOMPOSITION AND MODE SCALINGS

We study compressible magnetohydrodynamic turbulence, which holds the key to many astrophysical processes, including star formation and cosmic-ray propagation. To account for the variations of the magnetic field in the strongly turbulent fluid, we use wavelet decomposition of the turbulent velocity field into Alfven, slow, and fast modes, which presents an extension of the Cho & Lazarian decomposition approach based on Fourier transforms. The wavelets allow us to follow the variations of the local direction of the magnetic field and therefore improve the quality of the decomposition compared to the Fourier transforms, which are done in the mean field reference frame. For each resulting component, we calculate the spectra and two-point statistics such as longitudinal and transverse structure functions as well as higher order intermittency statistics. In addition, we perform a Helmholtz-Hodge decomposition of the velocity field into incompressible and compressible parts and analyze these components. We find that the turbulence intermittency is different for different components, and we show that the intermittency statistics depend on whether the phenomenon was studied in the global reference frame related to the mean magnetic field or in the frame defined by the local magnetic field. The dependencies of the measures we obtained are different for different components of the velocity; for instance, we show that while the Alfven mode intermittency changes marginally with the Mach number, the intermittency of the fast mode is substantially affected by the change.

3-D reconstruction by extended depth-of-field in failure analysis - Case study II: Fractal analysis of interlaminar fracture in carbon/epoxy composites

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

MRT letter: Extended depth from focus reconstruction method for stretch zone measurement in 15-5PH steel

The stretch zone width (SZW) data for 15-5PH steel CTOD specimens fractured at -150 degrees C to + 23 degrees C temperature were measured based on focused images and 3D maps obtained by extended depth-of-field reconstruction from light microscopy (LM) image stacks. This LM-based method, with a larger lateral resolution, seems to be as effective for quantitative analysis of SZW as scanning electron microscopy (SEM) or confocal scanning laser microscopy (CSLM), permitting to clearly identify stretch zone boundaries. Despite the worst sharpness of focused images, a robust linear correlation was established to fracture toughness (KC) and SZW data for the 15-5PH steel tested specimens, measured at their center region. The method is an alternative to evaluate the boundaries of stretched zones, at a lower cost of implementation and training, since topographic data from elevation maps can be associated with reconstructed image, which summarizes the original contrast and brightness information. Finally, the extended depth-of-field method is presented here as a valuable tool for failure analysis, as a cheaper alternative to investigate rough surfaces or fracture, compared to scanning electron or confocal light microscopes. Microsc. Res. Tech. 75:11551158, 2012. (C) 2012 Wiley Periodicals, Inc.

Black soliton in a quasi-one-dimensional trapped fermion-fermion mixture

Employing a time dependent mean-field-hydrodynamic model we study the generation of black solitons in a degenerate fermion-fermion mixture in a cigar-shaped geometry using variational and numerical solutions. The black soliton is found to be the first stationary vibrational excitation of the system and is considered to be a nonlinear continuation of the vibrational excitation of the harmonic oscillator state. We illustrate the stationary nature of the black soliton, by studying different perturbations on it after its formation.

Bound states of attractive Bose-Einstein condensates in shallow traps in two and three dimensions

Using variational and numerical solutions of the mean-field Gross-Pitaevskii equation for attractive interaction (with cubic or Kerr nonlinearity), we show that a stable bound state can appear in a Bose-Einstein condensate (BEC) in a localized exponentially screened radially symmetric harmonic potential well in two and three dimensions. We also consider an axially symmetric configuration with zero axial trap and a exponentially screened radial trap so that the resulting bound state can freely move along the axial direction like a soliton. The binding of the present states in shallow wells is mostly due to the nonlinear interaction with the trap playing a minor role. Hence, these BEC states are more suitable to study the effect of the nonlinear force on the dynamics. We illustrate the highly nonlinear nature of breathing oscillations of these states. Such bound states could be created in BECs and studied in the laboratory with present knowhow.

Dissipation-managed soliton in a quasi-one-dimensional Bose-Einstein condensate

We use the time-dependent mean-field Cross-Pitaevskii equation to study the formation of a dynamically-stabilized dissipation managed bright soliton in a quasi-one dimensional Bose-Einstein condensate (BEC). Because of three-body recombination of bosonic atoms to molecules, atoms are lost (dissipated) from a BEC. Such dissipation leads to the decay of a BEC soliton. We demonstrate by a perturbation procedure that an alimentation of atoms from an external source to the BEC may compensate for the dissipation loss and lead to a dynamically-stabilized soliton. The result of the analytical perturbation method is in excellent agreement with mean-field numerics. It seems possible to obtain such a dynamically stabilized BEC soliton without dissipation in laboratory.

Formation of bright solitons and soliton trains in a fermion-fermion mixture by modulational instability

We employ a time- dependent mean- field- hydrodynamic model to study the generation of bright solitons in a degenerate fermion - fermion mixture in a cigar- shaped geometry using variational and numerical methods. Due to a strong Pauli- blocking repulsion among identical spin- polarized fermions at short distances there cannot be bright solitons for repulsive interspecies interactions. Employing a linear stability analysis we demonstrate the formation of stable solitons due to modulational instability of a constant-amplitude solution of the model equations for a sufficiently attractive interspecies interaction. We perform a numerical stability analysis of these solitons and also demonstrate the formation of soliton trains by jumping the effective interspecies interaction from repulsive to attractive. These fermionic solitons can be formed and studied in laboratory with present technology.

Tightly bound gap solitons in a Fermi gas

Within the framework of the mean-field hydrodynamic model of a degenerate Fermi gas ( DFG), we study, by means of numerical methods and variational approximation ( VA), the formation of fundamental gap solitons ( FGSs) in a DFG ( or in a BCS superfluid generated by weak interaction between spin- up and spin- down fermions), which is trapped in a periodic optical- lattice ( OL) potential. An effectively one- dimensional ( 1D) con. guration is considered, assuming strong transverse confinement; in parallel, a proper 1D model of the DFG ( which amounts to the known quintic equation for the Tonks- Girardeau gas in the OL) is considered too. The FGSs found in the first two bandgaps of the OL- induced spectrum ( unless they are very close to edges of the gaps) feature a ( tightly bound) shape, being essentially confined to a single cell of the OL. In the second bandgap, we also find antisymmetric tightly bound subfundamental solitons ( SFSs), with zero at the midpoint. The SFSs are also confined to a single cell of the OL, but, unlike the FGSs, they are unstable. The predicted solitons, consisting of similar to 10(4) - 10(5) atoms, can be created by available experimental techniques in the DFG of Li-6 atoms.

Stabilization of bright solitons and vortex solitons in a trapless three-dimensional Bose-Einstein condensate by temporal modulation of the scattering length

Using variational and numerical solutions of the mean-field Gross-Pitaevskii equation we show that a bright soliton can be stabilized in a trapless three-dimensional attractive Bose-Einstein condensate (BEC) by a rapid periodic temporal modulation of scattering length alone by using a Feshbach resonance. This scheme also stabilizes a rotating vortex soliton in two dimensions. Apart from possible experimental application in BEC, the present study suggests that the spatiotemporal solitons of nonlinear optics in three dimensions can also be stabilized in a layered Kerr medium with sign-changing nonlinearity along the propagation direction.

Bright solitons and soliton trains in a fermion-fermion mixture

We use a time-dependent dynamical mean-field-hydrodynamic model to predict and study bright solitons in a degenerate fermion-fermion mixture in a quasi-one-dimensional cigar-shaped geometry using variational and numerical methods. Due to a strong Pauli-blocking repulsion among identical spin-polarized fermions at short distances there cannot be bright solitons for repulsive interspecies fermion-fermion interactions. However, stable bright solitons can be formed for a sufficiently attractive interspecies interaction. We perform a numerical stability analysis of these solitons and also demonstrate the formation of soliton trains. These fermionic solitons can be formed and studied in laboratory with present technology.

Miscibility in a degenerate fermionic mixture induced by linear coupling

We consider a one-dimensional mean-field-hydrodynamic model of a two-component degenerate Fermi gas in an external trap, each component representing a spin state of the same atom. We demonstrate that the interconversion between them (linear coupling), imposed by a resonant electromagnetic wave, transforms the immiscible binary gas into a miscible state, if the coupling constant, kappa, exceeds a critical value, kappa(cr). The effect is predicted in a variational approximation, and confirmed by numerical solutions. Unlike the recently studied model of a binary Bose-Einsten condensate with the linear coupling, the components in the immiscible phase of the binary fermion mixture never fill two separated domains with a wall between them, but rather form antilocked (pi-phase-shifted) density waves. Another difference from the bosonic mixture is spontaneous breaking of symmetry between the two components in terms of the numbers of atoms in them, N(1) and N(2). The latter effect is characterized by the parameter nu equivalent to(N(1)-N(2))/(N(1)+N(2)) (only N(1)+N(2) is a conserved quantity), the onset of miscibility at kappa >=kappa(cr) meaning a transition to nu equivalent to 0. At kappa

Gap solitons in a model of a superfluid fermion gas in optical lattices

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A aproximação de campo médio de Bethe-Peierls

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

ON THE INDEPENDENT-PARTICLE APPROXIMATION OF GAUGE-THEORIES - A SIMPLE EXAMPLE

In this work the independent particle model formulation is studied as a mean-field approximation of gauge theories using the path integral approach in the framework of quantum electrodynamics in 1 + 1 dimensions. It is shown how a mean-field approximation scheme can be applied to fit an effective potential to an independent particle model, building a straightforward relation between the model and the associated gauge field theory. An example is made considering the problem of massive Dirac fermions on a line, the so called massive Schwinger model. An interesting result is found, indicating a behaviour of screening of the charges in the relativistic limit of strong coupling. A forthcoming application of the method developed to confining potentials in independent quark models for QCD is in view and is briefly discussed.