Evolution of galaxies due to self-excitation

These lectures will cover methods for studying the evolution of galaxies since their formation. Because the properties of a galaxy depend on its history, an understanding of galaxy evolution requires that we understand the dynamical interplay between all components. The first part will emphasize n-body simulation methods which minimize sampling noise. These techniques are based on harmonic expansions and scale linearly with the number of bodies, similar to Fourier transform solutions used in cosmological simulations. Although fast, until recently they were only efficiently used for small number of geometries and background profiles. These same techniques may be used to study the modes and response of a galaxy to an arbitrary perturbation. In particular, I will describe the modal spectra of stellar systems and role of damped modes which are generic to stellar systems in interactions and appear to play a significant role in determining the common structures that we see. The general development leads indirectly to guidelines for the number of particles necessary to adequately represent the gravitational field such that the modal spectrum is resolvable. I will then apply these same excitation to understanding the importance of noise to galaxy evolution.Comment: 24 pages, 7 figures, using Sussp.sty (included). Lectures presented at the NATO Advanced Study Institute, "The Restless Universe: Applications of Gravitational N-Body Dynamics to Planetary, Stellar and Galactic Systems," Blair Atholl, July 200

Time-Dependent Reaction Coefficient Identification Problems with a Free Boundary

The determination of time-dependent reaction coefficients in free boundary heat transfer problems is numerically investigated. The additional data which provides a unique solution is given by the Stefan boundary condition and the heat moments. The finite difference method with the Crank-Nicolson scheme combined with a regularized nonlinear optimization is employed. The resulting nonlinear system of equations is solved numerically using the MATLAB toolbox routine lsqnonlin for minimizing the Tikhonov regularization functional. A discussion of the choice of regularization parameters is provided. Numerical results are presented and discussed

Finite temperature hydrodynamic modes of trapped quantum gases

The hydrodynamic equations of an ideal fluid formed by a dilute quantum gas in a parabolic trapping potential are studied analytically and numerically. Due to the appearance of internal modes in the fluid stratified by the trapping potential, the spectrum of low-lying modes is found to be dense in the high-temperature limit, with an infinitely degenerate set of zero-frequency modes. The spectrum for Bose-fluids and Fermi-fluids is obtained and discussed.Comment: 26 pages, Late

Bouncing trimer: a random self-propelled particle, chaos and periodical motions

A trimer is an object composed of three centimetrical stainless steel beads equally distant and is predestined to show richer behaviours than the bouncing ball or the bouncing dimer. The rigid trimer has been placed on a plate of a electromagnetic shaker and has been vertically vibrated according to a sinusoidal signal. The horizontal translational and rotational motions of the trimer have been recorded for a range of frequencies between 25 and 100 Hz while the amplitude of the forcing vibration was tuned for obtaining maximal acceleration of the plate up to 10 times the gravity. Several modes have been detected like e.g. rotational and pure translational motions. These modes are found at determined accelerations of the plate and do not depend on the frequency. By recording the time delays between two successive contacts when the frequency and the amplitude are fixed, a mapping of the bouncing regime has been constructed and compared to that of the dimer and the bouncing ball. Period-2 and period-3 orbits have been experimentally observed. In these modes, according to observations, the contact between the trimer and the plate is persistent between two successive jumps. This persistence erases the memory of the jump preceding the contact. A model is proposed and allows to explain the values of the particular accelerations for which period-2 and period-3 modes are observed. Finally, numerical simulations allow to reproduce the experimental results. That allows to conclude that the friction between the beads and the plate is the major dissipative process.Comment: 22 pages, 10 figure

Bilayer Membrane in Confined Geometry: Interlayer Slide and Steric Repulsion

We derived free energy functional of a bilayer lipid membrane from the first principles of elasticity theory. The model explicitly includes position-dependent mutual slide of monolayers and bending deformation. Our free energy functional of liquid-crystalline membrane allows for incompressibility of the membrane and vanishing of the in-plane shear modulus and obeys reflectional and rotational symmetries of the flat bilayer. Interlayer slide at the mid-plane of the membrane results in local difference of surface densities of the monolayers. The slide amplitude directly enters free energy via the strain tensor. For small bending deformations the ratio between bending modulus and area compression coefficient, Kb/KA, is proportional to the square of monolayer thickness, h. Using the functional we performed self-consistent calculation of steric potential acting on bilayer between parallel confining walls separated by distance 2d. We found that temperature-dependent curvature at the minimum of confining potential is enhanced four times for a bilayer with slide as compared with a unit bilayer. We also calculate viscous modes of bilayer membrane between confining walls. Pure bending of the membrane is investigated, which is decoupled from area dilation at small amplitudes. Three sources of viscous dissipation are considered: water and membrane viscosities and interlayer drag. Dispersion has two branches. Confinement between the walls modifies the bending mode with respect to membrane in bulk solution. Simultaneously, inter-layer slipping mode, damped by viscous drag, remains unchanged by confinement.Comment: 23 pages,3 figures, pd