Quantitative Study of the Nonlinearly Enhanced Photoacoustic/Photothermal Effect by Strong LSPR-Coupled Nanoassemblies

The extensive exploration of the collective optical and thermal effects for localized surface plasmon resonance (LSPR)-coupled nanoassemblies has propelled much recent research and development in fields of photoacoustic (PA) imaging and photothermal (PT) therapy, while the rational design and proper engineering of these assemblies under quantitative guidance is still a highly challenging task. In this work, by utilizing the finite element analysis (FEA) method and taking gold nanochains as example, the authors quantitatively studied the coupling optical/thermal response of the nanoassemblies and the associated nonlinearly enhanced PA/PT effect. Results show that compared with their individuals, the strong electromagnetic/thermal coupling between the individuals of the nanoassemblies results in a several-time enhancement of the per-particle-weighted optical absorption, consequential thermal field enhancement, and initial PA pressure, resulting in nonlinearly amplified energy conversion from incident light to heat and PA waves. The dependence of the nonlinear PA/PT enhancement on the assembly chain length, the size of the individuals, the interparticle distance, and the size uniformity of the building blocks is quantitatively discussed. PA experiments on gold nanochains and gold nanospheres are performed to validate the proposition, and the experiments well silhouetted the theoretical discussion. This work paves the way for the rational construction and optimization of plasmonic nanoassemblies with improved PA/PT conversion efficiency.


Introduction
Subwavelength metallic nanoparticles have been intensively researched for the past decades in fields of energy conversion [1], optical sensing [2], and biomedicine [3], with a unique host of new properties that are notably distinct from their bulk metals [4,5]. The most striking character of these metallic nanoparticles is their unique optical abilities to support localized surface plasmon resonance (LSPR) that arises from the free electrons oscillations [6]. When the metallic nanoparticles are illuminated by light with a wavelength much larger than the nanoparticle size, the high-density free electrons in the metal immediately respond to the external electromagnetic (EM) field and oscillate coherently at the same frequency as the incident light. This plasmon-medicated EM-nanoparticle interaction results in strong light scattering and local EM field enhancement in the vicinity of the nanoparticle surface, which has enabled a variety of applications such as LSPR-enhanced Raman spectroscopy and biological sensing [7][8][9]. As the optical properties of metallic nanoparticles can be fine-tuned by morphology engineering, the LSPR-induced appearance of intense light absorption bands at the tissue's

Results
The LSPR-mediated optical absorption of laser pulses and the resulting thermal responses, such as heat accumulation and thermoelastic PA wave production, are of particular interest for the fields of biological and biomedicine [24]. The light-matter interaction between the incident pulsed laser and GNCs involves a coupled physical process that spans multispatiotemporal scales, where a rigorous theoretical analysis requires a multiscale analysis that encompasses both quantum and continuum-level theory. Here, in order to simplify the problem, we readily assumed that the continuum theory can be rationally used for exploring the PA wave production process. Furthermore, the FEA computational method is employed as a numerical discretization and solution technique to provide Nanomaterials 2020, 10,1942 3 of 13 quantitative simulations of the collective LSPR effect, the temperature increase, and the PA signal production processes. Figure 1 shows the schematic diagram of the nonlinearly enhanced PA/PT effects for nanoassemblies. Nanomaterials 2020, 10, x FOR PEER REVIEW 3 of 13 provide quantitative simulations of the collective LSPR effect, the temperature increase, and the PA signal production processes. Figure 1 shows the schematic diagram of the nonlinearly enhanced PA/PT effects for nanoassemblies. The pronounced optical absorption of the GNCs arises from both the LSPR absorption of both the isolated gold nanosphere and the extremely strong interparticle coupling effects. When the GNCs are irradiated by the incident pulsed laser, free electrons in GNCs quickly jump to high energy levels above the Fermi level by photon-electron interactions within a sub-picosecond [25]. Then, the excited electrons strongly couple between individuals and eventually release their kinetic energy through nonradiative transition by electron-phonon interactions [26]. Assuming that the GNCs are irradiated by a time-harmonic electric field E with the frequency of ω and intensity of E0, the electric field distribution in and around the GNCs can be described using the Helmholtz equation [27,28]: where , , , and are the relative permeability, permittivity, conductivity, and the wavenumber, respectively. The strong interparticle coupling between the isolated nanoparticles in the GNCs can be theoretically described by the coupled dipole equation [29,30]: where Pj represents dipole moments of point dipoles located at positions rj; αj is the dipole polarizabilities; and Ej is the local electric field, which is the difference of the incident field and the retarded induced field by all other electric dipoles. Ajk is a 3 × 3 matrix for ≠ , which can be written as the following [30]: is the 3 × 3 identity matrix. By solving the coupled dipole equations, the absorption cross section of the interparticle coupled GNCs with the total dipole number of N can be evaluated as: The interaction between the electric dipoles is quite complicated, and it is hard to obtain an analytical solution for it. In this study, a FEA-based computational method is employed as a The pronounced optical absorption of the GNCs arises from both the LSPR absorption of both the isolated gold nanosphere and the extremely strong interparticle coupling effects. When the GNCs are irradiated by the incident pulsed laser, free electrons in GNCs quickly jump to high energy levels above the Fermi level by photon-electron interactions within a sub-picosecond [25]. Then, the excited electrons strongly couple between individuals and eventually release their kinetic energy through nonradiative transition by electron-phonon interactions [26]. Assuming that the GNCs are irradiated by a time-harmonic electric field E with the frequency of ω and intensity of E 0 , the electric field distribution in and around the GNCs can be described using the Helmholtz equation [27,28]: where µ r , ε r , σ, and k 0 are the relative permeability, permittivity, conductivity, and the wavenumber, respectively. The strong interparticle coupling between the isolated nanoparticles in the GNCs can be theoretically described by the coupled dipole equation [29,30]: where P j represents dipole moments of point dipoles located at positions r j ; α j is the dipole polarizabilities; and E j is the local electric field, which is the difference of the incident field and the retarded induced field by all other electric dipoles. A jk is a 3 × 3 matrix for j k, which can be written as the following [30]: Nanomaterials 2020, 10, 1942 4 of 13 with k = ω/c, and r jk = r j − r k ,r jk = r j − r k /r jk . M 3 is the 3 × 3 identity matrix. By solving the coupled dipole equations, the absorption cross section of the interparticle coupled GNCs with the total dipole number of N can be evaluated as: The interaction between the electric dipoles is quite complicated, and it is hard to obtain an analytical solution for it. In this study, a FEA-based computational method is employed as a quantitative technique to investigate the LSPR-mediated electric field coupling and the enhanced optical absorption. For gold nanoparticles, the permittivity is frequency-dependent and complex-valued, which can be calculated from an analytical model developed by Rioux et al. [31].
Here, in order to provide a systematic understanding of the resonance coupling between each isolated nanoparticle and their contribution to optical absorption, four basic coupling nanostructure cases, i.e., an isolated nanosphere, a nanosphere dimer, a three-nanosphere chain, and a six-nanosphere chain, are studied. The diameter of the gold nanosphere is set to be 15 nm, and the interparticle distance is set to be 2 nm. The incident laser is polarized parallel to the chain axis indicated by red arrows. As shown in Figure 2a-d, compared with the isolated nanosphere, the interaction of the incident electromagnetic radiation with the coupling nanostructures induces strong polarization of the free electrons and dipolar oscillations of the free electrons. For the cases in Figure 2b,c, when the interparticle distance is much less than the nanoparticle diameter, the near-field interaction arising from evanescent fields is the main coupling mechanism [32]. For the case in Figure 2d, where more nanospheres on a longer GNC are coupled, both the near-field and far-field interactions (arising from the propagating dipolar field) are involved and contribute to the electric field enhancement as well as the optical absorption.
Therefore, the six-nanosphere chain displays the highest LSPR coupling effect and electric field enhancement. The collective resonance of the electric dipoles then converts the EM energy into heat through electron-phonon interactions, where the produced thermal power can be quantitatively calculated in terms of resistive heating [33]: with J as the current density resulting from the free electron oscillations. Here, a plane electric wave with the amplitude of E 0 = 4.6 × 10 4 V/m is induced to irradiate the GNCs and to quantitatively calculate their optical absorption by the FEA method. The calculated incident laser power is about 0.15 mW/µm 2 according to the irradiance expression I in = cε 0 n sur E 2 0 /2, where c, ε 0 , and n sur are the light speed, vacuum permittivity, and refractive index of the surrounding environment, respectively. In our simulations, we set the GNCs to be immersed in water. Figure 2e shows the calculated optical heating power per sphere as a function of the incident laser wavelength for the isolated nanosphere, the nanosphere dimer, the three-nanosphere chain, and the six-nanosphere chain, respectively, where the highest absorption wavelengths are defined as the LSPR wavelengths. Once the absorbed power per sphere is obtained, the laser energy absorbed by the different coupling nanostructures can be calculated by multiplying the laser pulse width. From Figure 2e we can find that, owing to the interparticle resonance coupling, the peak absorption per sphere gradually grows with the chain length, where the per-particle-weighted optical absorption of the six-nanosphere chain is about 3.5 times higher than that of the isolated nanosphere. Meanwhile, as the chain length grows, the LSPR wavelength shows an obvious red shift from about 505 to 545 nm. The results indicate that, by utilizing the interparticle LSPR coupling through constructing GNCs, the per-particle-weighted optical absorption can be greatly enhanced.
polarization of the free electrons and dipolar oscillations of the free electrons. For the cases in Figure  2b,c, when the interparticle distance is much less than the nanoparticle diameter, the near-field interaction arising from evanescent fields is the main coupling mechanism [32]. For the case in Figure  2d, where more nanospheres on a longer GNC are coupled, both the near-field and far-field interactions (arising from the propagating dipolar field) are involved and contribute to the electric field enhancement as well as the optical absorption.  After the GNCs convert the absorbed laser energy into heat, this induces a temperature increase of the GNCs and their surrounding environment owing to the heat diffusion from GNCs to the surroundings, where the temperature field can be described by the heat diffusion equation [34]: where ρ is the density, C p is the heat capacity, T is temperature, k is thermal conductivity, Q abs is the produced thermal power by optical energy, and f (t) is the temporal function of the incident laser pulse.
In this study, we assumed that the laser temporal function can be regarded as a Gauss function with a pulse width of τ = 5 ns and a standard deviation of σ = τ/6, as f (t) = exp −(t − τ/2) 2 /2σ 2 / √ 2πσ. We quantitatively simulated the temperature distribution for the nanosphere dimer at different time points after it was irradiated by a pulsed laser with a laser irradiance of 0.15 mW/µm 2 to simulate the nonlinear PA/PT effects for GNCs. As shown in Figure 2f, when the laser pulse starts to irradiate the GNCs, their temperature gradually increases and then quickly decreases after the laser pulse ends, owing to the thermal diffusion from the GNCs to their surrounding environment. The temperature of the GNCs as a function of time is presented in Figure 2g, where the inset is the laser temporal function.
The local heating of the GNCs by the pulsed laser then induces thermal expansion of the heated area, followed by the production of PA waves, where the produced PA signals can be described by the thermoelastic displacement [35]: which can be usually transformed into p = −ρc L ∂u/∂t. Here, p, c L , and u are the PA pressure, the sound velocity, and the thermoelastic displacement, respectively. Figure 2h shows the simulated thermal expansion and induced stress of the GNCs and their heated surroundings, where obvious thermal-induced volume expansion is observed. By quantitatively extracting the thermal expansion induced PA displacement as a function of time, we obtained the quantitative PA pressure according to Equation (7), shown in Figure 2i. Simulation indicates that under the laser irradiation with an irradiance of 0.15 mW/µm 2 and a laser pulse of 5 ns, the produced in situ PA pressure peak can reach up to~MPa.
To demonstrate the nonlinearly enhanced PT effect of the GNCs, we quantitatively simulated the temperature field at the end of the laser pulse for the isolated nanosphere, nanosphere dimer, three-nanosphere chain, and six-nanosphere chain, respectively, in order to simulate the GNCs with different chain lengths, as shown in Figure 3a-d. Under the same laser irradiation, the six-nanosphere chain GNCs exhibit much higher temperature increase compared with other coupling nanostructures, where two main aspects contribute to this phenomenon. One aspect originates in the relatively large per-particle-weighted optical absorption for the six-nanosphere chain compared with other coupling nanostructures. Meanwhile, owing to the heat diffusion from GNCs to the surroundings arising from the small-size effect, as the chain length grows, the interparticle temperature field overlap contributes significantly to the nonlinear local heating of the GNCs. Figure 3e provides spatial temperature distribution for the coupling nanostructures in Figure 3a-d along the black dotted lines. It can be found from Figure 3e that under a certain laser pulse irradiation, the strong LSPR-coupled GNCs exhibits much higher temperature increase. Furthermore, the corresponding thermal induced volume expansion and stress of the GNCs for the four kind of GNCs are quantitatively simulated as shown in Figure 3f-i. Results show that an obvious thermal expansion is observed for all kinds of GNCs, where the six-nanosphere GNCs exhibit much larger PA displacement and volume expansion compared with other lengths of GNCs. By quantitatively extracting the PA displacement as a function of time, the authors obtained the quantitative PA pressure for the four kinds of GNCs according to Equation (7). As shown in Figure 3j, as the length of the GNCs grows larger, the simulated PA pressure gradually increases, where the PA pressure for the six-nanosphere GNCs is about 3 times higher than that of the individuals.
In order to quantify the nonlinear PA effect for the four kinds of GNCs, we calculated the produced PA wave power for one laser pulse, where the total PA energy of the GNCs can be calculated according to the calculation formula of ultrasound power [36]: Here , p e f f is the effective value of the PA pressure, and S is the surface area of the thermal expansion volume, where we reasonably defined S as the surface area in FEA simulation. The per-particle-weighted PA energy W for one laser pulse is shown in Figure 3k, where the per-particle-weighted PA energy for the six-nanosphere GNCs is amplified about 4.4 times compared with that of the individuals, even under the same laser irradiance. The results indicate that the self-assembling of metallic nanoparticles can nonlinearly enhance the PA/PT effect, owing to the synergy of both optical and thermal coupling effects Nanomaterials 2020, 10, 1942 7 of 13 between the individuals in the GCNs, where with the length of the GNCs growing larger, the nonlinear enhancement increases.
corresponding thermal induced volume expansion and stress of the GNCs for the four kind of GNCs are quantitatively simulated as shown in Figure 3f-i. Results show that an obvious thermal expansion is observed for all kinds of GNCs, where the six-nanosphere GNCs exhibit much larger PA displacement and volume expansion compared with other lengths of GNCs. By quantitatively extracting the PA displacement as a function of time, the authors obtained the quantitative PA pressure for the four kinds of GNCs according to Equation (7). As shown in Figure 3j, as the length of the GNCs grows larger, the simulated PA pressure gradually increases, where the PA pressure for the six-nanosphere GNCs is about 3 times higher than that of the individuals. One of the key factors that influence the EM and thermal coupling among the GNCs is the gap distance between the plasmonic individuals [37]. As indicated by Equations (2)-(4), the EM coupling is strongly sensitive to the relative position of the plasmonic individuals, thus resulting in the intensive dependence for the optical absorption cross sections of the GNCs on the gap distance. As simulated in Figure 4a, by taking the Au dimer as the example for GNCs, with the gap distance d decreasing from 2 nm to 0.5 nm, the collective resonance of the electric dipoles and a strong LSPR-induced electric field enhancement are observed. We quantitatively simulated the optical heating power per sphere under a laser irradiance of 0.15 mW/µm 2 for the Au dimer with the gap distance to be 0.5 nm, 1 nm, and 2 nm, respectively. The results in Figure 4b show that with the gap distance decreasing, the per-particle-weighted peak optical absorption of the dimer grows nearly 47%, with an obvious red shift of the absorption peak. The enhanced optical absorption will certainly induce the increment of the temperature filed of the Au dimer. As quantitatively simulated in Figure 4c, when under the same laser irradiance, with the gap distance decreasing, the temperature of the Au dimer is nonlinearly amplified. Figure 4d shows the spatial temperature distribution for the Au dimer in Figure 4c. The temperature enhancement is about 59%, which is higher than that of the optical absorption enhancement, mainly Nanomaterials 2020, 10, 1942 8 of 13 contributed by the thermal coupling enhancement by decreasing the gap distance. These results indicate that by decreasing the gap distance between plasmonic individuals, the PA/PT effects can be nonlinearly enhanced. The size of the individual is also a critical factor that strongly affects the nonlinear PA/PT effect of the GNCs. As shown in Figure 5a,b, by taking the Au dimer with the gap distance to be 2 nm as example, when irradiated by a laser pulse with an irradiance of 0.15 mW/μm 2 , as the diameter of the individual D increases, the simulated EM field is slightly enhanced, where the hot-spot phenomenon for the larger Au dimer becomes obvious. The quantitatively calculated peak optical heating power per volume of the Au dimer increases with the size. By further quantitatively simulating the temperature field for the Au dimer with a different individual size, we find that even under the same laser irradiance, the temperature field is nonlinearly enhanced with the size increasing, as shown in Figure 5c. Figure 5d is the spatial temperature distribution for the Au dimer in Figure 5c for a different individual size. These results demonstrate that by increasing the individual size, the nonlinear PA/PT effect can be enhanced. While in the practical applications, this size effect should be balanced because it also influences the biocompatibility [38]. The above discussion is established on the premise that the size of the individual in the GNCs is the same, where we can define this case as dimensionally homotypic GNCs. In Figure 5e-h, we provide a demo to quantitatively compare the nonlinear enhancement of the PA/PT effect for dimensionally homotypic and heterotypic GNCs. The heterotypic Au dimer is composed of two Au spheres with a diameter of 30 and 60 nm. The quantitatively simulated EM field distribution and the optical heating power per volume is presented in Figure 5e,f. The result in Figure 5f indicates that compared with the homotypic Au dimer with a diameter of 30 nm, the optical absorption capability of the heterotypic Au dimer becomes weak, which results in the corresponding temperature field decreasing, as shown in Figure 5g,h. The size of the individual is also a critical factor that strongly affects the nonlinear PA/PT effect of the GNCs. As shown in Figure 5a,b, by taking the Au dimer with the gap distance to be 2 nm as example, when irradiated by a laser pulse with an irradiance of 0.15 mW/µm 2 , as the diameter of the individual D increases, the simulated EM field is slightly enhanced, where the hot-spot phenomenon for the larger Au dimer becomes obvious. The quantitatively calculated peak optical heating power per volume of the Au dimer increases with the size. By further quantitatively simulating the temperature field for the Au dimer with a different individual size, we find that even under the same laser irradiance, the temperature field is nonlinearly enhanced with the size increasing, as shown in Figure 5c. Figure 5d is the spatial temperature distribution for the Au dimer in Figure 5c for a different individual size. These results demonstrate that by increasing the individual size, the nonlinear PA/PT effect can be enhanced. While in the practical applications, this size effect should be balanced because it also influences the biocompatibility [38]. The above discussion is established on the premise that the size of the individual in the GNCs is the same, where we can define this case as dimensionally homotypic GNCs. In Figure 5e-h, we provide a demo to quantitatively compare the nonlinear enhancement of the PA/PT effect for dimensionally homotypic and heterotypic GNCs. The heterotypic Au dimer is composed of two Au spheres with a diameter of 30 and 60 nm. The quantitatively simulated EM field distribution and the optical heating power per volume is presented in Figure 5e,f. The result in Figure 5f indicates that compared with the homotypic Au dimer with a diameter of 30 nm, the optical absorption capability of the heterotypic Au dimer becomes weak, which results in the corresponding temperature field decreasing, as shown in Figure 5g,h. 9 of 13 Nanomaterials 2020, 10, x FOR PEER REVIEW 9 of 13 Figure 5. (a) LSPR-induced electric field enhancement for gold nanodimer with different particle sizes; D is the diameter of the plasmonic individuals; (b) Per-sphere-weighted optical heating power with a different particle size; the gap distance is 2 nm; (c) Corresponding quantitative temperature fields for the gold nanodimer with different particle size at t = 5 ns; the laser irradiance is 0.15 mW/μm 2 ; (d) Temperature distribution along the dotted lines in (c). The inset indicates the nonlinear temperature dependence on the size of the individuals; (e) LSPR-induced electric field enhancement for heterotypic gold nanodimer; (f) The per-sphere-weighted optical heating power for heterotypic gold nano-dimer; (g) Corresponding quantitative temperature fields for the heterotypic gold nanodimer; the laser irradiance is 0.15 mW/μm 2 ; (h) The quantitative temperature increment for homotypic and heterotypic gold nano-dimers.
To further demonstrate the nonlinear enhancement of the PA/PT effect for nanoassemblies, the PA performance of GNCs and the individual gold nanospheres are experimentally compared. The samples are brought from XFNANO Tech. Co., Ltd., Nanjing, China, and their concentrations are 2.4 nM. The two samples were firstly both diluted to 2.4 nM. Then, the two diluted samples were characterized by a transmission electron microscope (TEM) to show their morphologies. As shown in the inset of Figure 6a, TEM images shows that the gold nanospheres are uniformly distributed, and the GNCs are long chains, which are quite suitable to be used to validate our purpose. The size of the gold nanospheres and the individuals in the GNCs is about 17 nm. The per-particle-weighted optical absorption spectra for the two samples are shown in Figure 6a. It shows that the massweighted optical absorption of the GNCs is much larger than that of the gold nanospheres, with an obvious red shift of the absorption peak, which quite coincides with the simulations. Then, in order to experimentally demonstrate the nonlinear enhancement of the PA/PT effect for nanoassemblies, Figure 5. (a) LSPR-induced electric field enhancement for gold nanodimer with different particle sizes; D is the diameter of the plasmonic individuals; (b) Per-sphere-weighted optical heating power with a different particle size; the gap distance is 2 nm; (c) Corresponding quantitative temperature fields for the gold nanodimer with different particle size at t = 5 ns; the laser irradiance is 0.15 mW/µm 2 ; (d) Temperature distribution along the dotted lines in (c). The inset indicates the nonlinear temperature dependence on the size of the individuals; (e) LSPR-induced electric field enhancement for heterotypic gold nanodimer; (f) The per-sphere-weighted optical heating power for heterotypic gold nano-dimer; (g) Corresponding quantitative temperature fields for the heterotypic gold nanodimer; the laser irradiance is 0.15 mW/µm 2 ; (h) The quantitative temperature increment for homotypic and heterotypic gold nano-dimers.
To further demonstrate the nonlinear enhancement of the PA/PT effect for nanoassemblies, the PA performance of GNCs and the individual gold nanospheres are experimentally compared. The samples are brought from XFNANO Tech. Co., Ltd., Nanjing, China, and their concentrations are 2.4 nM. The two samples were firstly both diluted to 2.4 nM. Then, the two diluted samples were characterized by a transmission electron microscope (TEM) to show their morphologies. As shown in the inset of Figure 6a, TEM images shows that the gold nanospheres are uniformly distributed, and the GNCs are long chains, which are quite suitable to be used to validate our purpose. The size of the gold nanospheres and the individuals in the GNCs is about 17 nm. The per-particle-weighted optical absorption spectra for the two samples are shown in Figure 6a. It shows that the mass-weighted optical absorption of the GNCs is much larger than that of the gold nanospheres, with an obvious red shift of the absorption peak, which quite coincides with the simulations. Then, in order to experimentally demonstrate the nonlinear enhancement of the PA/PT effect for nanoassemblies, the two samples were prepared in three different concentrations (0.8, 1.2, 2.4 nM) to perform the PA experiments. As shown in Figure 6b, the PA image intensity for the two samples grows with the concentrations, where Figure 6c provides the linear dependence of the PA signal amplitude as a function of concentration for the GNCs. Importantly, even under the same mass concentrations, the GNCs show much higher image intensity compared with that of the gold nanospheres; this is mainly attributed to the nonlinear enhanced PA effect arising from the nanoassemblies. The experimental results coincide well with the simulations. Furthermore, in order to demonstrate the bioapplication of the nonlinearly enhanced PA effect by the nanoassemblies, gold nanospheres and nanochains are injected to in vivo mouse model, and then the PA experiment is immediately performed. The result in Figure 6d shows that the PA imaging intensity for the area of GNCs is much lighter than that of the gold nanospheres, indicating the great potential of the LSPR-mediated nanoassemblies for PA and PT applications.  Figure 6b, the PA image intensity for the two samples grows with the concentrations, where Figure 6c provides the linear dependence of the PA signal amplitude as a function of concentration for the GNCs. Importantly, even under the same mass concentrations, the GNCs show much higher image intensity compared with that of the gold nanospheres; this is mainly attributed to the nonlinear enhanced PA effect arising from the nanoassemblies. The experimental results coincide well with the simulations. Furthermore, in order to demonstrate the bioapplication of the nonlinearly enhanced PA effect by the nanoassemblies, gold nanospheres and nanochains are injected to in vivo mouse model, and then the PA experiment is immediately performed. The result in Figure 6d shows that the PA imaging intensity for the area of GNCs is much lighter than that of the gold nanospheres, indicating the great potential of the LSPR-mediated nanoassemblies for PA and PT applications.

Discussion
The self-assembly of nanoparticles has been of great interest to science and technology, both of which benefit from its ability to provide effective building blocks for physical, chemical, and biological applications. The rational design and construction of the nanoassemblies under theoretical guidance can effectively promote their applications in biological labeling, biochemical sensing, and biomedical imaging/therapy techniques [37]. Moreover, these nanoassemblies are also of tremendous potential to boost the applications in fields of accurate and reliable nanoscale manipulation and micromechanics, such as nanotweezers, microelectromechanical systems, and information storage devices [39].

Discussion
The self-assembly of nanoparticles has been of great interest to science and technology, both of which benefit from its ability to provide effective building blocks for physical, chemical, and biological applications. The rational design and construction of the nanoassemblies under theoretical guidance can effectively promote their applications in biological labeling, biochemical sensing, and biomedical imaging/therapy techniques [37]. Moreover, these nanoassemblies are also of tremendous potential to boost the applications in fields of accurate and reliable nanoscale manipulation and micromechanics, such as nanotweezers, microelectromechanical systems, and information storage devices [39]. This work demonstrates the nonlinearly enhanced PA/PT effect through constructing nanoassemblies by taking GNCs as example, where the strategy is also appliable for other metallic nanoparticles with other shapes [40]. Even though only a one-dimensional optical/thermal coupled nanochains is shown, this strategy can also be readily extended to two-or three-dimensional coupled nanoassemblies with much higher coupling efficiency, and a nonlinear PA/PT effect can be expected. With the number of the nanoparticles in the nanoassemblies increases, the coupling efficiency can become stronger; however, the size of the nanoassemblies should be considered in practical biological applications. Meanwhile, the structural stability of the nanoassemblies should be also considered in practical applications.

Conclusions
In conclusion, we quantitatively simulated the nonlinearly enhanced PA/PT effect by strong LSPR-coupled nanoassemblies. The strong optical and thermal coupling between the individual nanoparticles resulted in a several-time enhancement of the per-particle-weighted optical absorption, consequential thermal field enhancement, and initial PA pressure, compared with their individuals. Quantitative discussion about the dependence of the nonlinear PA/PT enhancement on the assembly chain length, the size of the individuals, the interparticle distance, and the size uniformity of the building blocks were provided. PA experiments on GNCs and individual gold nanospheres were performed to validate the proposition, and experimental results well silhouetted the theoretical discussion. Our work paves the way for the rational design and engineering of plasmonic nanoassemblies with improved PA/PT conversion efficiency.