Analytical optimization of spread and change of exponential decay of generalized Wannier functions in one dimension

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

Transient non-Darcy forced convective heat transfer from a flat plate embedded in a fluid-saturated porous medium

Transient non-Darcy forced convection on a flat plate embedded in a porous medium is investigated using the Forchheimer-extended Darcy law. A sudden uniform pressure gradient is applied along the flat plate, and at the same time, its wall temperature is suddenly raised to a high temperature. Both the momentum and energy equations are solved by retaining the unsteady terms. An exact velocity solution is obtained and substituted into the energy equation, which then is solved by means of a quasi-similarity transformation. The temperature field can be divided into the one-dimensional transient (downstream) region and the quasi-steady-state (upstream) region. Thus the transient local heat transfer coefficient can be described by connecting the quasi-steady-state solution and the one-dimensional transient solution. The non-Darcy porous inertia works to decrease the velocity level and the time required for reaching the steady-state velocity level. The porous-medium inertia delays covering of the plate by the steady-state thermal boundary layer. © 1990.

An improved three-dimensional model for growth of oxide films induced by laser heating

In this paper we consider a three-dimensional heat diffusion model to explain the growth of oxide films which takes place when a laser beam is shined on and heats a metallic layer deposited on a glass substrate in a normal atmospheric environment. In particular, we apply this model to the experimental results obtained for the dependence of the oxide layer thickness on the laser density power for growth of TiO2 films grown on Ti-covered glass slides. We show that there is a very good agreement between the experimental results and the theoretical predictions from our proposed three-dimensional model, improving the results obtained with the one-dimensional heat diffusion model previously reported. Our theoretical results also show the occurrence of surface cooling between consecutive laser pulses, and that the oxide track surface profile closely follows the spatial laser profile indicating that heat diffusive effects can be neglected in the growth of oxide films by laser heating. © 2001 Elsevier Science B.V.

Localization of a Bose-Einstein condensate in a bichromatic optical lattice

By direct numerical simulation of the time-dependent Gross-Pitaevskii equation, we study different aspects of the localization of a noninteracting ideal Bose-Einstein condensate (BEC) in a one-dimensional bichromatic quasiperiodic optical-lattice potential. Such a quasiperiodic potential, used in a recent experiment on the localization of a BEC, can be formed by the superposition of two standing-wave polarized laser beams with different wavelengths. We investigate the effect of the variation of optical amplitudes and wavelengths on the localization of a noninteracting BEC. We also simulate the nonlinear dynamics when a harmonically trapped BEC is suddenly released into a quasiperiodic potential, as done experimentally in a laser speckle potential. We finally study the destruction of the localization in an interacting BEC due to the repulsion generated by a positive scattering length between the bosonic atoms. © 2009 The American Physical Society.

Modelagem numérica da influência do eletrojato equatorial em dados magnetotelúricos produzidos por estruturas tridimensionais

A América do Sul apresenta várias peculiaridades geomagnéticas, uma delas, é a presença do Eletrojato Equatorial, o qual se estende de leste para oeste no Brasil ao longo de aproximadamente 3500 km. Considerando-se o fato de que a influência do Eletrojato Equatorial pode ser detectada a grandes distâncias do seu centro, isto suscita o interesse em se estudar os seus efeitos na exploração magnetotelúrica no Brasil. A influência do eletrojato equatorial na prospecção magnetotelúrica tem sido modelada para meios geológicos uni e bidimensionais valendo-se para isto de soluções analíticas fechadas e de técnicas numéricas tais como elementos finitos e diferenças finitas. Em relação aos meios geológicos tridimensionais, eles tem sido modelados na forma de "camadas finas", usando o algoritmo "thin sheet". As fontes indutoras utilizadas para simular o eletrojato equatorial nestes trabalhos, tem sido linhas de corrente, eletrojatos gaussianos e eletrojatos ondulantes. Por outro lado, o objetivo principal da nossa tese foi o modelamento dos efeitos que o eletrojato equatorial provoca em estruturas tridimensionais próprias da geofísica da prospecção. Com tal finalidade, utilizamos o esquema numérico da equação integral, com as fontes indutoras antes mencionadas. De maneira similar aos trabalhos anteriores, os nossos resultados mostram que a influência do eletrojato equatorial somente acontece em frequências menores que 10^-1 Hz. Este efeito decresce com a distância, mantendo-se até uns 3000 km do centro do eletrojato. Assim sendo, a presença de grandes picos nos perfis da resistividade aparente de um semi-espaço homogêneo, indica que a influência do eletrojato é notável neste tipo de meio. Estes picos se mostram com diferente magnitude para cada eletrojato simulado, sendo que a sua localização também muda de um eletrojato para outro. Entretanto, quando se utilizam modelos geo-elétricos unidimensionais mais de acordo com a realidade, tais como os meios estratificados, percebe-se que a resposta dos eletrojatos se amortece significativamente e não mostra muitas diferenças entre os diferentes tipos de eletrojato. Isto acontece por causa da dissipação da energia eletromagnética devido à presença da estratificação e de camadas condutivas. Dentro do intervalo de 3000 km, a resposta eletromagnética tridimensional pode ser deslocada para cima ou para baixo da resposta da onda plana, dependendo da localização do corpo, da frequência, do tipo de eletrojato e do meio geológico. Quando a resposta aparece deslocada para cima, existe um afastamento entre as sondagens uni e tridimensionais devidas ao eletrojato, assim como um alargamento da anomalia dos perfis que registra a presença da heterogeneidade tridimensional. Quando a resposta aparece deslocada para baixo, no entanto, há uma aproximação entre estes dois tipos de sondagens e um estreitamento da anomalia dos perfis. Por outro lado, a fase se mostra geralmente, de uma forma invertida em relação à resistividade aparente. Isto significa que quando uma sobe a outra desce, e vice-versa. Da mesma forma, comumente nas altas frequências as respostas uni e tridimensionais aparecem deslocadas, enquanto que nas baixas frequências se mostram com os mesmos valores, com exceção dos eletrojatos ondulantes com parâmetros de ondulação α = —2 e —3. Nossos resultados também mostram que características geométricas próprias das estruturas tridimensionais, tais como sua orientação em relação à direção do eletrojato e a dimensão da sua direção principal, afetam a resposta devido ao eletrojato em comparação com os resultados da onda plana. Desta forma, quando a estrutura tridimensional é rotacionada de 90°, em relação à direção do eletrojato e em torno do eixo z, existe uma troca de polarizações nas resistividades dos resultados, mas não existem mudanças nos valores da resistividade aparente no centro da estrutura. Ao redor da mesma, porém, se percebe facilmente alterações nos contornos dos mapas de resistividade aparente, ao serem comparadas com os mapas da estrutura na sua posição original. Isto se deve à persistência dos efeitos galvânicos no centro da estrutura e à presença de efeitos indutivos ao redor do corpo tridimensional. Ao alongar a direção principal da estrutura tridimensional, as sondagens magnetotelúricas vão se aproximando das sondagens das estruturas bidimensionais, principalmente na polarização XY. Mesmo assim, as respostas dos modelos testados estão muito longe de se considerar próximas das respostas de estruturas quase-bidimensionais. Porém, os efeitos do eletrojato em estruturas com direção principal alongada, são muito parecidos com aqueles presentes nas estruturas menores, considerando-se as diferenças entre as sondagens de ambos tipos de estruturas. Por outro lado, os mapas de resistividade aparente deste tipo de estrutura alongada, revelam um grande aumento nos extremos da estrutura, tanto para a onda plana como para o eletrojato. Este efeito é causado pelo acanalamento das correntes ao longo da direção principal da estrutura. O modelamento de estruturas geológicas da Bacia de Marajó confirma que os efeitos do eletrojato podem ser detetados em estruturas pequenas do tipo "horst" ou "graben", a grandes distâncias do centro do mesmo. Assim, os efeitos do eletrojato podem ser percebidos tanto nos meios estratificados como tridimensionais, em duas faixas de freqüência (nas proximidades de 10^-1 Hz e para freqüências menores que 10^-3 Hz), possivelmente influenciados pela presença do embasamento cristalino e a crosta inferior, respectivamente. Desta maneira, os resultados utilizando o eletrojato como fonte indutora, mostram que nas baixas freqüências as sondagens magnetotelúricas podem ser fortemente distorcidas, tanto pelos efeitos galvânicos da estrutura tridimensional como pela presença da influência do eletrojato. Conseqüêntemente, interpretações errôneas dos dados de campo podem ser cometidas, se não se corrigirem os efeitos do eletrojato equatorial ou, da mesma forma, não se utilisarem algoritmos tridimensionais para interpretar os dados, no lugar do usual modelo unidimensional de Tikhonov - Cagniard.

Funções de Wannier generalizadas para aplicações em Ciência dos Materiais

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

A new approach to the electronic structure of two dimensional disordered alloys

We propose a novel method to calculate the electronic Density of States (DOS) of a two dimensional disordered binary alloy. The method is highly reliable and numerically efficient, and Short Range Order (SRO) correlations can be included with no extra computational cost. The approach devised rests on one dimensional calculations and is applied to very long stripes of finite width, the bulk regime being achieved with a relatively small number of chains in the disordered case. Our approach is exact for the pure case and predicts the correct DOS structure in important limits, such as the segregated, random, and ordered alloy regimes. We also suggest important extensions of the present work. © 1995.

Is the meander growth in the Brazil Current system off Southeast Brazil due to baroclinic instability?

Temporally-growing frontal meandering and occasional eddy-shedding is observed in the Brazil Current (BC) as it flows adjacent to the Brazilian Coast. No study of the dynamics of this phenomenon has been conducted to date in the region between 22 degrees S and 25 degrees S. Within this latitude range, the flow over the intermediate continental slope is marked by a current inversion at a depth that is associated with the Intermediate Western Boundary Current (IWBC). A time series analysis of 10-current-meter mooring data was used to describe a mean vertical profile for the BC-IWBC jet and a typical meander vertical structure. The latter was obtained by an empirical orthogonal function (EOF) analysis that showed a single mode explaining 82% of the total variance. This mode structure decayed sharply with depth, revealing that the meandering is much more vigorous within the BC domain than it is in the IWBC region. As the spectral analysis of the mode amplitude time series revealed no significant periods, we searched for dominant wavelengths. This search was done via a spatial EOF analysis on 51 thermal front patterns derived from digitized AVHRR images. Four modes were statistically significant at the 95% confidence level. Modes 3 and 4, which together explained 18% of the total variance, are associated with 266 and 338-km vorticity waves, respectively. With this new information derived from the data, the [Johns, W.E., 1988. One-dimensional baroclinically unstable waves on the Gulf Stream potential vorticity gradient near Cape Hatteras. Dyn. Atmos. Oceans 11, 323-350] one-dimensional quasi-geostrophic model was applied to the interpolated mean BC-IWBC jet. The results indicated that the BC system is indeed baroclinically unstable and that the wavelengths depicted in the thermal front analysis are associated with the most unstable waves produced by the model. Growth rates were about 0.06 (0.05) days(-1) for the 266-km (338-km) wave. Moreover, phase speeds for these waves were low compared to the surface BC velocity and may account for remarks in the literature about growing standing or stationary meanders off southeast Brazil. The theoretical vertical structure modes associated with these waves resembled very closely to the one obtained for the current-meter mooring EOF analysis. We interpret this agreement as a confirmation that baroclinic instability is an important mechanism in meander growth in the BC system. (C) 2008 Elsevier B.V. All rights reserved.

Magnetic Tuning of All-Organic Binary Alloys between Two Stable Radicals

Mixtures of 2-(4,5,6,7-tetrafluorobenzimidazol-2-yl)-4,4,5,5-tetramethyl-4,5-dihydro-1H-imidazole-3-oxide-1-oxyl (F4BImNN) and 2-(benzi-midazol-2-yl)-4,4,5,5-tetramethyl-4,5-dihydro-1H-imidazole-3-oxide-1-oxyl (BImNN.) crystallize as solid solutions (alloys) across a wide range of binary compositions. (F4BImNN)(x)(BImNN)((1-x)) with x < 0.8 gives orthorhombic unit cells, while x >= 0.9 gives monoclinic unit cells. In all crystalline samples, the dominant intermolecular packing is controlled by one-dimensional (1D) hydrogen-bonded chains that lead to quasi-1D ferromagnetic behavior. Magnetic analysis over 0.4-300 K indicates ordering with strong 1D ferromagnetic exchange along the chains (J/k = 12-22 K). Interchain exchange is estimated to be 33- to 150-fold weaker, based on antiferromagnetic ordered phase formation below Neel temperatures in the 0.4-1.2 K range for the various compositions. The ordering temperatures of the orthorhombic samples increase linearly as (1 - x) increases from 0.25 to 1.00. The variation is attributed to increased interchain distance corresponding to decreased interchain exchange, when more F4BImNN is added into the orthorhombic lattice. The monoclinic samples are not part of the same trend, due to the different interchain arrangement associated with the phase change.

Präparation von Bismutsulfid-Absorberschichten für Dünnschichtsolarzellen mittels physikalischer Gasphasenabscheidung

Angesichts der sich abzeichnenden Erschöpfung fossiler Ressourcen ist die Erforschung alternativer Energiequellen derzeit eines der meistbeachteten Forschungsgebiete. Durch ihr enormes Potential ist die Photovoltaik besonders im Fokus der Wissenschaft. Um großflächige Beschichtungsverfahren nutzen zu können, wird seit einigen Jahren auf dem Gebiet der Dünnschichtphotovoltaik intensiv geforscht. Jedoch sind die gegenwärtigen Solarzellenkonzepte allesamt durch die Verwendung giftiger (Cd, As) oder seltener Elemente (In, Ga) oder durch eine komplexe Phasenbildung in ihrem Potential beschränkt. Die Entwicklung alternativer Konzepte erscheint daher naheliegend.rnAufgrund dessen wurde in einem BMBF-geförderten Verbundprojekt die Abscheidung von Dünnschichten des binären Halbleiters Bi2S3 mittels physikalischer Gasphasenabscheidung mit dem Ziel der Etablierung als quasi-intrinsischer Absorber in Solarzellenstrukturen mit p-i-n-Schichtfolge hin untersucht.rnDurch sein von einem hochgradig anisotropen Bindungscharakter geprägtes Kristallwachstum war die Abscheidung glatter, einphasiger und für die Integration in eine Multischichtstruktur geeigneter Schichten mit Schichtdicken von einigen 100 nm eine der wichtigsten Herausforderungen. Die Auswirkungen der beiden Parameter Abscheidungstemperatur und Stöchiometrie wurden hinsichtlich ihrer Auswirkungen auf die relevanten Kenngrößen (wie Morphologie, Dotierungsdichte und Photolumineszenz) untersucht. Es gelang, erfolgreich polykristalline Schichten mit geeigneter Rauigkeit und einer Dotierungsdichte von n ≈ 2 1015cm-3 auf anwendungsrelevanten Substraten abzuscheiden, wobei eine besonders starke Abhängigkeit von der Gasphasenzusammensetzung ermittelt werden. Es konnten weiterhin die ersten Messungen der elektronischen Zustandsdichte unter Verwendung von Hochenergie-Photoemissionsspektroskopie durchgeführt werden, die insbesondere den Einfluss variabler Materialzusammensetzungen offenbarten.rnZum Nachweis der Eignung des Materials als Absorberschicht standen innerhalb des Projektes mit SnS, Cu2O und PbS prinzipiell geeignete p-Kontaktmaterialien zur Verfügung. Es konnten trotz der Verwendung besonders sauberer Abscheidungsmethoden im Vakuum keine funktionstüchtigen Solarzellen mit Bi2S3 deponiert werden. Jedoch war es unter Verwendung von Photoemissionspektroskopie möglich, die relevanten Grenzflächen zu spektroskopieren und die Ursachen für die Beobachtungen zu identifizieren. Zudem konnte erfolgreich die Notwendigkeit von Puffermaterialien bei der Bi2S3-Abscheidung nachgewiesen werden, um Oberflächenreaktionen zu unterbinden und die Transporteigenschaften an der Grenzfläche zu verbessern.rn

STOCHASTIC CHARGE TRANSPORT IN MULTI-ISLAND SINGLE-ELECTRON TUNNELING DEVICES

The physics of the operation of singe-electron tunneling devices (SEDs) and singe-electron tunneling transistors (SETs), especially of those with multiple nanometer-sized islands, has remained poorly understood in spite of some intensive experimental and theoretical research. This computational study examines the current-voltage (IV) characteristics of multi-island single-electron devices using a newly developed multi-island transport simulator (MITS) that is based on semi-classical tunneling theory and kinetic Monte Carlo simulation. The dependence of device characteristics on physical device parameters is explored, and the physical mechanisms that lead to the Coulomb blockade (CB) and Coulomb staircase (CS) characteristics are proposed. Simulations using MITS demonstrate that the overall IV characteristics in a device with a random distribution of islands are a result of a complex interplay among those factors that affect the tunneling rates that are fixed a priori (e.g. island sizes, island separations, temperature, gate bias, etc.), and the evolving charge state of the system, which changes as the source-drain bias (VSD) is changed. With increasing VSD, a multi-island device has to overcome multiple discrete energy barriers (up-steps) before it reaches the threshold voltage (Vth). Beyond Vth, current flow is rate-limited by slow junctions, which leads to the CS structures in the IV characteristic. Each step in the CS is characterized by a unique distribution of island charges with an associated distribution of tunneling probabilities. MITS simulation studies done on one-dimensional (1D) disordered chains show that longer chains are better suited for switching applications as Vth increases with increasing chain length. They are also able to retain CS structures at higher temperatures better than shorter chains. In sufficiently disordered 2D systems, we demonstrate that there may exist a dominant conducting path (DCP) for conduction, which makes the 2D device behave as a quasi-1D device. The existence of a DCP is sensitive to the device structure, but is robust with respect to changes in temperature, gate bias, and VSD. A side gate in 1D and 2D systems can effectively control Vth. We argue that devices with smaller island sizes and narrower junctions may be better suited for practical applications, especially at room temperature.

Momentan verfügbare Selbstkontrollkraft: Vorstellung eines Messinstruments und erste Befunde aus pädagogisch-psychologischen Kontexten

Selbstkontrolle spielt in pädagogischen Kontexten oft eine entscheidende Rolle, beispielsweise wenn Lernende attraktiven Handlungsalternativen widerstehen müssen. Etliche Studien zeigen, dass die erfolgreiche Ausübung von Selbstkontrolle von der momentanen Verfügbarkeit an Selbstkontrollkraft abhängt. Zur Erfassung der momentanen Selbstkontrollkraft in deutschsprachigen Stichproben, wurde die State Self-Control Capacity Scale ins Deutsche adaptiert. In der vorliegenden Arbeit wird über fünf Teilstudien berichtet, in denen die psychometrische Güte dieser Skala überprüft wurde. Über alle Studien hinweg erwies sich die eindimensionale Skala als sehr reliabel (Cronbachs Alphas ≥ .93). Bezüglich ihrer Validität zeigten sich erwartete Zusammenhänge mit Validitätskriterien (z. B. Angst, Stress, am gegenwärtigen Tag bereits in Lehrveranstaltungen verbrachte Zeit); zusätzlich gelang die empirische Trennung von verschiedenen Stimmungsaspekten. Des Weiteren bildete die neue Skala (quasi-) experimentell im Labor und im Feld (Hochschule) induzierte Unterschiede in der Selbstkontrollkraft ab. Auch eine aus zehn Items bestehende Kurzskala erwies sich als hinreichend reliabel und valide.

On the similarity of Sturm-Liouville operators with non-Hermitian boundary conditions to self-adjoint and normal operators

We consider one-dimensional Schrödinger-type operators in a bounded interval with non-self-adjoint Robin-type boundary conditions. It is well known that such operators are generically conjugate to normal operators via a similarity transformation. Motivated by recent interests in quasi-Hermitian Hamiltonians in quantum mechanics, we study properties of the transformations and similar operators in detail. In the case of parity and time reversal boundary conditions, we establish closed integral-type formulae for the similarity transformations, derive a non-local self-adjoint operator similar to the Schrödinger operator and also find the associated “charge conjugation” operator, which plays the role of fundamental symmetry in a Krein-space reformulation of the problem.

Low-dimensional modeling of streaks in a wedge flow boundary layer

This paper is concerned with the low dimensional structure of optimal streaks in a wedge flow boundary layer, which have been recently shown to consist of a unique (up to a constant factor) three-dimensional streamwise evolving mode, known as the most unstable streaky mode. Optimal streaks exhibit a still unexplored/unexploited approximate self-similarity (not associated with the boundary layer self-similarity), namely the streamwise velocity re-scaled with their maximum remains almost independent of both the spanwise wavenumber and the streamwise coordinate; the remaining two velocity components instead do not satisfy this property. The approximate self-similar behavior is analyzed here and exploited to further simplify the description of optimal streaks. In particular, it is shown that streaks can be approximately described in terms of the streamwise evolution of the scalar amplitudes of just three one-dimensional modes, providing the wall normal profiles of the streamwise velocity and two combinations of the cross flow velocity components; the scalar amplitudes obey a singular system of three ordinary differential equations (involving only two degrees of freedom), which approximates well the streamwise evolution of the general streaks.

Systematic derivation of partition functions for ligand binding to two-dimensional lattices

The Ising problem consists in finding the analytical solution of the partition function of a lattice once the interaction geometry among its elements is specified. No general analytical solution is available for this problem, except for the one-dimensional case. Using site-specific thermodynamics, it is shown that the partition function for ligand binding to a two-dimensional lattice can be obtained from those of one-dimensional lattices with known solution. The complexity of the lattice is reduced recursively by application of a contact transformation that involves a relatively small number of steps. The transformation implemented in a computer code solves the partition function of the lattice by operating on the connectivity matrix of the graph associated with it. This provides a powerful new approach to the Ising problem, and enables a systematic analysis of two-dimensional lattices that model many biologically relevant phenomena. Application of this approach to finite two-dimensional lattices with positive cooperativity indicates that the binding capacity per site diverges as Na (N = number of sites in the lattice) and experiences a phase-transition-like discontinuity in the thermodynamic limit N → ∞. The zeroes of the partition function tend to distribute on a slightly distorted unit circle in complex plane and approach the positive real axis already for a 5×5 square lattice. When the lattice has negative cooperativity, its properties mimic those of a system composed of two classes of independent sites with the apparent population of low-affinity binding sites increasing with the size of the lattice, thereby accounting for a phenomenon encountered in many ligand-receptor interactions.