Effect of dust grain porosity on the appearance of protoplanetary disks

We theoretically analyze protoplanetary disks consisting of porous dust grains. In the analysis of observations of protoplanetary disks the dust phase is often assumed to consist of spherical grains, allowing one to apply the Mie scattering formalism. However, in reality, the shape of dust grains is expected to deviate strongly from that of a sphere. We investigate the influence of porous dust grains on the temperature distribution and observable appearance of protoplanetary disks for dust grain porosities of up to 60 %. We performed radiative transfer modeling to simulate the temperature distribution, spectral energy distribution, and spatially resolved intensity and polarization maps. The optical properties of porous grains were calculated using the method of discrete dipole approximation. We find that the flux in the optical wavelength range is for porous grains higher than for compact, spherical grains. The profile of the silicate peak at 9.7 um strongly depends on the degree of grain porosity. The temperature distribution shows significant changes in the direction perpendicular to the midplane. Moreover, simulated polarization maps reveal an increase of the polarization degree by a factor of about four when porous grains are considered, regardless of the disk inclination. The polarization direction is reversed in selected disk regions, depending on the wavelength, grain porosity, and disk inclination. We discuss several possible explanations of this effect and find that multiple scattering explains the effect best. Porosity influences the observable appearance of protoplanetary disks. In particular, the polarization reversal shows a dependence on grain porosity. The physical conditions within the disk are altered by porosity, which might have an effect on the processes of grain growth and disk evolution.Comment: 12 pages, 18 figure

Detecting planets in protoplanetary disks: A prospective study

We investigate the possibility to find evidence for planets in circumstellar disks by infrared and submillimeter interferometry. We present simulations of a circumstellar disk around a solar-type star with an embedded planet of 1 Jupiter mass. The three-dimensional (3D) density structure of the disk results from hydrodynamical simulations. On the basis of 3D radiative transfer simulations, images of this system were calculated. The intensity maps provide the basis for the simulation of the interferometers VLTI (equipped with the mid-infrared instrument MIDI) and ALMA. While MIDI/VLTI will not provide the possibility to distinguish between disks with or without a gap on the basis of visibility measurements, ALMA will provide the necessary basis for a direct gap detection.Comment: 5 page

Planar oscillatory stirring apparatus

The present invention is directed to an apparatus for stirring materials using planar orthogonal axes oscillations. The apparatus has a movable slide plate sandwiched between two fixed parallel support plates. Pressurized air is supplied to the movable slide plate which employs a tri-arm air bearing vent structure which allows the slide plate to float and to translate between the parallel support plates. The container having a material to be stirred is secured to the upper surface of the slide plate through an aperture in the upper support plate. A motor driven eccentric shaft loosely extends into a center hole bearing of the slide plate to cause the horizontal oscillations. Novelty lies in the combination of elements which exploits the discovery that low frequency, orthogonal oscillations applied horizontally to a Bridgman crucible provides a very rigorous stirring action, comparable with and more effective by an order of magnitude than the accelerated crucible rotation technique

Hybrid Pathwise Sensitivity Methods for Discrete Stochastic Models of Chemical Reaction Systems

Stochastic models are often used to help understand the behavior of intracellular biochemical processes. The most common such models are continuous time Markov chains (CTMCs). Parametric sensitivities, which are derivatives of expectations of model output quantities with respect to model parameters, are useful in this setting for a variety of applications. In this paper, we introduce a class of hybrid pathwise differentiation methods for the numerical estimation of parametric sensitivities. The new hybrid methods combine elements from the three main classes of procedures for sensitivity estimation, and have a number of desirable qualities. First, the new methods are unbiased for a broad class of problems. Second, the methods are applicable to nearly any physically relevant biochemical CTMC model. Third, and as we demonstrate on several numerical examples, the new methods are quite efficient, particularly if one wishes to estimate the full gradient of parametric sensitivities. The methods are rather intuitive and utilize the multilevel Monte Carlo philosophy of splitting an expectation into separate parts and handling each in an efficient manner.Comment: 30 pages. The numerical example section has been extensively rewritte

An investigation of radiative heat transfer in absorbing, emitting, and scattering media

Radiative heat transfer in nonisothermal emitting, absorbing, and scattering medi

Action Potential Onset Dynamics and the Response Speed of Neuronal Populations

The result of computational operations performed at the single cell level are coded into sequences of action potentials (APs). In the cerebral cortex, due to its columnar organization, large number of neurons are involved in any individual processing task. It is therefore important to understand how the properties of coding at the level of neuronal populations are determined by the dynamics of single neuron AP generation. Here we analyze how the AP generating mechanism determines the speed with which an ensemble of neurons can represent transient stochastic input signals. We analyze a generalization of the $\theta$-neuron, the normal form of the dynamics of Type-I excitable membranes. Using a novel sparse matrix representation of the Fokker-Planck equation, which describes the ensemble dynamics, we calculate the transmission functions for small modulations of the mean current and noise noise amplitude. In the high-frequency limit the transmission function decays as $\omega^{-\gamma}$, where $\gamma$ surprisingly depends on the phase $\theta_{s}$ at which APs are emitted. In a physiologically plausible regime up to 1kHz the typical response speed is, however, independent of the high-frequency limit and is set by the rapidness of the AP onset, as revealed by the full transmission function. In this regime modulations of the noise amplitude can be transmitted faithfully up to much higher frequencies than modulations in the mean input current. We finally show that the linear response approach used is valid for a large regime of stimulus amplitudes.Comment: Submitted to the Journal of Computational Neuroscienc

Two New Species of Leafblight Fungi on Kalmia Latifolia

The evergreen shrub, Kalmia latifolia L., commonly known as mountain laurel, calico bush, or sheep-kill, grows widely on rocky, acid soils in the eastern United States. Whether growing in its natural habit or in cultivation, mountain laurel appears to be equally subject to attack by fungi. The following account characterizes and discusses two of these fungi. One of them has not been described previously and additional observations have been made regarding the developmental morphology of the other one. Both pathogens are Pyrenomycetes, one a Physalospora and the other a Diaporthe. Each produces a leafblight disease. Tiny brown discolorations on young leaves characterize the early stages of attack by both organisms. These small lesions gradually enlarge and become irregular brown spots that may encompass the major portion of the leaf surface. The invaded tissues are darkest near the margins of the lesions, but a reddish zone lies between the darker border and the surrounding green tissues. Severely attacked leaves are deformed and shed prematurely. The reproductive structures of the Physalospora occur on the lower surface and begin to develop before the leaves are shed. The pycnidial stromata of the Diaporthe elevate the epidermis and caticle, and consequently produce grayish spots on the leaf surface. Both fungi continue to develop after the leaves have fallen, and since the mycelia extend beyond the margins of the lesions, perithecia ultimately may occupy most of the leaf surface. [excerpt

Mie-scattering function

Computer program for calculating normalized Mie scattering function

The location of the UK cotton textiles industry in 1838 : a quantitative analysis

We examine the geography of cotton textiles in Britain in 1838 to test claims about why the industry came to be so heavily concentrated in Lancashire. Our analysis considers both first and second nature aspects of geography including the availability of water power, humidity, coal prices, market access and sunk costs. We show that some of these characteristics have substantial explanatory power. Moreover, we exploit the change from water to steam power to show that the persistent effect of first nature characteristics on industry location can be explained by a combination of sunk costs and agglomeration effects