SABBATICAL LEAVE REQUEST
Prepared October 13, 2008

Joel E. Tohline, Alumni Professor

Louisiana State University
Department of Physics & Astronomy, 202 Nicholson Hall,
Baton Rouge, LA 70803-4001 U.S.A.

Through this document, I am requesting permission from the LSU Board of Supervisors to take one full semester of sabbatical leave at full pay during the Spring semester of the year 2010.


1. Purpose and Objectives of the Leave

This sabbatical leave will provide an opportunity for me to become immersed in several specific activities within a research institute at an educational institution that offers a strong tradition of academic excellence and a stimulating research environment, but where I will be uninterrupted by regular classroom teaching duties and the numerous regular committee assignments that have become associated with my daily routine at LSU. My primary objectives are:

  1. To learn how to take full advantage of the rapidly expanding capabilities of modern scientific visualization tools to quantitatively analyze the results of my research group's large-scale, high-performance-computer simulations of astrophysical fluid flows.

  2. To expand considerably the content, capabilities and utility of the online, graduate-level textbook that I have been developing over the past decade at LSU, incorporating interactive scientific visualization tools and the ability to trace the provenance of numerical simulations that further our understanding of astrophysical systems.

  3. To learn how to modify my group's primary numerical simulation algorithms so that they will execute efficiently on future petascale computer architectures whose processing power will be drawn from hardware accelerators that resemble today's graphics processing units (GPUs).

Over the past 30 years (26 of which I have spent at LSU), the research activities of my group have largely been focused on gaining a better understanding of dynamically evolving astrophysical systems, such as: interacting and merging binary stars; neutron stars that form in conjunction with supernova explosions; star formation processes in our Galaxy's interstellar medium; and the gaseous disks of remote galaxies. In most of these studies, large-scale computational fluid dynamic (CFD) techniques have been used to model these astrophysical systems and scientific visualization techniques have been used to analyze the results of each CFD simulation. The lessons that I expect to learn while pursuing each of the objectives itemized above will enable my research group to (a) carry out more ambitious numerical simulations over the coming decade, (b) more thoroughly analyze the numerical results of each simulation, and (c) publish reproducible simulation results in a venue that embraces digital technologies much more fully than do existing archival journals.

In the past, my group has invested a considerable amount of effort developing its own visualization tools, borrowing from scattered sources various techniques that have adequately addressed a limited set of our data analysis needs [1]. Over the past decade, the computer science community has significantly accelerated its research efforts in the area of scientific visualization. As a result of this activity, a broad array of visualization tools has been developed and is being supported via an active "open source" community of developers and programmers. In my judgment, the time is ripe for my group to begin to take full advantage of the computer science community's new, extensive "visualization toolkit." Objective #1 has grown out of this strategic realization.

Objective #2 is closely related to the first. Since the dawn of the era of web browsers, I have been extremely frustrated by the slow pace at which the scientific community has been willing to embrace digital technologies in support of archival journal and textbook publications. When coupled with the growing popularity and versatility of "wikis," the computer science community's extensive visualization toolkit offers an excellent opportunity to significantly expand the content and value of scientific publications, especially publications that are intended to report on the results of large-scale numerical simulations. The requested leave will provide me with an opportunity to lead the computational science community into this new arena.

As alluded to in objective #3, developments in the visualization community are simultaneously touching the research efforts of my group in a third important way. The GPUs that reside in virtually all laptop and desktop PCs — primarily intended to support the execution of computationally demanding video games — are extraordinarily fast compute engines which, for certain classes of problems, can significantly outperform the PC's general-purpose, central processing unit. Over the past few years, the manufacturers of GPUs have begun to make it possible for researchers outside of the video-game industry to write programs that are tuned to a GPU's capabilities. I want to explore to what extend the most computationally demanding segments of my group's CFD simulation algorithms can be migrated to run efficiently on GPUs.


2. Outline of Proposed Activities

a.) Background

Throughout my career at LSU, a significant portion of my research activities have been focused on the development of tools (primarily in the form of efficient numerical algorithms) that will permit astronomers to accurately model the structure, stability, and dynamical evolution of rapidly rotating, (Newtonian) self-gravitating astrophysical fluid systems. Through continuous funding from the astronomy division of the National Science Foundation (NSF) over the first 15 years of my LSU career, much of this work was conducted with the expressed purpose of gaining a better understanding of the processes by which binary star systems form in our Galaxy; see my invited, 2001 review article on this topic [2]. Ten years ago, I requested [3] and was granted permission to spend one semester of sabbatical leave at the California Institute of Technology (Caltech). One stated objective accompanying that 1998 sabbatical leave request was,

"to better understand the connection between the models of (Newtonian) self-gravitating fluid system that have been constructed ... by students in my group ... and the models of relativistic systems that may be used effectively to decipher the properties of gravitational wave signals that are expected to be detected by LIGO instrumentation." The Spring 2000 semester that I spent at Caltech while on sabbatical leave was extraordinarily successful. The two research papers [4, 5] that resulted directly from a collaborative research project that was spawned in Kip Thorne's research group that spring proved to be just the tip of the iceburg. Among the collection of referreed publications that I have published since 2000, twelve have been devoted to studies of astrophysical systems that are potential sources of gravitational-wave radiation. The reputation that my group quickly established regarding its ability to produce quantitatively accurate models of such, often quite complex, astrophysical sources has permitted me to maintain a continuous stream of funding from NSF's astronomy program while switching fields (from having recognized expertise in "star formation" to having recognized expertise in the area of "gravitational-wave sources"). This also put me in an opportunistic position to collaborate effectively with Luis Lehner in landing a 4-year, multi-million dollar grant from the NSF's Information Technology Research (ITR) program in 2003. Other NSF and NASA grants have followed; over the past 5.5 years, I have either led or been a significant co-investigator on federal awards totaling just over $5.5 million.

My success in crossing from one field into another, and in nurturing successful multi-disciplinary research collaborations has been recognized in several additional ways outside of LSU:

b.) Objective #1

Observational data from arrays of millimeter-wave radio telescopes now provide sufficient spatial resolution and signal-to-noise to permit mapping of the structure and dynamical properties of star-forming gas clouds with linear scales approaching the size of our own solar system (cf., Sargent & Welch 1993; Ohashi et al. 1997). Hence, there is an opportunity to directly compare the properties of my group's most recent models with the observed properties of star-forming clouds having comparable scales. One of the most productive facilities of this type is the Owens Valley Millimeter Array which is operated by CalTech as part of the larger Owens Valley Radio Observatory (OVRO). As Director of the OVRO, Anneila Sargent leads a strong radio astronomy group at CalTech whose efforts are largely focused on studies of star forming regions in our Galaxy.

By spending an extended period of time interacting with this radio astronomy group at CalTech, I expect to gain a much better appreciation of the variety and quality of data that is being collected in connection with ongoing star formation processes in the solar neighborhood of our Galaxy as well as in extragalactic systems. While on leave, I plan to modify my group's current (generally adiabatic) modeling algorithms to incorporate radiative transfer techniques that will permit us to produce spatially resolved surface brightness and velocity maps from our dynamical models for comparison with the observations. Having gained a much better appreciation of the type (eg., mass, number density, and linear scale) of system that can be well studied with maturing millimeter-wave arrays, I will understand better on what types of systems my group should focus its future dynamical modeling efforts upon my return to LSU.

c.) Objective #2

Although directed primarily toward a better understanding of how binary star systems are formed, most of my group's recent modeling efforts have been of a sufficiently general nature that they also can provide insight into the structural and stability properties of other self-gravitating systems, such as galaxies or highly evolved binary systems with or without accompanying accretion disks. In connection with highly evolved, compact stars or compact binary systems that might give rise to observable levels of gravitational radiation, for example, my group recently has made two significant contributions ( New et al. 1995; New & Tohline 1997; see also New 1996). The second of these papers, in particular, reports on an extensive set of numerical calculations in which the relative stability of close, equal-mass binary stars having a wide variety of different equations of state was carefully examined.

As mentioned in § 2a, above, my group recently has demonstrated that it is possible to construct dynamically stable, self-gravitating configurations with highly nonaxisymmetric structures out of compressible gases. If scaled to the size of compact stars, such long-lived nonaxisymmetric configurations would produce significant, luminous sources of gravitational radiation. We are anxious to gain a better understanding of the connection between (Newtonian) models of this type and models of relativistic systems that may be detectable by LIGO instrumentation. Under Kip Thorne's direction, CalTech is widely recognized as the home base for one of the world's leading theory groups whose efforts are directed largely toward models of LIGO sources. By participating in the activities of this group on a regular basis while on leave, I should be able to achieve this second objective. Upon my return to LSU, I should be in a good position to decide whether an expansion of my group's activities in this direction will be of value to the LIGO scientific community.

d.) Objective #3

Over the past few years I have been investing a considerable amount of time developing an extensive, technical web-based document that covers the extremely broad topic of "The Structure, Stability, and Dynamics of Self-Gravitating Systems." (http://www.phys.lsu.edu/astro/H_Book.current/H_Book.shtml) My objective is to gather into a single source: much of the classical literature on the structural and stability properties of rotating, incompressible fluids that has been reviewed by Chandrasekhar (1969); much of the work on the stability of rotating stars that was presented by Tassoul (1978) two decades ago; and an overview of numerical models that have contributed to our understanding of the nonlinear development of instabilities in such systems over the past 20 years. A web-based document (as opposed to a printed textbook) is particularly appealing when presenting recent numerical developments because the reader can be given direct access to existing numerical tools, and animation sequences (or Java applets) can be called upon to illustrate clearly the behavior of time-evolving, dynamical systems.

Although I have already invested a great deal of time laying out the structural elements of this web-based document and have written a good deal of the introductory material, the "meat" of the subject has yet to be presented in a satisfactory manner. (Surprisingly, the document already has received some publicity; it was featured on the cover and in an "Internet Goldmine" article in the 1998 July/August issue of Computers in Physics.) During this proposed sabbatical leave, I expect to find many uninterrupted hours during which the primary text of this web-based document can be written and supporting utilities developed. Although my hope is that the astrophysical community at large will find this technical document to be of value in a variety of different contexts, my primary motivation for developing the document is to provide a modern resource to graduate students at LSU who wish to push forward the frontiers of knowledge in this field. I also expect the document to serve as a primary resource for the graduate-level astrophysics course (ASTR 7741) that is offered on a regular basis to physics students at LSU.


3. Location of Leave

My plans are to base my sabbatical leave activities on the campus of the California Institute of Technology (CalTech) in Pasadena, California. In concert with my planned research activities, described above, I will have two primary academic contacts at CalTech:

As documented by the accompanying letters from Professors Sargent and Thorne, I have been invited to participate in activities associated with both of their research programs in an effort to develop a more complete understanding of how the numerical, fluid dynamic simulations currently being performed by my group at LSU relate to (a) high-resolution radio observations of newly forming stars both inside and outside the solar neighborhood of our Galaxy that are being acquired with OVRO instruments, and (b) efforts that are underway to predict the luminosity and spectral signature of strong sources of gravitational radiation that will be detectable by LIGO. Professor Sargent have assured me that office space will be provided on the CalTech campus from which my research activities may comfortably be based, and that I will have opportunities to visit the OVRO.


4. Compensation

I do not expect to receive any compensation from sources other than the LSU System while on leave.


5. References
  1. Tohline, J. E., "Scientific Visualization: A necessary chore," Computing in Science & Engineering, 9, no. 6, 76-81 (2007)
    [www.phys.lsu.edu/~tohline/CiSE/CiSE2007.Vol9No6.pdf]

  2. Tohline, J. E., "The Origin of Binary Stars," Annual Reviews of Astronomy & Astrophysics, 40, 349-385 (2002)

  3. Tohline, J. E., "Sabbatical Leave Request," (1998)
    [www.phys.lsu.edu/~tohline/sabbatical/leave.brief.html]

  4. Lindblom, L., Tohline, J. E., & Vallisneri, M., "Non-Linear Evolution of the r-Modes in Neutron Stars," Physical Review Letters, 86, 1152-1155 (2001)

  5. Lindblom, L., Tohline, J. E., & Vallisneri, M., "Numerical Evolutions of Nonlinear r-Modes in Neutron Stars," Physical Review D, 65, 084039 (2002)



Andalib, S. W., J.E. Tohline, & D.M. Christodoulou. (1997). The Astrophysical Journal Supplement, 108, 471-487.

Chandrasekhar, S. (1969), Ellipsoidal Figures of Equilibrium. Yale University Press, New Haven, CT, U.S.A.

Mathieu, R.D. (1994), Ann. Rev. Astr. Ap., 32, p. 465

Ohashi, N., Hayashi, M., Ho, P.T.P., Momose, M., Tamura, M., Hirano, N., and Sargent A.I. (1997). The Astrophysical Journal, 488, 317.

New, K.C.B. (1996), Instabilities in and Gravitational Radiation from Compact Stars and Compact Binary Systems, Ph.D. Dissertation, Louisiana State University

New, K.B.C., Chanmugam, G., Johnson, W.W., and Tohline, J.E. (1995), The Astrophysical Journal, 450, p. 757

New, K.B.C. and Tohline, J.E. (1997), The Astrophysical Journal, 490, p. 311

Sargent, A.i., and Welch, W.J (1993), Ann. Rev. Astr. Ap., 31, p. 297

Tassoul, J.-L. (1978), Theory of Rotating Stars. Princeton University Press, Princeton, NJ, U.S.A.

Tohline, J.E., Cazes, J.E., and Cohl, H.S (1998), in Proceedings of the Conference, Numerical Astrophysics 1998, in press.

Woodward, J. W., J. E. Tohline, and I. Hachisu. (1994). The Astrophysical Journal, 420, 247-267.


6. Resume

Curriculum Vitae [http://www.phys.lsu.edu/faculty/tohline/resume.html]
Publications [http://www.phys.lsu.edu/faculty/tohline/ref_ref.html]