Amphiphilic compounds including phosphate group such as phospholipids are fundamental components of many creatures. It is well known that these phospholipids are closely related with biological phenomena such as substance transportation and molecular recognition in living bodies  . Monoalkyl phosphate (MAP, Figure 1(a)) is an important candidate for investigation and elucidation of the above phenomena similar to two alkyl chains phospholipid such as dipalmitoyl phosphatidyl choline (DPPC) commonly and commercially known as lecithin. Moreover, the MAP has a characteristic of a dibasic acid, thus it is considered in the perspective of an anionic surfactant      . High yield production of monoalkyl phosphate was done by Kurosaki  and Imokawa et al.  and they conducted safety evaluation to promote and utilize basic detergents for body washer, metal corrosion inhibitor, catalyst for chemical reactions, and antistatic agents in resin treatment.
Among various surfactants, “gemini” surfactant essentially has a structure in which two general or common category monosurfactants having one hydrophobic and one hydrophilic moiety are connected by a spacer group such as an alkyl chain or a peptide group    . In aqueous solution of a gemini surfactant, Krafft point (Kp) and critical micelle concentration (CMC) are lower and the interfacial activity is higher compared to the corresponding monosurfactant. Therefore, gemini surfactants are more superior surfactants, thus several researchers     have investigated them thoroughly. Moreover, gemini surfactants have two hydrophobic alkyl chains, thus they are better candidates to prepare and investigate model membranes for applications in drug delivery systems (DDS)    .
In this research, to increase the scope and applications of gemini surfactants, we have synthesized a novel anionic gemini surfactant, dihexadecyl phosphates (DHGP-n (n = 3, 4, 5, and 6)) as shown in Figure 1(b). DHGP-n has acted as biological substances for model biomembranes. We prepared DHGP-n monolayers on water surface and studied morphology of monolayer using surface tension measurement (STm) and Brewster angle microscopy (BAM). Monolayers were prepared by dropping method  . Structural optimization of DHGP-n molecules was performed by density functional theory (DFT) method and results were discussed at molecular level to justify the morphology of these monolayers.
2. Materials and Methods
Chemicals: 1-Hexadecanol (C16H33OH, 99%), phosphoryl chloride (POCl3, >99%), diethylether (C2H5-O-C2H5, >99%), tetramethyl ammonium hydroxide-5H2O ((N(CH3)4OH-5H2O), >97%)), (methanol, >99%), (Chloroform, >99%) were obtained from Wako Pure Chemical Industries Ltd., Osaka,
Figure 1. Molecular structure of (a): monoalkyl phosphate (MAP) and (b) dihexadecyl Gemini phosphates (DHGP-n, n = 3, 4, 5, 6).
Japan. Purified water with conductance < 0.07 μS/cm was used to prepare monolayer. Water was purified using super water purifying system (WL-21P; Yamato Scientific Corp. Ltd., Tokyo, Japan). Di-bromoalkane (Br-CnH2n-Br, n = 3, 4, 5, and 6, >97%) of mentioned spacer lengths were from Tokyo Chemical Industry Corp., Ltd., Tokyo, Japan.
2.2. Synthesis of Compounds
2.2.1. Synthesis of Gemini Dihexadecyl Phosphates (DHGP-n,) with n = 3, 4, 5, and 6
Novel DHGP-n (n = 3, 4, 5, and 6) were synthesized by modifying synthesis methods of anionic phosphate surfactants  and dialkyl diphosphate  .
2.2.2. Synthesis of 1-Hexadecyl Phosphate ((1-HP, C16H33-O-PO(OH)2)
At first, starting material 1-hexadecyl phosphate (1-HP, Figure 1(a)) was synthesized. At room temperature, 1-hexadecanol (C16H33OH) was reacted to phosphoryl chloride (POCl3) in solvent diethylether (C2H5-O-C2H5) thus product 1-hexadecyl dichloro phosphate (C16H33OPOCl2) was obtained. The product 1-hexadecyl dichloro phosphate then reacted with water in the same solvent to obtain 1-HP. The product was confirmed by nuclear magnetic resonance (NMR) and infrared spectroscopy (IR) recording.
2.2.3. Synthesis of DHGP-n (n = 3, 4, 5, and 6)
Final products of dihexadecyl phosphates gemini surfactants series (DHGP-n,n = 3, 4, 5, and 6) (Figure 1(b)) were synthesized as following. N+(CH3)4)2(C16H33- O-P(O)(O−)2 was obtained when 1-HP was reacted with tetramethyl ammonium hydroxide-5H2O in methanol at room temperature. Di-bromoalkanes (Br-CnH2n-Br, n = 3, 4, 5, and 6) was reacted with N+(CH3)4)2(C16H33-O-P(O)(O−)2 in the same solvent mixture and N+(CH3)4)2(C16H33-O-P(O)(O−)-O-CnH2n-O-P(O)(O−)-O-C16H33 was obtained. Finally, 36% hydrochloric acid was added to the solution to prepare C16H33-O-P(O)(OH)-O-CnH2n-O-P(O)(OH)-O-C16H33 (DHGP-n) series compounds. Products were again confirmed by NMR and IR spectroscopic methods.
2.3. Apparatus Employed and Software Used
STm of monolayer was measured with Surface Tensiometer (CBVP-A3, Kyowa Interface Science Corp. Ltd., Tokyo). Brewster angle microscope used was of the model EMM633K; Filgen Inc., Nagoya, Japan. Image analysis of BAM data was carried out using USB-CAP type (SD-USB2CUP3; AREA Co. Ltd, Tokyo, Japan) imaging analysis software. A CCD camera, C5948-70; Hamamatsu Photonics, Hamamatsu, Japan, was used to detect reflected light from BAM microscope. Molecular area of surfactant (DHGP-n) molecules was calculated for experimental comparison by using a molecular model kit (MOL-TALOU Models, basic-A set, Koriyama, Japan).
2.4. Monolayer Formation
Details of monolayer formation by dropping method have been reported previously  . In brief, 1 mM DHGP-n solution was prepared in chloroform. The solution was dropped and spread on the purified water surface using 100 μl microsyringe (Ge-0583-04; Hamilton Co., Nevada, USA). The dropping was done at 1 μl/min intervals. DHGP-n-chloroform droplet was dropped gently on water surface, next dropping was done only after the previous droplet was well expanded on the surface. The completion of monolayer formation was confirmed by the formation of lens of droplets observed on the water surface and the stable value of surface tension reached as described in the following section.
2.5. Surface Tension Measurement (STm)
STm for each DHGP-n monolayer was measured with a Surface Tensiometer. The droplet of (DHGP-n)-chloroform solution was spread on the purified water surface. Wilhelmy plate dipped in the solution records the value of surface pressure in Surface Tensiometer when droplets were spread on the water surface. The surface pressure of each monolayer was recorded as a function of spreading volume at 26.0˚C ± 0.2˚C.
2.6. Brewster Angle Microscopy (BAM)
The observation of BAM for visualizing the morphology of each DHGP-n monolayer was carried out using a Brewster angle microscope. A p-polarized light of 632.8 nm wavelength was irradiated at the Brewster angle of 53.1˚ on each monolayer surface and the reflected light was detected by a CCD camera. The lateral resolution of BAM was about 1 μm. The formation process of each monolayer was directly observed in situ at arbitrary surface pressures: 2, 12, and 25 mN/m. The BAM measurement was conducted at 26˚C ± 1˚C and measurements were repeated three times for reproducibility.
2.7. Density Functional Theory (DFT) Method Calculation with B3LYP/6-31G(d) Basis Set
DFT calculations by B3LYP/6-31G (d) method have been carried out by using the GAUSSIAN 03 program  for optimized molecular structure of DHGP-n. The initial structural parameters for geometry optimization were assumed to be 60˚ (−60˚) (the G (G’) conformation) torsional angles around two phosphate groups (C16H33-O-(P-O-)3, C16H33-O-(P-O-)3 spacer), all-trans conformation for the other structures (alkyl chain and spacer), and the same direction of two alkyl chains.
3. Results and Discussion
Figure 2 shows surface pressure (SP, π) versus molecular area (A) (π-A) isotherm curves of four DHGP-n (spacer length n = 3, 4, 5, 6) monolayers obtained at 26.0˚C ± 0.2˚C by using dropping method on the water surface. The horizontal axis represents A, calculated from molecular numbers in the dropping volume. Whereas vertical axis represents π plot when equilibrium reached after dropping DHGP-n-chloroform droplet and the SP value becomes constant. Isotherm curve of 1-octadecyl phosphate (1-OP, C18H35-O-PO(OH)2), synthesized by the scheme similar to 1-HP as in Section 2.2.2., has also been shown in Figure 2. Monolayer of 1-OP was obtained by using dropping method for comparison with DHGP-n surfactants. Isotherm of the monolayer formed on the water surface, can be explained in order as shown in the drawing (Figure 2) from the lower right to the upper left.
π-A isotherm curves formed on the monolayer of insoluble amphiphilic molecule at water surface is classified mainly into three categories   . First one shows an extremely gradual increase of surface pressure and followed by sudden and steep increase in π with completion of monolayer formation. The shape of isotherm resembles with a typical titration curve. Such a monolayer is
Figure 2. Surface pressure versus molecular area (π-A) isotherm curve of 1-OP and DHGP-n (n = 3, 4, 5, 16) at 26˚C. ■: DHGP-3, ◆: DHGP-4, ▲: DHGP-5, ▼: DHGP-6, ●: 1-OP.
an indication of existence of single liquid condensed (LC) state. Second one shows a gradual increase in π with monolayer formation. The shape of π-A isotherm is a monotonous curve without any break. Such a monolayer is an indication of existence of single liquid expanded (LE) state. Third one shows a mixed state of the above two changes in a π-A isotherm, at lower π values gradual increase in π is noticed, then there exists a LE-LC transition state with a plateau range in π, and followed by a steep increase in surface pressure, with monolayer formation. The shape of isotherm is a two steps curve.
Compared to the characteristics of the above general π-A isotherm curves, each DHGP-n compound studied in this research showed isotherm curve similar to a LE type monolayer. π increased gradually and monotonously up to 28 mN/m with completion of monolayer formation. The shapes of these isotherm curves were somewhat smooth but did not always resemble with a typical LE state such as dimyristoyl phosphatidyl choline (DMPC) monolayer   . 1- OP had an intermediate isotherm curve between LE and LC states monolayers, a gradual increase in the beginning and subsequently a steep increase in π was noticed with monolayer formation at around 15 mN/m (Figure 2).
LE type of monolayer isotherms of DHGP-n molecules was also confirmed by the observation of BAM. Figure 3 shows BAM images of DHGP-6 monolayer at three πs (recorded in Figure 2): (a) 2 mN/m, (b) 12 mN/m, and (c) 25 mN/m. Figure 3 also shows the image of water surface as comparison. At 2 mN/m (Figure 3(a)), the image is in slight brighter contrast compared to that of water surface. Increasing surface pressure to 12 mN/m (Figure 3(b)), the image became brighter compared to Figure 3(a) and did not have a characteristic of line- and circle-type domain states as observed in the cholesterol and dipalmitoyl
Figure 3. BAM images of DHGP-6 monolayer on various surface pressure at 26˚C. (a) 2 mN/m, (b) 12 mN/m, (c) 25 mN/m, including water surface.
phosphatidyl choline (DPPC) monolayers  . The DHGP-6 monolayer was formed homogeneously. At 25 mN/m pressure (Figure 3(c)), the brightness of image increased further supporting homogeneous formation of DHGP-6 monolayer on water surface. Moreover, small circular-type 3D aggregates appeared in the image around high π (not shown). It was probably because at around the monolayer completion some splatter small lenses (chloroform solution) evaporated on the spot. Furthermore, almost same BAM images were obtained when other DHGP-n (n = 3, 4, 5) monolayers formation were investigated.
In LE state monolayer, it is generally known that the alkyl chain in the monolayer molecule contains either gauche structure or double bond, regardless of the number of alkyl chain in the molecule    . In case of gauche structure, especially the amphiphilic molecule with the carbon number of fourteen or less in alkyl chain forms LE state monolayer at room temperature. Actually, DMPC and 3,4-bis-tetradecylcarbonyl-hexanedioic acid (Gemsurf-14) both having carbon number fourteen formed LE state monolayers at 26˚C, from our recent researches   . Whereas the amphiphilic molecule with the carbon number sixteen or more in alkyl chain such as DPPC or Gemsurf-16 forms LC state monolayer around room temperature    . DHGP-n monolayer of the carbon number sixteen in alkyl chain was different from the above general tendency of molecules because if formed a LE type monolayer (Figure 2). It was also confirmed that two alkyl chains in DHGP-n molecules form almost all-trans structure at the completion of monolayer at π = 25 mN/m as measured by infrared external reflection spectroscopy (IERS)  , CH2 anti-symmetric stretching vibration band (vas(CH2)) for the methylene groups was observed around 2921 cm−1. Spatially free motion of two alkyl chains due to the existence of alkyl chain spacer contributes to the formation of LE type monolayer of DHGP-n molecules. For the same reason, LE state was observed in DPPC monolayer formation in the early stage of monolayer formation  .
A remarkable point noticed in Figure 2, was the limiting molecular area (A0) of DHGP-n monolayers. Under the present experimental conditions, the value of each A0 was 0.91 nm2, 0.71 nm2, 1.00 nm2, and 1.12 nm2 for the compounds DHGP-3, DHGP-4, DHGP-5, and DHGP-6, respectively. Reproducibility of A0 was ±0.05 nm2, experiments were repeated twice for STm measurement. It is expected that the molecular area becomes large with increasing the length of alkyl chain spacer. Approximate calculation done by using molecular model kit described in section 2.3.also showed that the apparent molecular area of hydrophilic group (two phosphates and an alkyl chain spacer) of DHGP-n (Ah-g) was estimated to be 0.40 nm2, 0.47 nm2, 0.54 nm2, and 0.61 nm2 for DHGP-3, DHGP-4, DHGP-5, and DHGP-6, respectively. Whereas A0 of DHGP-n was becoming large by the order of 4, 3, 5, and 6. A0 of DHGP-4 was extremely small among all DHGP-n. A0 of 1-OP was 0.23 nm2 and close to the value of LC monolayer of mono-alkyl amphiphilic molecule such as stearic acid  . Ah-g of 1-OP was also 0.25 nm2 and close to A0. A0 of all DHGP-n molecules was more than three times larger than that of 1-OP.
The specificity of A0 in each DHGP-n monolayer may have close relationship with the conformation of two hydrophilic phosphate groups and hydrophobic interaction between two alkyl chains (-C16H33-, -C16H33-) in DHGP-n. DFT calculations by B3LYP/6-31G(d) method have been carried out for the structural optimization of DHGP-n. These calculations have been performed with consideration that two alkyl chains in DHGP-n molecules orientate toward same direction similar to the molecular structure of monolayer on the water surface   . Figure 4 shows the conformation of each DHGP-n molecule optimized by DFT calculations. The upper picture set (A) is the front of DHGP-n molecule, the middle (B) is the side of molecule, and the lower most (C) is the phosphate groups’ side appearance. In Figure 4, three axes: a, b, and c defined are of DHGP-n molecules. From the initial condition that two alkyl chains orientate toward same direction allowing the interaction between two alkyl chains, the calculation results showed the possibilities: 1) two alkyl chains open outside on the a − b plane and these chain axes are not parallel (Figure 4(a)); 2) the alkyl chains are twisted and open opposite direction reciprocally from a-b plane (Figure 4(b)) and the degree of opening in DHGP-3 and -5 (odd spacers) is larger than DHGP-4 and -6 (even spacers) (Figure 4(b) and Figure 4(c)); 3) a zigzag plane of alkyl chain spacer is perpendicular to the a-b plane and parallel to the b-c plane (Figure 4(a) and Figure 4(c)). Moreover, torsional angles around two phosphate groups were GG (one side) and G’G’ (another side) for DHGP-3 and -5 (odd spacers), and GG (one side) and GG (another side) for DHGP-4 and -6 (even spacers), respectively. Table 1 shows each value (height, width, and depth) of DHGP-n as inserted in DHGP-3 (Figure 4(a) and Figure 4(b)), assuming that the structure of DHGP-n forms a rectangular parallelepiped. Calculated molecular area (Ac-0) was 1.88, 1.28, 1.91, 1.75 nm2 for DHGP- 3, DHGP-4, DHGP-5, and DHGP-6, respectively. The order of Ac-0 and size of
Figure 4. Conformation of DHGP-n (n = 3, 4, 5, 6) molecule optimized by DFT calculations. (a): front of DHGP-n molecule, (b): side of molecule, (c): phosphate groups side.
Table 1. Each value (height, width, and depth) of DHGP-n as inserted in DHGP-3 of Figure 3. These values were measured by assuming that the structure of DHGP-n forms a rectangular parallelepiped.
molecular area did not necessarily agree with the experimental results because the calculation was performed on the condition of one molecule and DHGP-n molecule was in an ideal vacuum. In addition, initial calculation condition was also somewhat specific at high twist state near the hydrophilic group. Nevertheless, DFT calculation showed that only Ac-0 of DHGP-4 was small compared to those of other DHGP-n.
In the case of carbon number of alkyl chain spacer is odd, each oxygen atom in two phosphates connecting to spacer both ends arranges the same side, thereby the degree of twist in two alkyl chains becomes larger for structural stabilization. Even number of spacer, on the other hand, makes the degree of twist smaller. The length of spacer also has an influence on the interaction between two alkyl chains. The longer is the length of spacer, the weaker is the interaction. These multiple effects would show the unique carbon number dependence of alkyl chain spacer on Ac-0. Therefore spacers’ carbon number dependence on Ac-0 is the sum of these effects. The length of alkyl chain spacers in DHGP-n molecules forming monolayer on water surface is responsible for molecule’s structure between all-trans and gauche forms   . When the length of spacer becomes longer, the degree of twist of the spacer would become smaller for forming a monolayer. Due to restricted twist of spacer, the spacer would shrink and thus it strengthens the interaction between two alkyl chains leading to stable monolayer formation. Actually, shrinking of spacer was supported by recorded π-A isotherm, as experimental A0 of each DHGP-n was 30% to 50% smaller than those of calculated Ac-0. This anomaly, which A0 of DHGP-4 surfactant is smallest in the series though DHGP-3 spacer is shortest, has been interpreted considering the balance of the following factors: 1) conformation of each hydrophilic phosphate group restricted between alkyl chain and alkyl chain spacer; 2) length and structure of spacer and its even-odd effect; and 3) interaction between two alkyl chains.
Recent research by the electric conductivity measurement has shown that sodium-DHGP-4 salt (C16H33-O-P(O)(ONa)-O-C4H8-O-P(O)(ONa)-O-C16H33) synthesized from DHGP-4 had critical micellar concentration (CMC) of 0.071 mmol/kg in its water solution at 40˚C. The CMC value was the same order concentration as other Gemini-type surfactant     and also extremely small compared to general single chain anionic surfactants  . With the above result, we would like to expect that this novel DHGP-n would be not only utilized as new detergents and industrial products but also applied as model biomembrane substance and the functional substance such as solubilization surfactant for proteins.
We have investigated the surface properties of the novel DHGP-n (n = 3, 4, 5, 6) compound’s monolayer formed on the water surface at 26˚C, using STm and BAM. It was found that all DHGP-n forms LE type monolayers. The limiting molecular area (A0) was increased in the order n = 4, 3, 5, and 6. Molecular area of compound n = 4 (four methyl spacer) was smallest in the series. DFT calculation also supported the smallest A0 of DHGP-4 among the series. DHGP-n surfactant compounds showed surface properties similar to other commercial and biosurfactants. These first synthesized novel gemini surfactants can be used as commercial detergents and model biomembrane forming substances.
 Brown, D.A., Malkin, T. and Maliphant, G.K. (1955) An X-Ray Examination of Long-Chain Alkyl Dihydrogen Phosphates and Dialkyl Hydrogen Phosphates and Their Sodium Salts. Journal of the Chemical Society, 0, 1584-1588.
 Arakawa, J. and Pethica, B.A. (1980) Micellization in Aqueous Solutions of Monoalkyl Phosphate Salts. Journal of Colloid and Interface Science, 75, 441-450.
 Kurosaki, T., Imokawa, G. and Ishida, A. (1987) The Development of Monoalkyl Phosphate as a Low Skin Irritating Anionic Surfactant. Japan Oil Chemists’ Society, 36, 629-637.
 Imokawa, G. and Tsutsumi, H. (1978) Surface Activity and Cutaneous Effects of Monoalkyl Phosphate Surfactants. Journal of the American Oil Chemist’s Society, 55, 839-843.
 Umezaki, K., Sakai, S., Koeda, S., Yamamoto, Y., Kondo, M., Ikeda, A., Dewa, T., Taga, K., Tanaka, T. and Mizuno, T. (2012) Formation of Planar Bilayer Membranes on Solid Supports Using Peptide Gemini Surfactants. Chemistry Letters, 41, 1430-1432.
 Sommerdijk, N.A.J.M., Hoeks, T.H.L., Synak, M., Feiters, M.C., Nolte, R.J.M. and Zwanenburg, B. (1997) Stereodependent Fusion and Fission of Vesicles: Calcium Binding of Synthetic Gemini Phospholipids Containing Two Phosphate Groups. Journal of American Chemical Society, 119, 4338-4344.
 McGregor, C., Perrin, C., Monck, M., Camilleri, P. and Kirby, A.J. (2001) Rational Approaches to the Design of Cationic Gemini Surfactants for Gene Delivery. Journal of American Chemical Society, 123, 6215-6220.
 Yoshida, D., Yokoyama, T., Shimoaki, T., Tomita, T., Yoshida, T., Yamamoto, Y., Taga, K., Sumino, A., Dewa, T., Nango, M., Yamamoto, M. and Shervani, Z. (2013) Morphology Observation of Dipalmitoyl Phosphatidyl Choline (DPPC) Monolayer on Water Surface by Dropping Method. Journal of Biophysical Chemistry, 4, 115-121.
 Tyagi, P. and Tyagi, R. (2011) Synthesis and Physico-Chemical Properties of Novel Dialkyl Diphosphate Gemini Surfactants Based on Octadecanol. European Journal of Lipid Science and Technology, 113, 848-855.
 Frisch, M.J., Trucks, G.W., Schlegel, H.B., Scuseria, G.E., Robb, M.A., Cheeseman, J.R., Scalmani, G., Barone, V., Mennucci, B., Petersson, G.A, Nakatsuji, H., Caricato, M., Li, X., Hratchian, H.P., Izmaylov, A.F., Bloino, J., Zheng, G., Sonnenberg, J.L., Hada, M., Ehara, M., Toyota, K., Fukuda, R., Hasegawa, J., Ishida, M., Nakajima, T., Honda, Y., Kitao, O., Nakai, H., Vreven, T., Montgomery, J.A., Peralta, J.E., Ogliaro, F., Bearpark, M., Heyd, J.J., Brothers, E., Kudin, K.N., Staroverov, V.N., Kobayashi, R., Normand, J., Raghavachari, K., Rendell, A., Burant, J.C., Iyengar, S.S., Tomasi, J., Cossi, M., Rega, N., Millam, N.J., Klene, M., Knox, J.E., Cross, J.B., Bakken, V., Adamo, C., Jaramillo, J., Gomperts, R., Stratmann, R.E., Yazyev, O., Austin, A.J., Cammi, R., Pomelli, C., Ochterski, J.W., Martin, R.L., Morokuma, K., Zakrzewski, V.G., Voth, G.A., Salvador, P., Dannenberg, J.J., Dapprich, S., Daniels, A.D., Farkas, O., Foresman, J.B., Ortiz, J.V., Cioslowski, J. and Fox, D.J. (2009) Gaussian 09, Revision D.01. Gaussian, Inc., Wallingford.
 Yamamoto, Y., Yokoyama, T., Yoshida, D., Mori, H., Sekiguchi, K., Shimoaki, T., Yoshino, A., Taga, K., Shervani, Z. and Yamamoto, M. (2015) Interactions between Phospholipid Monolayers (DPPC and DMPC) and Anesthetic Isoflurane Observed by Quartz Crystal Oscillator. Journal of Biophysical Chemistry, 6, 42-53.
 Yoshida, T., Yamamoto, Y., Taga, K., Kamaya, H. and Ueda, I. (2003) Effects of Water-Soluble Alcohols on the Surface Conductance of Lipid Monolayers. Journal of Physical Chemistry B, 107, 3196-3198.
 Kobayashi, Y., Hirai, R., Ito, D., Yamamoto, Y., Taga, K. and Shervani, Z. (2016) Morphology Observation of Novel 3,4-Bis-Alkyloxycarbonyl-Hexanedioic Acid Monolayer on Water Surface by Dropping Method. Journal of Biophysical Chemistry, 7, 43-48.
 Yamamoto, M., Suzuki, M., Kimura, T. and Itoh, K. (2008) Molecular Structure at Free Surfaces of Liquid n-Heptadecane Investigated by Infrared External Reflection Spectroscopy. Journal of Physical Chemistry B, 112, 13232-13239.
 Kato, R., Nakahara, H. and Shibata, O. (2015) Effect of Spacer Length I Fluorinated Amphiphilies on DPPC Monolayers. Proceedings of 66th Division of Meeting Colloids and Surface Chemistry, P-192.