Numerical constraints on the model of stochastic excitation of solar-type oscillations

Analyses of a 3D simulation of the upper layers of a solar convective envelope provide constraints on the physical quantities which enter the theoretical formulation of a stochastic excitation model of solar p modes, for instance the convective velocities and the turbulent kinetic energy spectrum. These constraints are then used to compute the acoustic excitation rate for solar p modes, P. The resulting values are found ~5 times larger than the values resulting from a computation in which convective velocities and entropy fluctuations are obtained with a 1D solar envelope model built with the time-dependent, nonlocal Gough (1977) extension of the mixing length formulation for convection (GMLT). This difference is mainly due to the assumed mean anisotropy properties of the velocity field in the excitation region. The 3D simulation suggests much larger horizontal velocities compared to vertical ones than in the 1D GMLT solar model. The values of P obtained with the 3D simulation constraints however are still too small compared with the values inferred from solar observations. Improvements in the description of the turbulent kinetic energy spectrum and its depth dependence yield further increased theoretical values of P which bring them closer to the observations. It is also found that the source of excitation arising from the advection of the turbulent fluctuations of entropy by the turbulent movements contributes ~ 65-75 % to the excitation and therefore remains dominant over the Reynolds stress contribution. The derived theoretical values of P obtained with the 3D simulation constraints remain smaller by a factor ~3 compared with the solar observations. This shows that the stochastic excitation model still needs to be improved.Comment: 11 pages, 9 figures, accepted for publication in A&

The Order of Phase Transitions in Barrier Crossing

A spatially extended classical system with metastable states subject to weak spatiotemporal noise can exhibit a transition in its activation behavior when one or more external parameters are varied. Depending on the potential, the transition can be first or second-order, but there exists no systematic theory of the relation between the order of the transition and the shape of the potential barrier. In this paper, we address that question in detail for a general class of systems whose order parameter is describable by a classical field that can vary both in space and time, and whose zero-noise dynamics are governed by a smooth polynomial potential. We show that a quartic potential barrier can only have second-order transitions, confirming an earlier conjecture [1]. We then derive, through a combination of analytical and numerical arguments, both necessary conditions and sufficient conditions to have a first-order vs. a second-order transition in noise-induced activation behavior, for a large class of systems with smooth polynomial potentials of arbitrary order. We find in particular that the order of the transition is especially sensitive to the potential behavior near the top of the barrier.Comment: 8 pages, 6 figures with extended introduction and discussion; version accepted for publication by Phys. Rev.

Magnetohydrodynamic turbulence in warped accretion discs

Warped, precessing accretion discs appear in a range of astrophysical systems, for instance the X-ray binary Her X-1 and in the active nucleus of NGC4258. In a warped accretion disc there are horizontal pressure gradients that drive an epicyclic motion. We have studied the interaction of this epicyclic motion with the magnetohydrodynamic turbulence in numerical simulations. We find that the turbulent stress acting on the epicyclic motion is comparable in size to the stress that drives the accretion, however an important ingredient in the damping of the epicyclic motion is its parametric decay into inertial waves.Comment: to appear in the proceedings of the 20th Texas Symposium on Relativistic Astrophysics, J. C. Wheeler & H. Martel (eds.

The response of a turbulent accretion disc to an imposed epicyclic shearing motion

We excite an epicyclic motion, whose amplitude depends on the vertical position, $z$, in a simulation of a turbulent accretion disc. An epicyclic motion of this kind may be caused by a warping of the disc. By studying how the epicyclic motion decays we can obtain information about the interaction between the warp and the disc turbulence. A high amplitude epicyclic motion decays first by exciting inertial waves through a parametric instability, but its subsequent exponential damping may be reproduced by a turbulent viscosity. We estimate the effective viscosity parameter, $\alpha_{\rm v}$, pertaining to such a vertical shear. We also gain new information on the properties of the disc turbulence in general, and measure the usual viscosity parameter, $\alpha_{\rm h}$, pertaining to a horizontal (Keplerian) shear. We find that, as is often assumed in theoretical studies, $\alpha_{\rm v}$ is approximately equal to $\alpha_{\rm h}$ and both are much less than unity, for the field strengths achieved in our local box calculations of turbulence. In view of the smallness ($\sim 0.01$) of $\alpha_{\rm v}$ and $\alpha_{\rm h}$ we conclude that for $\beta = p_{\rm gas}/p_{\rm mag} \sim 10$ the timescale for diffusion or damping of a warp is much shorter than the usual viscous timescale. Finally, we review the astrophysical implications.Comment: 12 pages, 18 figures, MNRAS accepte

Simulations of Oscillation Modes of the Solar Convection Zone

We use the three-dimensional hydrodynamic code of Stein and Nordlund to realistically simulate the upper layers of the solar convection zone in order to study physical characteristics of solar oscillations. Our first result is that the properties of oscillation modes in the simulation closely match the observed properties. Recent observations from SOHO/MDI and GONG have confirmed the asymmetry of solar oscillation line profiles, initially discovered by Duvall et al. In this paper we compare the line profiles in the power spectra of the Doppler velocity and continuum intensity oscillations from the SOHO/MDI observations with the simulation. We also compare the phase differences between the velocity and intensity data. We have found that the simulated line profiles are asymmetric and have the same asymmetry reversal between velocity and intensity as observed. The phase difference between the velocity and intensity signals is negative at low frequencies and jumps in the vicinity of modes as is also observed. Thus, our numerical model reproduces the basic observed properties of solar oscillations, and allows us to study the physical properties which are not observed.Comment: Accepted for publication in ApJ Letter

Presentations: from Kac-Moody groups to profinite and back

We go back and forth between, on the one hand, presentations of arithmetic and Kac-Moody groups and, on the other hand, presentations of profinite groups, deducing along the way new results on both

Numerical MHD Simulations of Solar Magnetoconvection and Oscillations in Inclined Magnetic Field Regions

The sunspot penumbra is a transition zone between the strong vertical magnetic field area (sunspot umbra) and the quiet Sun. The penumbra has a fine filamentary structure that is characterized by magnetic field lines inclined toward the surface. Numerical simulations of solar convection in inclined magnetic field regions have provided an explanation of the filamentary structure and the Evershed outflow in the penumbra. In this paper, we use radiative MHD simulations to investigate the influence of the magnetic field inclination on the power spectrum of vertical velocity oscillations. The results reveal a strong shift of the resonance mode peaks to higher frequencies in the case of a highly inclined magnetic field. The frequency shift for the inclined field is significantly greater than that in vertical field regions of similar strength. This is consistent with the behavior of fast MHD waves.Comment: 9 pages, 6 figures, Solar Physics (in press

Numerical 3D constraints on convective eddy time-correlations : consequences for stochastic excitation of solar p modes

A 3D simulation of the upper part of the solar convective zone is used to obtain information on the frequency component, chi_k, of the correlation product of the turbulent velocity field. This component plays an important role in the stochastic excitation of acoustic oscillations. A time analysis of the solar simulation shows that a gaussian function does not correctly reproduce the nu-dependency of chi_k inferred from the 3D simuation in the frequency range where the acoustic energy injected into the solar p modes is important (nu ~ 2 - 4 mHz). The nu-dependency of chi_k is fitted with different analytical functions which can then conveniently be used to compute the acoustic energy supply rate P injected into the solar radial oscillations. With constraints from a 3D simulation, adjustment of free parameters to solar data is no longer necessary and is not performed here. The result is compared with solar seismic data. Computed values of P obtained with the analytical function which fits best chi_k are found ~ 2.7 times larger than those obtained with the gaussian model and reproduce better the solar seismic observations. This non-gaussian description also leads to a Reynolds stress contribution of the same order as the one arising from the advection of the turbulent fluctuations of entropy by the turbulent motions. Some discrepancy between observed and computed P values still exist at high frequency and possible causes for this discrepancy are discussed.Comment: 11 pages; 4 figures, accepted for publication in A&