Abstract
In this paper, we develop an analytic theory for describing the photoacoustic wave generation from a spheroidal droplet and derive the first complete analytic solution. Our derivation is based on solving the photoacoustic Helmholtz equation in spheroidal coordinates with the separation-of-variables method. As the verification, besides carrying out the asymptotic analyses which recover the standard solutions for a sphere, an infinite cylinder and an infinite layer, we also confirm that the partial transmission and reflection model previously demonstrated for these three geometries still stands. We expect that this analytic solution will find broad practical uses in interpreting experiment results, considering that its building blocks, the spheroidal wave functions (SWFs), can be numerically calculated by the existing computer programs.
© 2014 Optical Society of America
1. Introduction
As fundamental problems in photoacoustics, the photoacoustic wave generation in one dimension from an infinite layer, in two dimensions from an infinite cylinder, and in three dimensions from a sphere have all been extensively studied [1] and their analytical solutions under the optically-thin-body condition are well-known [2,3]. The analytical solutions dutifully convey the characteristic spectral and temporal information embedded in the photoacoustic radiation, which are tightly correlated to the specific geometries of these bodies.
In this paper, we embark on the derivation of the complete analytic solution in three dimensions for a droplet with a more general shape other than sphere—a spheriodal droplet which can be either a prolate spheroid or an oblate spheroid. To the best of our knowledge, such theoretical investigation has not been taken account before. While such an analytic solution should be broadly applicable, we anticipate that it will first find important applications in biomedical photoacoustic imaging, considering the rapid progress in this research field to detect individual cells or organelles [4–8] and that these light absorbing biological capsules such as red blood cells and cell nuclei are close to spheroids in shape.
Our derivation is based on solving the photoacoustic Helmholtz equation by employing the separation-of-variables method in spheroidal coordinates. Working with the spheroidal coordinate systems (see Fig. 1, including the prolate spheroidal coordinate system and the oblate spheroidal coordinate system) is a natural choice because a spheroidal droplet can be fitted into such a typical coordinate system and the boundary conditions can be easily defined and applied. The resulted analytic formulas are expressed as the sum of a serials of with the weighting coefficients to be determined by the boundary conditions, where and are defined respectively as the angular and radial spheroidal wave functions (SWFs) [9].
The method of solving a Helmholtz equation in spheroidal coordinates has its roots in the problems of sound scattering as well as light scattering by a spheroidal particle [10,11], and these researches still undergo some updates [12–15]. Due to the specific rotational symmetry of the oblate spheroidal coordinates, this method also incurs new applications in describing paraxial Gaussian-Laguerre laser beam [16,17] and in calculating the first-order correctness of Fabry-Perot-resonator eigenmodes [18]. We find that introducing this approach into photoacoustics is very attractive because only the modes contribute to the analytic solution due to the fact that the photoacoustic wave generation does not depend on the laser incident direction relative to the spheroidal droplet under the conventional optically-thin-body assumption. This is contrary to the cases of sound or light scattering where this simplicity only exists in the axisymmetric solution when the sound or light incident direction is restricted to being along the rotational axis of the spheroid.
In order to validate the analytic solution we reached, we present two extended theoretical examinations. The first is the asymptotic analyses of the analytic solution under three extreme conditions: when the interfocal distance of a prolate spheroid and an oblate spheroid diminishes to zero, when the interfocal distance of a prolate spheroid stretches to infinity, and when the interfocal distance of an oblate spheroid expands to infinity. It is found that the analytic solution asymptotically degenerates respectively to the standard formulas for a sphere, an infinite cylinder and an infinite layer. The second is the decomposition analysis on the time domain solution from the perspective of photoacoustic wave partial transmission and reflections at the droplet boundary. We demonstrate that the photoacoustic pulse outside the spheroidal droplet can be decomposed to successively transmitted waves which experience multiple internal reflections before passing through the boundary, while the photoacoustic pulse formed inside can be decomposed to successively reflected waves plus the source terms. This type of decomposition is consistent with the physics picture previously revealed in the studies of the three special geometries [2].
Around these major interrogations, this paper is organized as follows. Section 2 details the derivation of the general solution, first in the frequency domain and then in the time domain. Section 3 then analyzes the asymptotic behaviors of the solution by using the asymptotic properties of angular and radial SWFs. Section 4 along with Appendix provides the proof for how to fit the time domain solution into the frame of the partial transmission and reflection model where the Wronskian relation of radial SWFs plays an essential role. Finally, Section 5 summarizes our major findings and also directs the future works.
2. General solution
The problem is to resolve the photoacoustic waves generated by a spheroidal droplet due to laser illumination, as illustrated in Fig. 1. The light absorbing droplet is either a prolate spheroid or an oblate spheroid which is immersed in a fluid absorbing little light. The laser beam can be in the form of intensity-modulated continuous wave (cw) or short pulse. Because of the laser heating, a temperature change and subsequently a pressure variation will be introduced both inside the droplet and in the surrounding fluid. Putting into formulas, this process is governed by the following two equations for nonviscous fluid [1]:
where is the thermal conductivity, is the mass density, is the specific heat capacity at constant pressure, is the heat energy per unit volume and time deposited by the laser, is the sound speed (Only the longitudinal sound speed is considered here since for nonvisous fluid the effect of the shear sound speed can be neglected.), and is the pressure expansion coefficient at constant volume defined as . The pressure variation is the photoacoustic (PA) wave. In Fig. 1, the photoacoustic wave outside of the droplet is shown.Given that the thermal conduction term is usually small compared to , its effect can be neglected. This treatment is valid not only for the nanosecond laser pulse typically employed [19] but also for the cw laser with its modulation frequency in the ultrasound range [20]. Therefore, the photoacoustic wave functions of Eq. (1) and Eq. (2) can be reduced to [1]
where is the thermal expansion coefficient related to as . To explicitly express that the laser heating happens only in the droplet and the sound speeds in the droplet and in the surrounding fluid are generally different, we divide Eq. (3) intowhere the subscript “s” denotes quantities of the spheroidal droplet while “f” denotes those of the surrounding fluid.The solution of in Eq. (4) can be further divided into two parts:
where is the special solution forand satisfiesThis solving method is appropriate because the quantity of is uniform across the whole droplet under the “optically-thin-body” assumption. This assumption also implies that the photoacoustic wave generation does not depend on the laser beam illumination angle relative the spheroidal droplet.Since it is trivial to solve in Eq. (6), the problem is thus focused on to solve the two Helmholtz equations: one in Eq. (4) for , and the other Eq. (7) for . We coin them together hereafter the “photoacoustic Helmholtz equation”.
2.1 Frequency domain solution
The frequency domain solution corresponds to the photoacoustic waves produced by a cw laser beam with the intensity function of where represents the intensity modulation frequency, under which the heating function is
where represents the light absorption coefficient of the droplet, and the percentage of the absorbed light energy being converted to heat which takes account the possibility of other energy-consumption routes such as fluorescent emission.The frequency domain solution of corresponding to Eq. (6) can be obtained easily by substituting of Eq. (8) into Eq. (6) and taking out the time dependent part , as
To get the frequency domain solution of from Eq. (7) and from Eq. (4), the standard separation-of-variables method of solving a Helmholtz equation in spheroidal coordinates can be applied [9], resulting in Here the dimensionless variables and are defined aswhere is the interfocal distance of the spheroid droplet (therefore the angular SWFs and radial SWFs all depend on ), and are the weighting coefficients to be determined by applying boundary conditions. We have used the branch brackets to work as the denotation symbol to differentiate the variables in angular SWFs and the radial SWFs and (where are called the first kind and the third kind) depending on whether these SWFs are considered in the prolate spheroidal coordinate system or in the oblate spheroidal coordinate system: when switching between the former system to the latter system, the replacements of ( represents either of and ) and are needed. We will use this type of expression throughout this paper.As plotted in Fig. 1, the boundary of the spheroidal droplet fits to a spheroidal surface with a constant , i.e. . With the conventional notation of and respectively as the semi-major and semi-minor axes of the revolving ellipse that forms the spheroid, the expression for determining will be different depending on the droplet is a prolate spheroid or an oblate spheroid, respectively as
These expressions can be directly extracted from the following coordinate transformation relations between spheroidal coordinates and rectangular coordinates (by looking at the z axis along which ):The boundary conditions that shall be applied here is identical to that employed in the classic solutions [1,2], and is just the same as that applied to solve the problem of ultrasound scattering from a penetrable spheroid [21], namely, that the pressure and the normal acceleration must both be continuous across the boundary, i.e.
As shown in Fig. 1, since the confocal hyperbolas are orthogonal to the confocal ellipses, the normal direction of the boundary is actually along the direction of the hyperbolas. It can be deduced from Eq. (15) and Eq. (16) thatwhere is the polar angle defined in the spherical coordinates. The spheroidal coordinate is defined in the conventional sense, as the azimuthal angle.Based on Eq. (17), before going for the solution of and needed in Eq. (10) and Eq. (11), we first arrive at the following constraint conditions
These constraints can be understood more intuitively from the symmetry consideration. The constraint condition corresponds to the requirement that and should be both rotationally invariant with respect to which meanswhile the constraint condition corresponds to the parity requirement of these quantities respect to which iswhich further requests that in Eq. (10) and Eq. (11). However, the angular SWFs are related to the Legendre functions as [9]where are called the spheroidal coefficients and the prime over the summation sign here indicates that the sum is over only even or odd depending on is even or odd, thereby finally resulting in the condition that must be an even number, i.e. . Therefore, the expressions for and in Eq. (10) and Eq. (11) are much simplified, being reduced toNow, it is the time to applying the boundary conditions by substituting the expressions in Eq. (9), Eq. (23) and Eq. (24) into Eq. (17). We further substitute the expansion for as shown in Eq. (22) and remove the dependence on by employing the orthogonal property of the Legendre function, subsequently obtaining
where and both denote the first derivative with respect to . Note that Eq. (25) actually contains a series of such linear equations with the total number determined by the dimension of . It will be much more convenient to cast these linear equations into a square-matrix form by setting the dimension of equal to that of , thereby arriving atHere, all of , , and are column-vectors where the elements of and are the weighting coefficients and to be determined, and the element of is simply defined as (here is the Kronecker delta symbol)On the other hand, all of , , , , , represent square-matrixes with their elements predetermined respectively as: All of , , , are diagonal-matrixes.In principle, , can be solved from Eq. (26) with matrix algebra, resulting in
To this end, by incorporating the results of Eq. (31) and Eq. (32) respectively back into Eq. (23) and Eq. (24), and further combining this result with Eqs. (9) and (5), the frequency domain solution for and in Eq. (4) is then finally determined.2.2 Time domain solution
Once the frequency domain solution has been settled, the time domain solution for short laser pulse excitation can be readily reached by taking Fourier transform based on the spectral analysis method [20].
Let us first consider decomposing the laser pulse intensity into its spectral components as
Then the time domain solution is just where has the form as shown in Eq. (9) and and have already been determined as shown in Eq. (31) and Eq. (32).For the ideal delta-pulse excitation, there is
under which Eq. (34) and Eq. (35) can be further simplified.3. Asymptotic analyses
As illustrated in the right panel of Fig. 1, when the interfocal distance , the prolate spheroid and the oblate spheroid will be both rounded respectively to a sphere; on the other hand, when , the prolate spheroid will be elongated to an infinitely long cylinder, while the oblate spheriod will be stretched to an infinitely wide circular-plate. Therefore, it is a natural conjecture that the general frequency domain solution derived in Section 2.1 (thus also the time domain solution) will asymptotically degenerate into the standard solutions for a sphere, an infinite cylinder, and an infinite layer respectively at these extreme limits. At follows, we shall provide the proofs. For the sake of brevity, here we only describe the analyses for the photoacoustic waves generated outside the droplet.
3.1 Sphere
For a prolate spheriod or an oblate spheroid, the condition of is equivalent to the requests that the dimensionless variables reach the limit of [refer to Eq. (12)]
Since there is [9]we can first deduce that [refer to Eq. (22)]Therefore, , in Eq. (28) are both reduced to the unit matrixBy substituting the results of Eq. (40) and Eq. (27) into Eq. (32), it is straightforward to reach
which means that only the single mode contributes.Now let us focus on determining the value of , , , which are related to the radial SWFs as defined in Eq. (29) and Eq. (30). We need to apply the following asymptotic relations [9]
wherehas been used, which comes from the coordinate transformation in Eq. (15) and Eq. (16) by setting . Here, is the spherical Bessel function and is the spherical Hankel function of the first kind. The asymptotic expression of and can be derived in the similar way, after which we can move a further step from Eq. (41) to getAs the final step, by substituting Eq. (45) into Eq. (24) and by noting the expression of in Eq. (38) and the expression of in Eq. (43), we arrive at
If we write out the expression for ,,, and, the standard solution for a sphere [1] will be recovered:where and .3.2 Infinite long cylinder
For a prolate spheroid, the condition of is equivalent to the requests that the dimensionless variables reach the limit of [refer to Eq. (12)]
This will lead towhich will in turn reduce the expression of in Eq. (32) toFor this time, we need to use the following asymptotic forms for the radial SWFs [22]:
where is zero-order Bessel function and is zero-order Hankel function of the first kind.Before moving forwards, we need first to analyze the asymptotic behavior of the angular SWFs which is [9]
where is the Hermite polynomials and presents constants only depending on . From Eq. (53), we can conclude that the solution of Eq. (24) with physical meaning will only exist fori.e. the photoacoustic waves are generated only along the direction perpendicular to the revolution axis (the z axis, refer to Fig. 1). This is amount to reducing the problem down to 2 dimensions. Let r denote the distance perpendicular to the z axis, we obtain from Eq. (15) thatNow by incorporating Eq. (55) into Eq. (51) and Eq. (52), and then by using them (also the expressions for and ) in Eq. (50), we reach that
Finally, by going back to Eq. (24) we can get
However, due tothe standard solution for an infinite cylinder [1] is therefore recovered with b representing the radius of the infinite cylinder:3.3 Infinite large layer
As the last asymptotic analysis, we consider the oblate spheroid with . In this case the dimensionless variables also satisfy
The results expressed same as Eq. (49) and Eq. (50) will be reached. However, it is now working under the oblate spheridal coordinate system rather than in the prolate spheroidal coordinate system. The asymptotic forms of the radial SWFs will not be the same as those in Eq. (51) and Eq. (52), but rather are [9]Here, the asymptotic expression for the angular SWFs to be used is [9]
where and are the Laguerre polynomials, and are constants only depending on . In order to have a solution with physical meaning, Eq. (63) requires thatThis is clearly different from that of Eq. (54). The constraint of Eq. (64) reduces the problem in the oblate spheroidal coordinate system down to 1 dimension (refer to Fig. 1), i.e. the photoacoustic waves are produced only along the direction of z axis. From the coordinate transformation of Eq. (16), it can be deduced thatNow, by incorporating the results of Eq. (61), Eq. (62) and Eq. (65) into the expression of [which is same as Eq. (50)], we can further arrive at
where denotes the diagonal matrix with the element .Finally, by substituting Eq. (66) and Eq. (62) into Eq. (24), and by processing the similar derivation of Eq. (58), as expected, we recover the standard solution for an infinite layer [1] with b representing here the half thickness of the infinite layer:
4. Decomposition analysis on time domain solution
As for the time domain solutions of a sphere, an infinite cylinder and an infinite layer, the model of photoacoustic wave partial transmission and reflection at the boundary has been applied successfully to explain the characteristic photoacoustic temporal responses, which appear as being composed by successive pulses [1,2]. After having demonstrated in Section 3 that the solution of a spheroidal droplet actually contains all of the standard solutions for these three special geometries as the asymptotic representations, we expect that this model should also be applicable to the time domain solution of a spheroidal droplet, which is illustrated in Fig. 2.
It is very interesting to find that for the ideal laser pulse excitation usually considered for this model [refer to Eq. (36)], the solution of and in Eq. (35) and Eq. (34) indeed can be decomposed respectively into (refer to Appendix for the derivation)
Here, , represent respectively the transmission matrix and the reflection matrix as defined in Appendix [Eq. (81)], and , are respectively defined as the initial vector and the initial vector which are calculated respectively from Eq. (32) and Eq. (31) under the conditions of (i.e. when the surrounding medium has the same acoustic properties as the spheroidal droplet under which will be a null matrix and ). The expressions of , areNow based on Eq. (68), the photoacoustic wave generation outside of the spheroidal droplet can be interpreted as follows. For each mode, the initial outgoing wave excited inside the droplet will undergo a series of successive partial transmissions and reflections. Integrating the processes of all of the modes completes the time domain solution. Note that since , are usually not diagonal matrixes, the mode hopping phenomenon will appear during partial transmission and reflection: one outgoing mode will smear into other modes. This is not unexpected, however, because in spheroidal coordinates the modes are not orthogonal with each other.
The photoacoustic wave generation inside the spheroidal droplet can be interpreted similarly based on Eq. (69). The first term and the second term on the right side of this equation assemble the initial photoacoustic pulse generated across the spheroidal droplet which is generated even when there is no reflection at the boundary. The third term shows that each mode will undergo multiple reflections. It is worth to point out that the factor of 2 inside the integral of the third term is important. This is because this factor needs to be combined with to fulfill
where the radial SWF and (the fourth kind) are necessary respectively for describing outgoing waves and ingoing waves ([9], also refer to Appendix). Only through this step, the interpretation from the perspective of partial transmission and reflection can be justified: for the photoacoustic wave generation inside the spheroidal droplet, we need to consider not only the successive reflections for an initially outgoing wave but also the successive reflections for an initially ingoing wave; the initially ingoing wave will change to an outgoing wave when reaching the inner most area corresponding to the minimum ( for a prolate spheroid and for an oblate spheroid). The mode hopping also prevails during the reflection.5. Conclusions and future works
We derived the analytic solution for describing the photoacoustic waves generated by a spheroidal droplet through two crucial steps. The first is to identify the photoacoustic Helmholtz equation, and the second is to solve it in spheroidal coordinates with separation-of-variable method which yields the expressions of Eq. (23) and Eq. (24). These expressions take the form of spheroidal modes summation with the weighting coefficients determined by Eq. (31) and Eq. (32), which are obtained by applying the boundary conditions as expressed in Eq. (17). We started our derivation in the frequency domain where the special solution of Eq. (9) also needs to be included into the complete solution. The time domain solution then can be established through Fourier transform as shown in Eq. (34) and Eq. (35).
Although the number of modes involved is theoretically unbounded and the compacted expressions of Eq. (31) and Eq. (32) in matrix format still seem complex, we were able to prove the analytic solution based on two approaches: one is to recover the standard solutions from the asymptotic analyses, and the other is to reconcile the partial transmission and reflection model by performing the decomposition analysis on the time domain solution. The former task was accomplished by exploiting the asymptotic behaviors of SWFs where the radial SWFs are asymptotically degenerated into spherical Bessel functions, Bessel functions, and sinusoidal functions which are required respectively in the solution for a sphere, an infinite cylinder and an infinite layer [23]. The key for accomplishing the latter task is to extensively apply the Wronskian relation for radial SWFs.
The information contained in the analytic solution is by no means limited by what we have demonstrated, and there is certainly more yet waiting to be explored. For example, the approximation solutions in three dimensions for a finite cylinder and for a finite plate should be reachable also through asymptotic analyses. Even more information can be extracted from the numerical realization of the analytic solution for practical uses. This is greatly promising because the hurdle of obtaining numerical results of SWFs has already been overcome with the ongoing development of the specialized computer programs [24–27]. We are currently working on this numerical realization.
Appendix: derivation for the formulas of photoacoustic wave partial transmission and reflection at the boundary of a spheroidal droplet
We perform the derivation in the frequency domain by considering an initial outgoing photoacoustic wave with the weighting coefficients of every mode all equal to 1, i.e.
where we use the vector to represent the weighting coefficients. Similar to establishing the partial transmission and reflection model for a sphere, an infinite cylinder, and an infinite layer [1,2], we seek the transmitted wave and reflected wave at the droplet boundary (refer to Fig. 2)which satisfying the boundary conditions [similar to Eq. (17)]:Here represent the radial SWFs of the fourth kinds and are necessary for describing an ingoing wave [9]. The relation between the different kinds of radial SWFs arewhere are the second kind.By using the method similar for obtaining Eq. (26), we find
where the matrixes and have been defined in Eq. (28) and the matrixes and have been defined in Eq. (30), and the matrixes ,,, are defined in the similar way, as Solving Eq. (77) with matrix algebra, we getwhere we define as the transmission matrix and the reflection matrix which have the forms as1. Expressing in terms of and
To construct the expected relation between in Eq. (32) and , we first derive the expression for . Assisted with the relation of Eq. (72), we can obtain from Eq. (81) that
Furthermore, assisted with the Wronskian relation [9]we can write asNext, by noting the interesting relation
where the left side has undergone an interchange between and and also an interchange between and (which can be easily verified also by applying the Wronskian relation), we can move a further step to reachNow, comparing Eq. (86) with Eq. (32), we can arrive at
However, noting that in Eq. (70) is just equal to (also using the Wronskian relation)we then finally achieve2.Expressing in terms of
For this task, we work on the expression of which can be obtained from Eq. (81) and Eq. (82) as that
We continue on by substituting the following expression (obtained by using the Wronskian relation)into Eq. (90), and can further reachNow, by comparing Eq. (92) with Eq. (31), and by noting that in Eq. (71) is equal to (using the Wronskian relation again)
we can finally arrive atAcknowledgments
This work is supported by the National Nature Science Foundation of China (Grant No. 11074134, 61138003), the Program for New Century Excellent Talents in University (No. NCET-10-0502), the Project Sponsored by the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry (No. 1139), and the Tianjin Municipal Science and Technology Commission under Grand number 14JCYBJC16600.
References and links
1. G. J. Diebold, “Photoacoustic monopole radiation: waves from objects with symmetry in one, two, and three dimensions,” in Photoacoustic Imaging and Spectroscopy, L. V. Wang, ed. (Taylor and Francis, 2009).
2. G. J. Diebold, M. I. Khan, and S. M. Park, “Photoacoustic “signatures” of particulate matter: optical production of acoustic monopole radiation,” Science 250(4977), 101–104 (1990). [CrossRef] [PubMed]
3. G. J. Diebold, T. Sun, and M. I. Khan, “Photoacoustic monopole radiation in one, two, and three dimensions,” Phys. Rev. Lett. 67(24), 3384–3387 (1991). [CrossRef] [PubMed]
4. L. V. Wang and S. Hu, “Photoacoustic tomography: in vivo imaging from organelles to organs,” Science 335(6075), 1458–1462 (2012). [CrossRef] [PubMed]
5. S. Hu and L. V. Wang, “Optical-resolution photoacoustic microscopy: auscultation of biological systems at the cellular level,” Biophys. J. 105(4), 841–847 (2013). [CrossRef] [PubMed]
6. E. I. Galanzha and V. P. Zharov, “Circulating tumor cell detection and capture by photoacoustic flow cytometry in vivo and ex vivo,” Cancers 5(4), 1691–1738 (2013). [CrossRef] [PubMed]
7. E. M. Strohm, E. S. L. Berndl, and M. C. Kolios, “Probing red blood cell morphology using high-frequency photoacoustics,” Biophys. J. 105(1), 59–67 (2013). [CrossRef] [PubMed]
8. E. M. Strohm, E. S. L. Berndl, and M. C. Kolios, “High frequency label-free photoacoustic microscopy of single cells,” Photoacoustics 1(3–4), 49–53 (2013). [CrossRef]
9. C. Flammer, Spheroidal Wave Functions (Stanford University, 1957).
10. R. D. Spence and S. Granger, “The scattering of sound from a prolate spheroid,” J. Acoust. Soc. Am. 23(6), 701–706 (1951). [CrossRef]
11. S. Asano and G. Yamamoto, “Light scattering by a spheroidal particle,” Appl. Opt. 14(1), 29–49 (1975). [CrossRef] [PubMed]
12. J. P. Barton, N. L. Wolff, H. Zhang, and C. Tarawneh, “Near-field calculations for a rigid spheroid with an arbitrary incident acoustic field,” J. Acoust. Soc. Am. 113(3), 1216–1222 (2003). [CrossRef] [PubMed]
13. B. R. Rapids and G. C. Lauchle, “Vector intensity field scattered by a rigid prolate spheroid,” J. Acoust. Soc. Am. 120(1), 38–48 (2006). [CrossRef]
14. M. J. Mendes, I. Tobias, A. Marti, and A. Luque, “Near-field scattering by dielectric spheroidal particles with sizes on the order of the illuminating wavelength,” J. Opt. Soc. Am. B 27(6), 1221–1231 (2010). [CrossRef]
15. M. J. Mendes, I. Tobías, A. Martí, and A. Luque, “Light concentration in the near-field of dielectric spheroidal particles with mesoscopic sizes,” Opt. Express 19(17), 16207–16222 (2011). [CrossRef] [PubMed]
16. K. T. McDonald, “Gaussian laser beams via oblate spheroidal waves,” arXiv:physics/0312024v1 (2003).
17. M. Zeppenfeld, “Solutions to Maxwell's equations using spheroidal coordinates,” New J. Phys. 11(7), 073007 (2009). [CrossRef]
18. M. Zeppenfeld and P. W. H. Pinkse, “Calculating the fine structure of a Fabry-Perot resonator using spheroidal wave functions,” Opt. Express 18(9), 9580–9591 (2010). [CrossRef] [PubMed]
19. L. V. Wang and H.-I. Wu, Biomedical Optics (Wiley, 2007), Chap. 12.
20. V. E. Gusev and A. A. Karabutov, Laser Optoacoustics (American Institute of Physics Press, 1993), Chap. 2.
21. A. D. Kotsis and J. A. Roumeliotis, “Acoustic scattering by a penetrable spheroid,” Acoust. Phys. 54(2), 153–167 (2008). [CrossRef]
22. W. M. John, “Asymptotic approximations for prolate spheroidal wave functions,” Stud. Appl. Math. 54, 315–349 (1975).
23. M. I. Khan, T. Sun, and G. J. Diebold, “Photoacoustic waves generated by absorption of laser radiation in optically thin layers,” J. Acoust. Soc. Am. 93(3), 1417–1425 (1993). [CrossRef]
24. S. Zhang and J. Jin, Computation of Special Functions (Wiley, 1997).
25. W. J. Thompson, “Spheroidal wave functions,” Comput. Sci. Eng. 1, 84–87 (1999).
26. L.-W. Li, X.-K. Kang, and M.-S. Leong, Spheroidal Wave Functions in Electromagnetic Theory (Wiley, 2002).
27. P. E. Falloon, P. C. Abbott, and J. B. Wang, “Theory and computation of spheroidal wave functions,” J. Phys. Math. Gen. 36(20), 5477–5495 (2003). [CrossRef]