Metastable chaos in the ammonia ring laser

We report experimental studies of metastable chaos in the far-infrared ammonia ring: laser. When the laser pump power is switched from above chaos threshold to slightly below, chaotic intensity pulsations continue for a varying time afterward before decaying to either periodic or cw emission. The behavior is in good qualitative agreement with that predicted by the Lorenz equations, previously used to describe this laser. The statistical distribution of the duration of the chaotic transient is measured and shown to be in excellent agreement with the Lorenz equations in showing a modified exponential distribution. We also give a brief numerical analysis and graphical visualization of the Lorenz equations in phase space illustrating the boundary between the metastable chaotic and the stable fixed point basins of attraction. This provides an intuitive understanding of the metastable dynamics of the Lorenz equations and the experimental system.

Spatiotemporal solitons in multidimensional optical media with a quadratic nonlinearity

We consider solutions to the second-harmonic generation equations in two-and three-dimensional dispersive media in the form of solitons localized in space and time. As is known, collapse does not take place in these models, which is why the solitons may be stable. The general solution is obtained in an approximate analytical form by means of a variational approach, which also allows the stability of the solutions to be predicted. Then, we directly simulate the two-dimensional case, taking the initial configuration as suggested by the variational approximation. We thus demonstrate that spatiotemporal solitons indeed exist and are stable. Furthermore, they are not, in the general case, equivalent to the previously known cylindrical spatial solitons. Direct simulations generate solitons with some internal oscillations. However, these oscillations neither grow nor do they exhibit any significant radiative damping. Numerical solutions of the stationary version of the equations produce the same solitons in their unperturbed form, i.e., without internal oscillations. Strictly stable solitons exist only if the system has anomalous dispersion at both the fundamental harmonic and second harmonic (SH), including the case of zero dispersion at SH. Quasistationary solitons, decaying extremely slowly into radiation, are found in the presence of weak normal dispersion at the second-harmonic frequency.

Controlling the complex Lorenz equations by modulation

We demonstrate that a system obeying the complex Lorenz equations in the deep chaotic regime can be controlled to periodic behavior by applying a modulation to the pump parameter. For arbitrary modulation frequency and amplitude there is no obvious simplification of the dynamics. However, we find that there are numerous windows where the chaotic system has been controlled to different periodic behaviors. The widths of these windows in parameter space are narrow, and the positions are related to the ratio of the modulation frequency of the pump to the average pulsation frequency of the output variable. These results are in good agreement with observations previously made in a far-infrared laser system.

Biomass equations for forest regrowth in the eastern Amazon using randomized branch sampling

Forest regrowth occupies an extensive and increasing area in the Amazon basin, but accurate assessment of the impact of regrowth on carbon and nutrient cycles has been hampered by a paucity of available allometric equations. We develop pooled and species-specific equations for total aboveground biomass for a study site in the eastern Amazon that had been abandoned for 15 years. Field work was conducted using randomized branch sampling, a rapid technique that has seen little use in tropical forests. High consistency of sample paths in randomized branch sampling, as measured by the standard error of individual paths (14%), suggests the method may provide substantial efficiencies when compared to traditional procedures. The best fitting equations in this study used the traditional form Y=a×DBHb, where Y is biomass, DBH is diameter at breast height, and a and b are both species-specific parameters. Species-specific equations of the form Y=a(BA×H), where Y is biomass, BA is tree basal area, H is tree height, and a is a species-specific parameter, fit almost as well. Comparison with previously published equations indicated errors from -33% to +29% would have occurred using off-site relationships. We also present equations for stemwood, twigs, and foliage as biomass components.

Sex-Specific Equations to Estimate Maximum Oxygen Uptake in Cycle Ergometry

AbstractBackground:Aerobic fitness, assessed by measuring VO2max in maximum cardiopulmonary exercise testing (CPX) or by estimating VO2max through the use of equations in exercise testing, is a predictor of mortality. However, the error resulting from this estimate in a given individual can be high, affecting clinical decisions.Objective:To determine the error of estimate of VO2max in cycle ergometry in a population attending clinical exercise testing laboratories, and to propose sex-specific equations to minimize that error.Methods:This study assessed 1715 adults (18 to 91 years, 68% men) undertaking maximum CPX in a lower limbs cycle ergometer (LLCE) with ramp protocol. The percentage error (E%) between measured VO2max and that estimated from the modified ACSM equation (Lang et al. MSSE, 1992) was calculated. Then, estimation equations were developed: 1) for all the population tested (C-GENERAL); and 2) separately by sex (C-MEN and C-WOMEN).Results:Measured VO2max was higher in men than in WOMEN: -29.4 ± 10.5 and 24.2 ± 9.2 mL.(kg.min)-1 (p < 0.01). The equations for estimating VO2max [in mL.(kg.min)-1] were: C-GENERAL = [final workload (W)/body weight (kg)] x 10.483 + 7; C-MEN = [final workload (W)/body weight (kg)] x 10.791 + 7; and C-WOMEN = [final workload (W)/body weight (kg)] x 9.820 + 7. The E% for MEN was: -3.4 ± 13.4% (modified ACSM); 1.2 ± 13.2% (C-GENERAL); and -0.9 ± 13.4% (C-MEN) (p < 0.01). For WOMEN: -14.7 ± 17.4% (modified ACSM); -6.3 ± 16.5% (C-GENERAL); and -1.7 ± 16.2% (C-WOMEN) (p < 0.01).Conclusion:The error of estimate of VO2max by use of sex-specific equations was reduced, but not eliminated, in exercise tests on LLCE.

Developing discontinuous Galerkin methods for the resolution of the incompressible Navier-Stokes equations

Informe de investigación elaborado a partir de una estancia en el Laboratorio de Diseño Computacional en Aeroespacial en el Massachusetts Institute of Technology (MIT), Estados Unidos, entre noviembre de 2006 y agosto de 2007. La aerodinámica es una rama de la dinámica de fluidos referida al estudio de los movimientos de los líquidos o gases, cuya meta principal es predecir las fuerzas aerodinámicas en un avión o cualquier tipo de vehículo, incluyendo los automóviles. Las ecuaciones de Navier-Stokes representan un estado dinámico del equilibrio de las fuerzas que actúan en cualquier región dada del fluido. Son uno de los sistemas de ecuaciones más útiles porque describen la física de una gran cantidad de fenómenos como corrientes del océano, flujos alrededor de una superficie de sustentación, etc. En el contexto de una tesis doctoral, se está estudiando un flujo viscoso e incompresible, solucionando las ecuaciones de Navier- Stokes incompresibles de una manera eficiente. Durante la estancia en el MIT, se ha utilizado un método de Galerkin discontinuo para solucionar las ecuaciones de Navier-Stokes incompresibles usando, o bien un parámetro de penalti para asegurar la continuidad de los flujos entre elementos, o bien un método de Galerkin discontinuo compacto. Ambos métodos han dado buenos resultados y varios ejemplos numéricos se han simulado para validar el buen comportamiento de los métodos desarrollados. También se han estudiado elementos particulares, los elementos de Raviart y Thomas, que se podrían utilizar en una formulación mixta para obtener un algoritmo eficiente para solucionar problemas numéricos complejos.

Brownian motion in short range random potentials

A numerical study of Brownian motion of noninteracting particles in random potentials is presented. The dynamics are modeled by Langevin equations in the high friction limit. The random potentials are Gaussian distributed and short ranged. The simulations are performed in one and two dimensions. Different dynamical regimes are found and explained. Effective subdiffusive exponents are obtained and commented on.

Stochastic semiclassical equations for weakly inhomogeneous cosmologies

Semiclassical Einstein-Langevin equations for arbitrary small metric perturbations conformally coupled to a massless quantum scalar field in a spatially flat cosmological background are derived. Use is made of the fact that for this problem the in-in or closed time path effective action is simply related to the Feynman-Vernon influence functional which describes the effect of the ``environment,'' the quantum field which is coarse grained here, on the ``system,'' the gravitational field which is the field of interest. This leads to identify the dissipation and noise kernels in the in-in effective action, and to derive a fluctuation-dissipation relation. A tensorial Gaussian stochastic source which couples to the Weyl tensor of the spacetime metric is seen to modify the usual semiclassical equations which can be veiwed now as mean field equsations. As a simple application we derive the correlation functions of the stochastic metric fluctuations produced in a flat spacetime with small metric perturbations due to the quantum fluctuations of the matter field coupled to these perturbations.

Multisoliton solutions to Einstein's equations

We discuss a multisoliton solution to Einsteins equations in vacuum. The solution is interpreted as many gravitational solitons propagating and colliding on a homogeneous cosmological background. Following a previous letter, we characterize the solitons by their localizability and by their peculiar properties under collisions. Furthermore, we define an associated frame-dependent velocity field which illustrates the solitonic character of these gravitational solitons in the classical sense.

Soliton collision in general relativity

An exact solution of the Einstein equations in vacuum representing two pairs of gravitational solitons propagating on an expanding universe is given and studied. It is suggested that the solitons evolve from quasiparticles to pure gravitational waves. Two of the four solitons collide and the focusing produced on null rays is studied. Although the spacetime following the collision is highly distorted, null rays do not focus to a singularity.

The problem of physical coordinates in predictive Hamiltonian systems

In the Hamiltonian formulation of predictive relativistic systems, the canonical coordinates cannot be the physical positions. The relation between them is given by the individuality differential equations. However, due to the arbitrariness in the choice of Cauchy data, there is a wide family of solutions for these equations. In general, those solutions do not satisfy the condition of constancy of velocities moduli, and therefore we have to reparametrize the world lines into the proper time. We derive here a condition on the Cauchy data for the individuality equations which ensures the constancy of the velocities moduli and makes the reparametrization unnecessary.

Dynamics of orientational fluctuations in a homeotropic Fréedericksz transition under a rotating magnetic field

We study the problem of the Fréedericksz transition under a rotating magnetic field by using a dynamical model which incorporates thermal fluctuations into the whole set of nematodynamic equations. In contrast to other geometries, nonuniform textures in the plane of the sample do not appear favored. The proper consideration of thermal noise enables us to describe the dynamics of orientational fluctuations both below and above the shifted instability.

Planar Topological Defects in Unconventional Superconductors

In this work, we consider the properties of planar topological defects in unconventional superconductors. Specifically, we calculate microscopically the interaction energy of domain walls separating degenerate ground states in a chiral p-wave fermionic superfluid. The interaction is mediated by the quasiparticles experiencing Andreev scattering at the domain walls. As a by-product, we derive a useful general expression for the free energy of an arbitrary nonuniform texture of the order parameter in terms of the quasiparticle scattering matrix. The thesis is structured as follows. We begin with a historical review of the theories of superconductivity (Sec. 1.1), which led the way to the celebrated Bardeen-Cooper- Schrieffer (BCS) theory (Sec. 1.3). Then we proceed to the treatment of superconductors with so-called "unconventional pairing" in Sec. 1.4, and in Sec. 1.5 we introduce the specific case of chiral p-wave superconductivity. After introducing in Sec. 2 the domain wall (DW) model that will be considered throughout the work, we derive the Bogoliubov-de Gennes (BdG) equations in Sec. 3.1, which determine the quasiparticle excitation spectrum for a nonuniform superconductor. In this work, we use the semiclassical (Andreev) approximation, and solve the Andreev equations (which are a particular case of the BdG equations) in Sec. 4 to determine the quasiparticle spectrum for both the single- and two-DW textures. The Andreev equations are derived in Sec. 3.2, and the formal properties of the Andreev scattering coefficients are discussed in the following subsection. In Sec. 5, we use the transfer matrix method to relate the interaction energy of the DWs to the scattering matrix of the Bogoliubov quasiparticles. This facilitates the derivation of an analytical expression for the interaction energy between the two DWs in Sec. 5.3. Finally, to illustrate the general applicability our method, we apply it in Sec. 6 to the interaction between phase solitons in a two-band s-wave superconductor.