Diversification and the Taxation of Capital Gains and Losses

Current U.S. law nets the total portfolio of realized capital gains and losses to compute capital gains taxes. Prior research, however, typically ignores the implication of this provision, i.e., the marginal tax rate for a specific gain or loss depends on the taxpayer's total portfolio of realized gains and losses. We find that these nettings introduce complexity into the relation between share values and capital gains taxes, creating an incentive to diversify. For firms with stock returns that are positively (negatively) correlated with those of the overall market, share values generally are decreasing (increasing) in the capital gains tax rate.

Multi-frequency fine resolution imaging radar instrumentation and data acquisition

Development of a dual polarized L-band radar imaging system to be used in conjunction with the present dual polarized X-band radar is described. The technique used called for heterodyning the transmitted frequency from X-band to L-band and again heterodyning the received L-band signals back to X-band for amplification, detection, and recording

Low Mach Number Modeling of Type Ia Supernovae

We introduce a low Mach number equation set for the large-scale numerical simulation of carbon-oxygen white dwarfs experiencing a thermonuclear deflagration. Since most of the interesting physics in a Type Ia supernova transpires at Mach numbers from 0.01 to 0.1, such an approach enables both a considerable increase in accuracy and savings in computer time compared with frequently used compressible codes. Our equation set is derived from the fully compressible equations using low Mach number asymptotics, but without any restriction on the size of perturbations in density or temperature. Comparisons with simulations that use the fully compressible equations validate the low Mach number model in regimes where both are applicable. Comparisons to simulations based on the more traditional anelastic approximation also demonstrate the agreement of these models in the regime for which the anelastic approximation is valid. For low Mach number flows with potentially finite amplitude variations in density and temperature, the low Mach number model overcomes the limitations of each of the more traditional models and can serve as the basis for an accurate and efficient simulation tool.Comment: Accepted for publication in the Astrophysical Journal 31 pages, 5 figures (some figures degraded in quality to conserve space

MAESTRO, CASTRO, and SEDONA -- Petascale Codes for Astrophysical Applications

Performing high-resolution, high-fidelity, three-dimensional simulations of Type Ia supernovae (SNe Ia) requires not only algorithms that accurately represent the correct physics, but also codes that effectively harness the resources of the most powerful supercomputers. We are developing a suite of codes that provide the capability to perform end-to-end simulations of SNe Ia, from the early convective phase leading up to ignition to the explosion phase in which deflagration/detonation waves explode the star to the computation of the light curves resulting from the explosion. In this paper we discuss these codes with an emphasis on the techniques needed to scale them to petascale architectures. We also demonstrate our ability to map data from a low Mach number formulation to a compressible solver.Comment: submitted to the Proceedings of the SciDAC 2010 meetin

Low Mach Number Modeling of Type Ia Supernovae. II. Energy Evolution

The convective period leading up to a Type Ia supernova (SN Ia) explosion is characterized by very low Mach number flows, requiring hydrodynamical methods well-suited to long-time integration. We continue the development of the low Mach number equation set for stellar scale flows by incorporating the effects of heat release due to external sources. Low Mach number hydrodynamics equations with a time-dependent background state are derived, and a numerical method based on the approximate projection formalism is presented. We demonstrate through validation with a fully compressible hydrodynamics code that this low Mach number model accurately captures the expansion of the stellar atmosphere as well as the local dynamics due to external heat sources. This algorithm provides the basis for an efficient simulation tool for studying the ignition of SNe Ia.Comment: 30 pages; accepted to the Astrophysical Journa

Direct Numerical Simulations of Type Ia Supernovae Flames II: The Rayleigh-Taylor Instability

A Type Ia supernova explosion likely begins as a nuclear runaway near the center of a carbon-oxygen white dwarf. The outward propagating flame is unstable to the Landau-Darrieus, Rayleigh-Taylor, and Kelvin-Helmholtz instabilities, which serve to accelerate it to a large fraction of the speed of sound. We investigate the Rayleigh-Taylor unstable flame at the transition from the flamelet regime to the distributed-burning regime, around densities of $10^7$ g/cc, through detailed, fully resolved simulations. A low Mach number, adaptive mesh hydrodynamics code is used to achieve the necessary resolution and long time scales. As the density is varied, we see a fundamental change in the character of the burning--at the low end of the density range the Rayleigh-Taylor instability dominates the burning, whereas at the high end the burning suppresses the instability. In all cases, significant acceleration of the flame is observed, limited only by the size of the domain we are able to study. We discuss the implications of these results on the potential for a deflagration to detonation transition.Comment: submitted to ApJ, some figures degraded due to size constraint

On the Hydrodynamic Interaction of Shock Waves with Interstellar Clouds. II. The Effect of Smooth Cloud Boundaries on Cloud Destruction and Cloud Turbulence

The effect of smooth cloud boundaries on the interaction of steady planar shock waves with interstellar clouds is studied using a high-resolution local AMR technique with a second-order accurate axisymmetric Godunov hydrodynamic scheme. A 3D calculation is also done to confirm the results of the 2D ones. We consider an initially spherical cloud whose density distribution is flat near the cloud center and has a power-law profile in the cloud envelope. When an incident shock is transmitted into a smooth cloud, velocity gradients in the cloud envelope steepen the smooth density profile at the upstream side, resulting in a sharp density jump having an arc-like shape. Such a ``slip surface'' forms immediately when a shock strikes a cloud with a sharp boundary. For smoother boundaries, the formation of slip surface and therefore the onset of hydrodynamic instabilities are delayed. Since the slip surface is subject to the Kelvin-Helmholtz and Rayleigh-Taylor instabilities, the shocked cloud is eventually destroyed in $\sim 3-10$ cloud crushing times. After complete cloud destruction, small blobs formed by fragmentation due to hydrodynamic instabilities have significant velocity dispersions of the order of 0.1 $v_b$, where $v_b$ is the shock velocity in the ambient medium. This suggests that turbulent motions generated by shock-cloud interaction are directly associated with cloud destruction. The interaction of a shock with a cold HI cloud should lead to the production of a spray of small HI shreds, which could be related to the small cold clouds recently observed by Stanimirovic & Heiles (2005). The linewidth-size relation obtained from our 3D simulation is found to be time-dependent. A possibility for gravitational instability triggered by shock compression is also discussed.Comment: 62 pages, 16 figures, submitted to Ap

The Physics of Flames in Type Ia Supernovae

We extend a low Mach number hydrodynamics method developed for terrestrial combustion, to the study of thermonuclear flames in Type Ia supernovae. We discuss the differences between 2-D and 3-D Rayleigh-Taylor unstable flame simulations, and give detailed diagnostics on the turbulence, showing that the kinetic energy power spectrum obeys Bolgiano-Obukhov statistics in 2-D, but Kolmogorov statistics in 3-D. Preliminary results from 3-D reacting bubble calculations are shown, and their implications for ignition are discussed.Comment: To appear in the Proceedings of the SciDAC 2005 meeting, IOP press (http://www.iop.org). Some figures degraded in quality to conserve spac