Spectral evolution of angularly resolved third-order harmonic generation by infrared femtosecond laser-pulse filamentation in air

We experimentally investigate the evolution of an angularly resolved spectrum of third harmonic generated by infrared femtosecond laser pulse filamentation in air. We show that at low pump intensity, phase matching between the fundamental and third-harmonic waves dominates the nonlinear optical effect and induces a ring structure of the third-harmonic beam, whereas at high pump intensity, the dispersion properties of air begin to affect the angular spectrum, leading to the formation of a nonlinear X wave at third harmonic.

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.

Characterization and modeling of premixed turbulent n-heptane flames in the thin reaction zone regime

n-heptane/air premixed turbulent flames in the high-Karlovitz portion of the thin reaction zone regime are characterized and modeled in this thesis using Direct Numerical Simulations (DNS) with detailed chemistry. In order to perform these simulations, a time-integration scheme that can efficiently handle the stiffness of the equations solved is developed first. A first simulation with unity Lewis number is considered in order to assess the effect of turbulence on the flame in the absence of differential diffusion. A second simulation with non-unity Lewis numbers is considered to study how turbulence affects differential diffusion. In the absence of differential diffusion, minimal departure from the 1D unstretched flame structure (species vs. temperature profiles) is observed. In the non-unity Lewis number case, the flame structure lies between that of 1D unstretched flames with "laminar" non-unity Lewis numbers and unity Lewis number. This is attributed to effective Lewis numbers resulting from intense turbulent mixing and a first model is proposed. The reaction zone is shown to be thin for both flames, yet large chemical source term fluctuations are observed. The fuel consumption rate is found to be only weakly correlated with stretch, although local extinctions in the non-unity Lewis number case are well correlated with high curvature. These results explain the apparent turbulent flame speeds. Other variables that better correlate with this fuel burning rate are identified through a coordinate transformation. It is shown that the unity Lewis number turbulent flames can be accurately described by a set of 1D (in progress variable space) flamelet equations parameterized by the dissipation rate of the progress variable. In the non-unity Lewis number flames, the flamelet equations suggest a dependence on a second parameter, the diffusion of the progress variable. A new tabulation approach is proposed for the simulation of such flames with these dimensionally-reduced manifolds.

Emissões de compostos orgânicos voláteis de um aterro controlado e o potencial formador de ozônio

O presente estudo teve como objetivo quantificar as emissões de Compostos Orgânicos Voláteis do Aterro Controlado Morro do Céu localizado na cidade de Niterói, Rio de Janeiro, Brasil. Para tanto, vinte amostras foram coletadas, usando uma bomba de ar operada a bateria durante dois dias de dezembro de 2009. Uma câmara de fluxo cilíndrica de PVC de 30L foi inserida 5 cm no solo do aterro, e as amostras foram coletadas através de uma válvula na parte superior da câmera. Os resultados indicaram um valor de 1.980 Kg Km-2 h-1. O modelo Gaussiano de dispersão atmosférica ISCST3 foi utilizado para calcular a difusão e transporte dos poluentes a fim de estimar as concentrações de COV no bairro, usando dados topográficos, meteorológicos e de emissões. Valores de 525 μg m-3 de COV foram encontrados a 500 metros do aterro. As emissões do aterro foram usadas em conjunto com dados meteorológicos, utilizando o modelo de trajetória OZIPR e o mecanismo químico SAPRC para demonstrar o impacto na formação do ozônio troposférico na região. É conhecido que o ozônio é formado pela reação entre COV, NOx e luz solar. A contribuição de valores elevados de COV provenientes das emissões do aterro conduzirá a uma nova situação com valores mais elevados de ozônio na região. Os resultados da modelagem indicaram um aumento maior que 1000% nos níveis de ozônio na região do aterro, se comparado com a modelagem do ozônio para a região metropolitana do Rio de Janeiro. Os resultados mostram ser necessário que maior atenção seja dada à política de gerenciamento de RSU no Brasil, incluindo a escolha adequada para o local de instalação, o monitoramento da área durante e após o período de operação e técnicas mais adequadas de disposição dos resíduos sólidos urbanos

Emissões atmosféricas e implicações potenciais sobre a biota terrestre devido as atividades antrópicas na Baía do Almirantado/ Ilha Rei George Antártica

Neste trabalho foi realizado um estudo sobre o impacto sobre a atmosfera e a biota terrestre devido às emissões antrópicas na Baia do Almirantado/Ilha Rei George Antártica. Foram monitoradas as emissões dos compostos orgânicos voláteis e semi-voláteis, tanto nas fontes emissoras como no entorno da Estação Antártica Comandante Ferraz e estimadas as emissões dos navios, dos geradores a diesel e da incineração de lixo. Na avaliação do entorno, coletaram-se amostras de ar, neve e penas de aves. Com os resultados das emissões, do estudo topográfico e da meteorologia, realizou-se uma modelagem de plumas gaussiana para avaliar os impactos. Quatro cenários foram avaliados: dois com a presença dos navios NApOc Ary Rongel e Maximiano apresentaram concentrações máximas de até 356 g m-3 de COV e 18 g m-3 de material particulado, enquanto os demais, sem a presença dos navios, apenas considerando as estações de pesquisa EACF e Arctowski, apresentaram concentrações máximas de até 2,5 g m-3 de COV e 1,3 g m-3 de material particulado. Amostras de COV coletadas foram compatíveis com o cenário mais crítico. O estudo de correlação para carbonilas e HPA atmosférico e Carbono elementar e HPA, depositados em neve, apontaram a EACF como a principal fonte de emissão. As concentrações de levoglucosano detectadas a aproximadamente 2 km da EACF apontaram para a prática de incineração de lixo da EACF. Todas as áreas de interesse biológicos, anteriormente mapeados, dentro da AAEG, são vulneráveis às emissões antrópicas, como sugeriu o modelo de dispersão e a sobreposição dos resultados encontrados.

Avaliação do impacto radiológico atmosférico de uma unidade de mineração e beneficiamento de urânio

Este trabalho teve como objetivo avaliar o impacto radiológico atmosférico da Unidade de Concentrado de Urânio URA, Caetité, BA, através da modelagem da dispersão de radionuclídeos e a estimativa da dose efetiva anual (em mSv.ano-1). Para tal, utilizou-se o programa MILDOS-AREA que foi desenvolvido pelo Argonne National Laboratory (ANL) em conjunto com a U.S. Nuclear Regulatory Commission (USNRC), para avaliar impacto radiológico ambiental atmosférico nas instalações de mineração e beneficiamento de urânio. O incremento de dose efetiva anual para três grupos críticos hipotéticos e oito grupos populacionais reais foi estimado com base na medida de fluxos de radônio e na estimativa das concentrações de radionuclídeos em particulados no ar dos principais termos fontes da URA (cava da mina, depósito de estéril e britador). Paralelamente, as medidas de concentração de radônio e taxa de kerma no ar, reportadas nos relatórios dos programas de monitoração ambiental pré-operacional (PMAPO) e operacional (PMAO) da URA, foram avaliadas. Os valores de dose efetiva anual estimados para os grupos críticos hipotéticos variaram de 1,78E-02 a 2,10E-02 mSv.ano-1, enquanto que para os grupos populacionais, variaram de 7,49E-05 a 1,56E-02 mSv.ano-1. A maior contribuição para o incremento da dose foi devida a inalação do radônio, sendo responsável por quase a totalidade da dose efetiva anual estimada. A média da concentração de atividade de radônio no entorno da URA foi 137,21 Bq m-3 e não sendo observada diferenças significativas entre as concentrações de radônio reportadas nos programas de monitoramento ambiental pré-operacional (valores de background) e operacional. Os valores médios de taxa de kerma no ar no entorno da URA foram de 0,136 μGy h-1. No entanto, em todos os pontos de monitoramento, os valores reportados no programa operacional foram inferiores aos valores reportados no programa pré-operacional (background), o que sugere problemas de medidas ou de coleta de dados durante a realização deste programa. O operador da URA utilizou para avaliação de impacto radiológico atmosférico, resultados apresentados em seus relatórios finais de análise de segurança (RFAS), um modelo próprio de simulação de dispersão, denominado Impacto Ambiental Radiológico (IAR7). Uma comparação entre o MILDOS-AREA e o IAR7, utilizando os mesmos parâmetros de entrada reportados no RFAS sugere que o IAR7 subestimou as concentrações de radônio no ar para os grupos críticos hipotéticos. Os resultados de simulação com o MILDOS-AREA mostram que as doses efetivas estimadas para os grupos críticos hipotéticos são inferiores a 0,3 mSv.ano-1 que é a restrição de dose estabelecida pela Comissão Nacional de Energia Nuclear. Recomenda-se que o código MILDOS-AREA seja utilizado no Brasil, para fins de licenciamento e controle, tendo em vista que o mesmo é um código validado e já utilizado em outros países para avaliar impacto radiológico ambiental atmosférico em instalações de mineração e beneficiamento de urânio

Modeling evaporation effects in conditional moment closure for spray autoignition

Simulations of an n-heptane spray autoigniting under conditions relevant to a diesel engine are performed using two-dimensional, first-order conditional moment closure (CMC) with full treatment of spray terms in the mixture fraction variance and CMC equations. The conditional evaporation term in the CMC equations is closed assuming interphase exchange to occur at the droplet saturation mixture fraction values only. Modeling of the unclosed terms in themixture fraction variance equation is done accordingly. Comparison with experimental data for a range of ambient oxygen concentrations shows that the ignition delay is overpredicted. The trend of increasing ignition delay with decreasing oxygen concentration, however, is correctly captured. Good agreement is found between the computed and measured flame lift-off height for all conditions investigated. Analysis of source terms in the CMC temperature equation reveals that a convective-reactive balance sets in at the flame base, with spatial diffusion terms being important, but not as important as in lifted jet flames in cold air. Inclusion of droplet terms in the governing equations is found to affect the mixture fraction variance field in the region where evaporation is the strongest, and to slightly increase the ignition delay time due to the cooling associated with the evaporation. Both flame propagation and stabilization mechanisms, however, remain unaffected. © 2011 Taylor & Francis.

Baroclinic waves in an air-filled thermally driven rotating annulus.

In this study an experimental investigation of baroclinic waves in air in a differentially heated rotating annulus is presented. Air has a Prandtl number of 0.707, which falls within a previously unexplored region of parameter space for baroclinic instability. The flow regimes encountered include steady waves, periodic amplitude vacillations, modulated amplitude vacillations, and either monochromatic or mixed wave number weak waves, the latter being characterized by having amplitudes less than 5% of the applied temperature contrast. The distribution of these flow regimes in parameter space are presented in a regime diagram. It was found that the progression of transitions between different regimes is, as predicted by recent numerical modeling results, in the opposite sense to that usually found in experiments with high Prandtl number liquids. No hysteresis in the flow type, with respect to variations in the rotation rate, was found in this investigation.