Magnetohydrodynamic turbulence

We simulate incompressible, MHD turbulence using a pseudo-spectral code. Our major conclusions are as follows.

1) MHD turbulence is most conveniently described in terms of counter propagating shear Alfvén and slow waves. Shear Alfvén waves control the cascade dynamics. Slow waves play a passive role and adopt the spectrum set by the shear Alfvén waves. Cascades composed entirely of shear Alfvén waves do not generate a significant measure of slow waves.

2) MHD turbulence is anisotropic with energy cascading more rapidly along k_⊥ than along k_∥, where k_⊥ and k_∥ refer to wavevector components perpendicular and parallel to the local magnetic field. Anisotropy increases with increasing k_⊥ such that excited modes are confined inside a cone bounded by k_∥ ∝ k^γ_⊥ where γ less than 1. The opening angle of the cone, θ(k_⊥) ∝ k^-(1-γ)_⊥, defines the scale dependent anisotropy.

3) MHD turbulence is generically strong in the sense that the waves which comprise it suffer order unity distortions on timescales comparable to their periods. Nevertheless, turbulent fluctuations are small deep inside the inertial range. Their energy density is less than that of the background field by a factor θ² (k_⊥)≪1.

4) MHD cascades are best understood geometrically. Wave packets suffer distortions as they move along magnetic field lines perturbed by counter propagating waves. Field lines perturbed by unidirectional waves map planes perpendicular to the local field into each other. Shear Alfvén waves are responsible for the mapping's shear and slow waves for its dilatation. The amplitude of the former exceeds that of the latter by 1/θ(k_⊥) which accounts for dominance of the shear Alfvén waves in controlling the cascade dynamics.

5) Passive scalars mixed by MHD turbulence adopt the same power spectrum as the velocity and magnetic field perturbations.

6) Decaying MHD turbulence is unstable to an increase of the imbalance between the flux of waves propagating in opposite directions along the magnetic field. Forced MHD turbulence displays order unity fluctuations with respect to the balanced state if excited at low k by δ(t) correlated forcing. It appears to be statistically stable to the unlimited growth of imbalance.

7) Gradients of the dynamic variables are focused into sheets aligned with the magnetic field whose thickness is comparable to the dissipation scale. Sheets formed by oppositely directed waves are uncorrelated. We suspect that these are vortex sheets which the mean magnetic field prevents from rolling up.

8) Items (1)-(5) lend support to the model of strong MHD turbulence put forth by Goldreich and Sridhar (1995, 1997). Results from our simulations are also consistent with the GS prediction γ = 2/3. The sole not able discrepancy is that the 1D power law spectra, E(k_⊥) ∝ k^-∝_⊥, determined from our simulations exhibit ∝ ≈ 3/2, whereas the GS model predicts ∝ = 5/3.

The Madden-Julian Oscillation : observation, modeling, and theory

The Madden-Julian Oscillation (MJO) is a pattern of intense rainfall and associated planetary-scale circulations in the tropical atmosphere, with a recurrence interval of 30-90 days. Although the MJO was first discovered 40 years ago, it is still a challenge to simulate the MJO in general circulation models (GCMs), and even with simple models it is difficult to agree on the basic mechanisms. This deficiency is mainly due to our poor understanding of moist convection—deep cumulus clouds and thunderstorms, which occur at scales that are smaller than the resolution elements of the GCMs. Moist convection is the most important mechanism for transporting energy from the ocean to the atmosphere. Success in simulating the MJO will improve our understanding of moist convection and thereby improve weather and climate forecasting.

We address this fundamental subject by analyzing observational datasets, constructing a hierarchy of numerical models, and developing theories. Parameters of the models are taken from observation, and the simulated MJO fits the data without further adjustments. The major findings include: 1) the MJO may be an ensemble of convection events linked together by small-scale high-frequency inertia-gravity waves; 2) the eastward propagation of the MJO is determined by the difference between the eastward and westward phase speeds of the waves; 3) the planetary scale of the MJO is the length over which temperature anomalies can be effectively smoothed by gravity waves; 4) the strength of the MJO increases with the typical strength of convection, which increases in a warming climate; 5) the horizontal scale of the MJO increases with the spatial frequency of convection; and 6) triggered convection, where potential energy accumulates until a threshold is reached, is important in simulating the MJO. Our findings challenge previous paradigms, which consider the MJO as a large-scale mode, and point to ways for improving the climate models.

Anomaly in photoelectron angular distributions

The recently observed anomaly in photoelectron angular distributions (PADs), the disappearance of the main lobes of PADs which should be usually in the direction of laser polarization, is reinterpreted as a minimum of generalized Bessel functions in the laser-polarization direction with the theory of nonperturbative quantum electrodynamics. The reinterpretation has no artificial fitting parameters and explains more features of the experimentally observed PADs, in contrast to the existing interpretation in which the anomaly is interpreted as a quantum interference of angular momentum partial waves. Some hierarchy anomalies are predicted for further experimental observations.

The power of quantum Fourier sampling

How powerful are Quantum Computers? Despite the prevailing belief that Quantum Computers are more powerful than their classical counterparts, this remains a conjecture backed by little formal evidence. Shor's famous factoring algorithm [Shor97] gives an example of a problem that can be solved efficiently on a quantum computer with no known efficient classical algorithm. Factoring, however, is unlikely to be NP-Hard, meaning that few unexpected formal consequences would arise, should such a classical algorithm be discovered. Could it then be the case that any quantum algorithm can be simulated efficiently classically? Likewise, could it be the case that Quantum Computers can quickly solve problems much harder than factoring? If so, where does this power come from, and what classical computational resources do we need to solve the hardest problems for which there exist efficient quantum algorithms?

We make progress toward understanding these questions through studying the relationship between classical nondeterminism and quantum computing. In particular, is there a problem that can be solved efficiently on a Quantum Computer that cannot be efficiently solved using nondeterminism? In this thesis we address this problem from the perspective of sampling problems. Namely, we give evidence that approximately sampling the Quantum Fourier Transform of an efficiently computable function, while easy quantumly, is hard for any classical machine in the Polynomial Time Hierarchy. In particular, we prove the existence of a class of distributions that can be sampled efficiently by a Quantum Computer, that likely cannot be approximately sampled in randomized polynomial time with an oracle for the Polynomial Time Hierarchy.

Our work complements and generalizes the evidence given in Aaronson and Arkhipov's work [AA2013] where a different distribution with the same computational properties was given. Our result is more general than theirs, but requires a more powerful quantum sampler.

Efficient methods for stochastic optimal control

The Hamilton Jacobi Bellman (HJB) equation is central to stochastic optimal control (SOC) theory, yielding the optimal solution to general problems specified by known dynamics and a specified cost functional. Given the assumption of quadratic cost on the control input, it is well known that the HJB reduces to a particular partial differential equation (PDE). While powerful, this reduction is not commonly used as the PDE is of second order, is nonlinear, and examples exist where the problem may not have a solution in a classical sense. Furthermore, each state of the system appears as another dimension of the PDE, giving rise to the curse of dimensionality. Since the number of degrees of freedom required to solve the optimal control problem grows exponentially with dimension, the problem becomes intractable for systems with all but modest dimension.

In the last decade researchers have found that under certain, fairly non-restrictive structural assumptions, the HJB may be transformed into a linear PDE, with an interesting analogue in the discretized domain of Markov Decision Processes (MDP). The work presented in this thesis uses the linearity of this particular form of the HJB PDE to push the computational boundaries of stochastic optimal control.

This is done by crafting together previously disjoint lines of research in computation. The first of these is the use of Sum of Squares (SOS) techniques for synthesis of control policies. A candidate polynomial with variable coefficients is proposed as the solution to the stochastic optimal control problem. An SOS relaxation is then taken to the partial differential constraints, leading to a hierarchy of semidefinite relaxations with improving sub-optimality gap. The resulting approximate solutions are shown to be guaranteed over- and under-approximations for the optimal value function. It is shown that these results extend to arbitrary parabolic and elliptic PDEs, yielding a novel method for Uncertainty Quantification (UQ) of systems governed by partial differential constraints. Domain decomposition techniques are also made available, allowing for such problems to be solved via parallelization and low-order polynomials.

The optimization-based SOS technique is then contrasted with the Separated Representation (SR) approach from the applied mathematics community. The technique allows for systems of equations to be solved through a low-rank decomposition that results in algorithms that scale linearly with dimensionality. Its application in stochastic optimal control allows for previously uncomputable problems to be solved quickly, scaling to such complex systems as the Quadcopter and VTOL aircraft. This technique may be combined with the SOS approach, yielding not only a numerical technique, but also an analytical one that allows for entirely new classes of systems to be studied and for stability properties to be guaranteed.

The analysis of the linear HJB is completed by the study of its implications in application. It is shown that the HJB and a popular technique in robotics, the use of navigation functions, sit on opposite ends of a spectrum of optimization problems, upon which tradeoffs may be made in problem complexity. Analytical solutions to the HJB in these settings are available in simplified domains, yielding guidance towards optimality for approximation schemes. Finally, the use of HJB equations in temporal multi-task planning problems is investigated. It is demonstrated that such problems are reducible to a sequence of SOC problems linked via boundary conditions. The linearity of the PDE allows us to pre-compute control policy primitives and then compose them, at essentially zero cost, to satisfy a complex temporal logic specification.

Velocity resolved - scalar modeled simulations of high Schmidt number turbulent transport

The objective of this thesis is to develop a framework to conduct velocity resolved - scalar modeled (VR-SM) simulations, which will enable accurate simulations at higher Reynolds and Schmidt (Sc) numbers than are currently feasible. The framework established will serve as a first step to enable future simulation studies for practical applications. To achieve this goal, in-depth analyses of the physical, numerical, and modeling aspects related to Sc>>1 are presented, specifically when modeling in the viscous-convective subrange. Transport characteristics are scrutinized by examining scalar-velocity Fourier mode interactions in Direct Numerical Simulation (DNS) datasets and suggest that scalar modes in the viscous-convective subrange do not directly affect large-scale transport for high Sc. Further observations confirm that discretization errors inherent in numerical schemes can be sufficiently large to wipe out any meaningful contribution from subfilter models. This provides strong incentive to develop more effective numerical schemes to support high Sc simulations. To lower numerical dissipation while maintaining physically and mathematically appropriate scalar bounds during the convection step, a novel method of enforcing bounds is formulated, specifically for use with cubic Hermite polynomials. Boundedness of the scalar being transported is effected by applying derivative limiting techniques, and physically plausible single sub-cell extrema are allowed to exist to help minimize numerical dissipation. The proposed bounding algorithm results in significant performance gain in DNS of turbulent mixing layers and of homogeneous isotropic turbulence. Next, the combined physical/mathematical behavior of the subfilter scalar-flux vector is analyzed in homogeneous isotropic turbulence, by examining vector orientation in the strain-rate eigenframe. The results indicate no discernible dependence on the modeled scalar field, and lead to the identification of the tensor-diffusivity model as a good representation of the subfilter flux. Velocity resolved - scalar modeled simulations of homogeneous isotropic turbulence are conducted to confirm the behavior theorized in these a priori analyses, and suggest that the tensor-diffusivity model is ideal for use in the viscous-convective subrange. Simulations of a turbulent mixing layer are also discussed, with the partial objective of analyzing Schmidt number dependence of a variety of scalar statistics. Large-scale statistics are confirmed to be relatively independent of the Schmidt number for Sc>>1, which is explained by the dominance of subfilter dissipation over resolved molecular dissipation in the simulations. Overall, the VR-SM framework presented is quite effective in predicting large-scale transport characteristics of high Schmidt number scalars, however, it is determined that prediction of subfilter quantities would entail additional modeling intended specifically for this purpose. The VR-SM simulations presented in this thesis provide us with the opportunity to overlap with experimental studies, while at the same time creating an assortment of baseline datasets for future validation of LES models, thereby satisfying the objectives outlined for this work.

Kinetic theory of normal quantum fluids

Close to equilibrium, a normal Bose or Fermi fluid can be described by an exact kinetic equation whose kernel is nonlocal in space and time. The general expression derived for the kernel is evaluated to second order in the interparticle potential. The result is a wavevector- and frequency-dependent generalization of the linear Uehling-Uhlenbeck kernel with the Born approximation cross section.

The theory is formulated in terms of second-quantized phase space operators whose equilibrium averages are the n-particle Wigner distribution functions. Convenient expressions for the commutators and anticommutators of the phase space operators are obtained. The two-particle equilibrium distribution function is analyzed in terms of momentum-dependent quantum generalizations of the classical pair distribution function h(k) and direct correlation function c(k). The kinetic equation is presented as the equation of motion of a two -particle correlation function, the phase space density-density anticommutator, and is derived by a formal closure of the quantum BBGKY hierarchy. An alternative derivation using a projection operator is also given. It is shown that the method used for approximating the kernel by a second order expansion preserves all the sum rules to the same order, and that the second-order kernel satisfies the appropriate positivity and symmetry conditions.

Rapid changes in the vegetation of a shallow pond in Epping Forest, related to recent droughts

Over much of Britain, 1995 and 1996 have been perceived as drought years. To evaluate the impact that local climatic conditions are having upon successional changes in higher vegetation (macrophytes), Speakmans Pond in Epping Forest was surveyed and mapped in 1996. The results are related to previous vegetation surveys carried out in 1989 and 1991. In 1989 the dominant marginal vegetation was floating sweet-grass Glyceria fluitans, which also covered a major part of the main body of the pond. Other abundant species included soft rush Juncus effusus, reed mace Typha latifolia and yellow flag Iris pseudocorus. A small (central) area of open water contained bladderwort Utricularia vulgaris and white water-lily Nymphaea alba. A similar plant coverage was found in 1991, with a dominance of floating sweet-grass along the shallow eastern edge. A marked change in the pond was found during the 1996 survey of vegetation in July, when the pool was dry. The major plant cover now consisted of creeping bent Agrostis stolonifera, with isolated clumps of Yorkshire fog Holcus lanatus around the edges; both are terrestrial grasses found on land surrounding the pond. Rushes (Juncus) had increased their distribution round the margins of the pond, and the patch of yellow flag noted in 1989 and 1991 was not found in 1996. The deeper trenches were also dry, but a small patch of white water-lily remained adjacent to one of the trenches.

Environmental variables and fisheries diversity of the Naaf River Estuary, Bangladesh

The Naaf River estuary is one of the large estuaries in the Bangladesh coastal region not to have been affected by extensive human disturbance. This research provides information about the fisheries diversity status by Estuarine Set Bag Net (ESBN) sampling relation to physicochemical variables in both spatio-temporal scales. About 25 km of the lower estuary was divided into six zones for sample collection by considering the accessibility and availability of the ESBN operation, fish landing centers and location of the fishing villages. In total 48 samples have been analyzed which were taken throughout March to October 2006. To quantify the species diversity, all fisheries data were analyzed by using EstimateS and EcoSim software which accounts the different diversity indices viz., species richness, Shannon–Wiener diversity Index, Dominance and Evenness index. The research results demonstrate that the Naaf River estuary is a habitat of 161 (species richness, Sobs=161, Choa 1=162±2.34, ACE=161.73) different species which belong to 98 fin fishes, 23 shrimps and prawns,13 crabs, 11 molluscs, 3 echinoderms, 4 other crustaceans;while 9 remain unidentified. Results on the aquatic environment,mainly salinity and turbidity were found to have a major influence on their occurrence and distribution. All the findings indicated that the Naaf River estuary is a highly productive system and provides a favourable environment for large variety of estuarine species assemblages.

Valley evolution by meandering rivers

Fluvial systems form landscapes and sedimentary deposits with a rich hierarchy of structures that extend from grain- to valley scale. Large-scale pattern formation in fluvial systems is commonly attributed to forcing by external factors, including climate change, tectonic uplift, and sea-level change. Yet over geologic timescales, rivers may also develop large-scale erosional and depositional patterns that do not bear on environmental history. This dissertation uses a combination of numerical modeling and topographic analysis to identify and quantify patterns in river valleys that form as a consequence of river meandering alone, under constant external forcing. Chapter 2 identifies a numerical artifact in existing, grid-based models that represent the co-evolution of river channel migration and bank strength over geologic timescales. A new, vector-based technique for bank-material tracking is shown to improve predictions for the evolution of meander belts, floodplains, sedimentary deposits formed by aggrading channels, and bedrock river valleys, particularly when spatial contrasts in bank strength are strong. Chapters 3 and 4 apply this numerical technique to establishing valley topography formed by a vertically incising, meandering river subject to constant external forcing—which should serve as the null hypothesis for valley evolution. In Chapter 3, this scenario is shown to explain a variety of common bedrock river valley types and smaller-scale features within them—including entrenched channels, long-wavelength, arcuate scars in valley walls, and bedrock-cored river terraces. Chapter 4 describes the age and geometric statistics of river terraces formed by meandering with constant external forcing, and compares them to terraces in natural river valleys. The frequency of intrinsic terrace formation by meandering is shown to reflect a characteristic relief-generation timescale, and terrace length is identified as a key criterion for distinguishing these terraces from terraces formed by externally forced pulses of vertical incision. In a separate study, Chapter 5 utilizes image and topographic data from the Mars Reconnaissance Orbiter to quantitatively identify spatial structures in the polar layered deposits of Mars, and identifies sequences of beds, consistently 1-2 meters thick, that have accumulated hundreds of kilometers apart in the north polar layered deposits.

Composição e abundância das espécies de anfíbios e répteis em uma paisagem fragmentada de Mata Atlântica no Estado do Rio de Janeiro

Neste estudo, procuramos avaliar a composição e abundância de espécies de anfíbios e répteis em uma paisagem fragmentada de Mata Atlântica, no município de Cachoeiras de Macacu, Estado do Rio de Janeiro. Amostramos a herpetofauna da região na área contínua de floresta da Reserva Ecológica de Guapiaçu (REGUA), em 12 fragmentos do entorno com tamanhos e graus de isolamento diferentes, e áreas de pasto (matriz). Utilizamos para amostragem destes animais as metodologias de buscas ativas visuais e armadilhas de interceptação e queda, além de encontros ocasionais. No total foram registradas 55 espécies de anfíbios anuros pertencentes a 12 famílias e entre os répteis foram registradas 26 espécies, sendo uma espécie de anfisbena, uma de jacaré, nove de lagartos e 15 de serpentes. Para os anfíbios, houve uma dominância de espécies da família Hylidae, que representaram mais da metade do total de espécies encontradas no estudo. Já entre os répteis, houve uma predominância de espécies de serpentes da família Dipsadidae. Considerando apenas os registros feitos pelas metodologias empregadas, a área contínua de floresta da REGUA possui uma riqueza de espécies de anfíbios (N = 30) e de lagartos (N = 4) menor do que o conjunto de fragmentos (N = 36 e N = 8), mas superior ao que foi encontrado na matriz (N = 25 e N = 1). Entretanto, para os anfíbios, mais de um terço das espécies (N = 11) que ocorreu na mata contínua não ocorreu nos fragmentos ou na matriz, o que sugere que estas espécies podem ser mais sensíveis a alterações do hábitat. A maior riqueza de espécies encontrada no conjunto de fragmentos pode ser parcialmente explicada pelo fato de muitas espécies tanto de anfíbios quanto de lagartos que são típicas de áreas abertas, terem sua ocorrência favorecida neste tipo de condição de ambiente relativamente menos fechado dos fragmentos. Quando avaliamos o efeito de métricas da paisagem, observamos diferentes respostas entre os anfíbios e os lagartos. Enquanto para os anfíbios houve uma tendência de fragmentos mais distantes terem uma menor riqueza de espécies e de modos reprodutivos associados a estas espécies, para os lagartos a área dos fragmentos parece ser uma importante variável na estruturação das comunidades. Entretanto, por particularidades das características fisiológicas e ecológicas de anfíbios e lagartos, é possível que outros fatores expliquem a distribuição diferenciada das espécies. De forma geral, as áreas de matriz amostradas pareceram ser hostis a espécies florestais tanto de anfíbios quanto de lagartos. Além disso, para anfíbios, a presença de ambientes reprodutivos pode ser fator crucial para a ocorrência de algumas espécies. De forma similar ao encontrado em outros estudos, para a manutenção da diversidade de anfíbios e lagartos na paisagem fragmentada tratada neste estudo é necessário preservar o grande bloco florestal e aumentar a conectividade deste com os fragmentos, o que poderia permitir a área contínua servir de área-fonte de dispersores para os remanescentes florestais.

First principles based multiparadigm modeling of electronic structures and dynamics

Electronic structures and dynamics are the key to linking the material composition and structure to functionality and performance.

An essential issue in developing semiconductor devices for photovoltaics is to design materials with optimal band gaps and relative positioning of band levels. Approximate DFT methods have been justified to predict band gaps from KS/GKS eigenvalues, but the accuracy is decisively dependent on the choice of XC functionals. We show here for CuInSe₂ and CuGaSe₂, the parent compounds of the promising CIGS solar cells, conventional LDA and GGA obtain gaps of 0.0-0.01 and 0.02-0.24 eV (versus experimental values of 1.04 and 1.67 eV), while the historically first global hybrid functional, B3PW91, is surprisingly the best, with band gaps of 1.07 and 1.58 eV. Furthermore, we show that for 27 related binary and ternary semiconductors, B3PW91 predicts gaps with a MAD of only 0.09 eV, which is substantially better than all modern hybrid functionals, including B3LYP (MAD of 0.19 eV) and screened hybrid functional HSE06 (MAD of 0.18 eV).

The laboratory performance of CIGS solar cells (> 20% efficiency) makes them promising candidate photovoltaic devices. However, there remains little understanding of how defects at the CIGS/CdS interface affect the band offsets and interfacial energies, and hence the performance of manufactured devices. To determine these relationships, we use the B3PW91 hybrid functional of DFT with the AEP method that we validate to provide very accurate descriptions of both band gaps and band offsets. This confirms the weak dependence of band offsets on surface orientation observed experimentally. We predict that the CBO of perfect CuInSe₂/CdS interface is large, 0.79 eV, which would dramatically degrade performance. Moreover we show that band gap widening induced by Ga adjusts only the VBO, and we find that Cd impurities do not significantly affect the CBO. Thus we show that Cu vacancies at the interface play the key role in enabling the tunability of CBO. We predict that Na further improves the CBO through electrostatically elevating the valence levels to decrease the CBO, explaining the observed essential role of Na for high performance. Moreover we find that K leads to a dramatic decrease in the CBO to 0.05 eV, much better than Na. We suggest that the efficiency of CIGS devices might be improved substantially by tuning the ratio of Na to K, with the improved phase stability of Na balancing phase instability from K. All these defects reduce interfacial stability slightly, but not significantly.

A number of exotic structures have been formed through high pressure chemistry, but applications have been hindered by difficulties in recovering the high pressure phase to ambient conditions (i.e., one atmosphere and room temperature). Here we use dispersion-corrected DFT (PBE-ulg flavor) to predict that above 60 GPa the most stable form of N₂O (the laughing gas in its molecular form) is a 1D polymer with an all-nitrogen backbone analogous to cis-polyacetylene in which alternate N are bonded (ionic covalent) to O. The analogous trans-polymer is only 0.03-0.10 eV/molecular unit less stable. Upon relaxation to ambient conditions both polymers relax below 14 GPa to the same stable non-planar trans-polymer, accompanied by possible electronic structure transitions. The predicted phonon spectrum and dissociation kinetics validate the stability of this trans-poly-NNO at ambient conditions, which has potential applications as a new type of conducting polymer with all-nitrogen chains and as a high-energy oxidizer for rocket propulsion. This work illustrates in silico materials discovery particularly in the realm of extreme conditions.

Modeling non-adiabatic electron dynamics has been a long-standing challenge for computational chemistry and materials science, and the eFF method presents a cost-efficient alternative. However, due to the deficiency of FSG representation, eFF is limited to low-Z elements with electrons of predominant s-character. To overcome this, we introduce a formal set of ECP extensions that enable accurate description of p-block elements. The extensions consist of a model representing the core electrons with the nucleus as a single pseudo particle represented by FSG, interacting with valence electrons through ECPs. We demonstrate and validate the ECP extensions for complex bonding structures, geometries, and energetics of systems with p-block character (C, O, Al, Si) and apply them to study materials under extreme mechanical loading conditions.

Despite its success, the eFF framework has some limitations, originated from both the design of Pauli potentials and the FSG representation. To overcome these, we develop a new framework of two-level hierarchy that is a more rigorous and accurate successor to the eFF method. The fundamental level, GHA-QM, is based on a new set of Pauli potentials that renders exact QM level of accuracy for any FSG represented electron systems. To achieve this, we start with using exactly derived energy expressions for the same spin electron pair, and fitting a simple functional form, inspired by DFT, against open singlet electron pair curves (H₂ systems). Symmetric and asymmetric scaling factors are then introduced at this level to recover the QM total energies of multiple electron pair systems from the sum of local interactions. To complement the imperfect FSG representation, the AMPERE extension is implemented, and aims at embedding the interactions associated with both the cusp condition and explicit nodal structures. The whole GHA-QM+AMPERE framework is tested on H element, and the preliminary results are promising.

Storm track response to perturbations in climate

This thesis advances our understanding of midlatitude storm tracks and how they respond to perturbations in the climate system. The midlatitude storm tracks are regions of maximal turbulent kinetic energy in the atmosphere. Through them, the bulk of the atmospheric transport of energy, water vapor, and angular momentum occurs in midlatitudes. Therefore, they are important regulators of climate, controlling basic features such as the distribution of surface temperatures, precipitation, and winds in midlatitudes. Storm tracks are robustly projected to shift poleward in global-warming simulations with current climate models. Yet the reasons for this shift have remained unclear. Here we show that this shift occurs even in extremely idealized (but still three-dimensional) simulations of dry atmospheres. We use these simulations to develop an understanding of the processes responsible for the shift and develop a conceptual model that accounts for it.

We demonstrate that changes in the convective static stability in the deep tropics alone can drive remote shifts in the midlatitude storm tracks. Through simulations with a dry idealized general circulation model (GCM), midlatitude storm tracks are shown to be located where the mean available potential energy (MAPE, a measure of the potential energy available to be converted into kinetic energy) is maximal. As the climate varies, even if only driven by tropical static stability changes, the MAPE maximum shifts primarily because of shifts of the maximum of near-surface meridional temperature gradients. The temperature gradients shift in response to changes in the width of the tropical Hadley circulation, whose width is affected by the tropical static stability. Storm tracks generally shift in tandem with shifts of the subtropical terminus of the Hadley circulation.

We develop a one-dimensional diffusive energy-balance model that links changes in the Hadley circulation to midlatitude temperature gradients and so to the storm tracks. It is the first conceptual model to incorporate a dynamical coupling between the tropical Hadley circulation and midlatitude turbulent energy transport. Numerical and analytical solutions of the model elucidate the circumstances of when and how the storm tracks shift in tandem with the terminus of the Hadley circulation. They illustrate how an increase of only the convective static stability in the deep tropics can lead to an expansion of the Hadley circulation and a poleward shift of storm tracks.

The simulations with the idealized GCM and the conceptual energy-balance model demonstrate a clear link between Hadley circulation dynamics and midlatitude storm track position. With the help of the hierarchy of models presented in this thesis, we obtain a closed theory of storm track shifts in dry climates. The relevance of this theory for more realistic moist climates is discussed.

A luta das enfermeiras por um espaço na FAB: a turma pioneira de oficiais (1981-1984)

Este estudo, de cunho histórico-social, tem como objeto a inserção de enfermeiras como oficiais da Força Aérea Brasileira (FAB) por meio do pioneiro Quadro Feminino de Oficiais (QFO). O marco inicial do estudo refere-se ao início do Estágio de Adaptação militar, em 02 de agosto de 1982 no Centro de Instrução Especializada da Aeronáutica (CIEAR), localizado na cidade do Rio de Janeiro. O marco final do estudo diz respeito ao término do período inicial obrigatório de dois anos de cumprimento de serviço ativo dessas enfermeiras, que culminou com a promoção das mesmas ao posto de 1Tenente (1984). Os objetivos do estudo são: descrever as circunstâncias de inserção das enfermeiras no processo seletivo do QFO, analisar o processo de incorporação do habitus militar durante o Estágio de Adaptação, e discutir as estratégias de luta das enfermeiras militares para ocuparem seus lugares devidos nos hospitais da FAB. A técnica de coleta de dados utilizada foi a entrevista e ocorreu no período de abril a maio de 2009 em hospitais da FAB da cidade do Rio de Janeiro. Foram entrevistadas cinco enfermeiras militares da primeira turma do QFO. O estudo foi cadastrado no SISNEP e aprovado pelo Comitê de Ética da FAB. Todos os sujeitos assinaram o Termo de consentimento livre e esclarecido e o Termo de doação de depoimento oral. O método utilizado foi o da História oral temática o referencial teórico do estudo foi baseado no pensamento do sociólogo francês Pierre Bourdieu, cujos conceitos de poder simbólico, habitus, campo, espaço social e violência simbólica sustentaram a construção desta dissertação. Para a análise e interpretação dos dados, seguimos os passos propostos por Maria Cecília Minayo de ordenação de dados, que compreendeu a transcrição na íntegra dos depoimentos; classificação cronológica e temática dos documentos escritos; classificação dos dados e a análise final. Evidenciou-se que diversos motivos incentivaram as enfermeiras a almejarem sua inserção na FAB como a boa remuneração, estabilidade financeira, progressão profissional, desbravamento de um novo campo de trabalho, clientela diferenciada, aposentadoria com salário integral e pioneirismo na FAB. O objetivo do Estágio de Adaptação militar foi inculcar do habitus militar nas candidatas a partir de ensinamentos baseados na hierarquia, disciplina, ética, dever e compromisso militar. Ao se inserirem nos hospitais da FAB, as enfermeiras receberam diversos cargos e funções, galgando um poder simbólico sobre a equipe de enfermagem. As inevitáveis lutas simbólicas dessas enfermeiras ocorreram com os médicos militares, com a equipe de enfermagem, com as enfermeiras civis e com a própria administração do hospital, e revelaram aspectos característicos de violência simbólica desencadeada por lutas de gênero e pela manutenção do poder, visto que as enfermeiras, dotadas de status de chefe e de militar, se inseriram num campo eminentemente masculino.

PP backward elastic scattering between 0.70 and 2.37 GeV/c

̄pp backward elastic scattering has been measured for the cos θ_cm region between – 1.00 and – 0.88 and for the incident ̄p laboratory momentum region between 0.70 and 2.37 GeV/c. These measurements, done in intervals of approximately 0.1 GeV/c, have been performed at the Alternating Gradient Synchrotron at Brookhaven National Laboratory during the winter of 1968. The measured differential cross sections, binned in cos θ_cm intervals of 0.02, have statistical errors of about 10%. Backward dipping exists below 0.95 GeV/c and backward peaking above 0.95 GeV/c. The 180˚ differential cross section extrapolated from our data shows a sharp dip centered at 0.95 GeV/c and a broad hump centered near 1.4 GeV/c. Our data have been interpreted in terms of resonance effects and in terms of diffraction dominance effects.