The molecular mechanisms of plasticity in crystal forms of theophylline

The molecular mechanisms of plasticity in crystal forms of theophylline

PDF

Hongji Liu^a, H.C. Stephen Chan^b^,^c^,^*, Li Zhang^d^,^**, Yang Lu^d, Jiaquan Li^a, Jianting Li^a, Liang Li^e, Zhengzheng Zhou^a^,^**

Chinese Chemical Letters | 2023, 34(8) : 108057

Less

Chinese Chemical Letters | 2023, 34(8): 108057

• Communication •

The molecular mechanisms of plasticity in crystal forms of theophylline

Full

Hongji Liu^a, H.C. Stephen Chan^b^,^c^,^*, Li Zhang^d^,^**, Yang Lu^d, Jiaquan Li^a, Jianting Li^a, Liang Li^e, Zhengzheng Zhou^a^,^**

Affiliations

^a NMPA Key Laboratory for Safety Evaluation of Cosmetics, Guangdong Provincial Key Laboratory of Tropical Disease Research, Department of Hygiene Inspection & Quarantine Science, School of Public Health, Southern Medical University, Guangzhou 510515, China

^b Actmol AG, Hagentalerstrasse 49, Basel CH-4055, Switzerland

^c Research Center for Computer-Aided Drug Discovery, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China

^d Beijing City Key Laboratory of Polymorphic Drugs, Institute of Materia Medica, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing 100050, China

^e Department of Forensic Toxicological Analysis, Zhongshan School of Medicine, Sun Yat-sen University, Guangzhou 510080, China

Published: 2023-08-15 doi: 10.1016/j.cclet.2022.108057

Outline

Abstract

Less

Plastic and elastic behaviors of organic crystals have profound influence on the processability of pharmaceutical substances. Analogous to metals, the identifications of molecular slip planes in organic crystals are regarded as a strategy for harnessing plasticity. In this work, we experimentally characterized the form Ⅱ anhydrous theophylline (THPa) and its monohydrate (THPm) for their distinct plastic and elastic behaviors. Extensive DFT calculations were performed to model the effects of increasing lattice strains on molecular packing. We discovered that the energy barrier associated with the strain-induced molecular rearrangement would link to the plasticity of THPa, and possibly other simple aromatic compounds. Meanwhile, water molecules in THPm disrupt the stacking architecture from THPm and effectively undermine the general mechanism for plasticity. Hydrate formation would therefore be an alternative strategy to engineer the mechanical property of organic crystalline materials.

Key words

Theophylline / Elasticity / Plastic crystal / Density functional theory / Mechanochemistry

Cite this Article

Hongji Liu, H.C. Stephen Chan, Li Zhang, Yang Lu, Jiaquan Li, Jianting Li, Liang Li, Zhengzheng Zhou. The molecular mechanisms of plasticity in crystal forms of theophylline[J]. Chinese Chemical Letters, 2023 , 34 (8) : 108057 - . DOI: 10.1016/j.cclet.2022.108057

Full Text

Less

In recent years, the elastic and plastic deformations of organic crystals are brought into the spotlight of both material science and crystal engineering communities [1-5]. Elastic crystal is thought to require weak isotropic interactions, such as π–π interaction and many others [3,6,7], but more importantly, an interlocking packing architecture that hinders the slippage of molecular layers as lacked in plastic crystals [3,7-9]. Hirshfeld surface calculation is often applied to locate important molecular contacts for maintaining the packing architecture [10,11]. In pharmaceutical industry, plastic crystals with lower stiffness may be pursued for better powder compaction and tabletability [12-14]. Crystal slip plane is a characteristic feature observed in plastic crystal [7,15,16]. Energy framework, parameterized from quantum mechanical calculations, is capable of identifying potential slip planes within a crystal, when the interaction within a certain molecular layer significantly outweighs that between the layers [17-20]. However, the underlying mechanism that restricts these layers from restoring their original positions has not been fully addressed. In this work, we focus on a naturally occurring caffeine derivative, theophylline (THP, Scheme 1), and discovered the plasticity behavior in its anhydrous crystal form Ⅱ. The molecular layers formed via π-π stacking become excellent targets for theoretical investigations. During our characterization experiments, the monohydrate of THP presents excellent elastic property, where π-π stacking also exists. The role of water in tuning the mechanical property therefore becomes another intriguing subject for investigation. DFT level calculations were performed to predict the molecular events in these two crystal forms, in response to various anisotropic deformations [21-24]. The resulting energy landscapes may shed light on their disparate mechanical behaviors.

Our experimental data demonstrate the distinct plastic and elastic behaviors in anhydrous theophylline (THPa) and the monohydrate (THPm) crystals respectively. THPa crystal exhibits super plasticity and undergoes irreversible bending on (200) face (Fig. 1a). By applying the load with a metallic needle on this face, three independent THPa crystals would be molded in the letters “SMU” respectively (Fig. 1b). Meanwhile, THPm crystal undergoes elastic deformation and shows excellent bending resilience when the load is repeatedly applied on its (002) face (Fig. 1c). The maximum tolerable elastic strain was measured to be approximately 1.22% for a 20 µm thick THPm crystal (Fig. S4 in Supporting information), lower than the exceptionally elastic celecoxib crystal (3.56%) [25] but 2–3 times higher than most crystalline alloys (0.50%) [26,27]. Moreover, the fractured THPm crystallites formed beyond its elastic limit do not remain in an irreversible plastic deformation, but restore their linear shapes (Fig. S3 in Supporting information).

Nanoindentation [17,28] coupled with atomic force microscopy (AFM) was used to investigate further the differences in mechanical behaviors between THPa and THPm (Fig. 2a). The load-displacement (P-h) curve of THPa shows a stronger linear relationship and a consistently higher resistance to deformation. In contrast, the P-h curve of THPm significantly deviates from linearity during loading and unloading. THPm is highly compliant to deformation at a small load and develops resistance gradually with the increasing penetration depth. Another key difference is that discrete “pop-in” events are observed in the P-h curve of THPa, but not of THPm. The differences in penetration depth between any two subsequent events are similar, suggesting that the plastic transformation in THPa could be related to certain periodically stacked layers. The images from the scanning electron microscope (SEM) [1,29,30] confirm that the THPa crystals indeed consist of a multilayer topology (Fig. S5a in Supporting information). Meanwhile, surface pile-up phenomenon is observed from the AFM results of THPa, indicating that the molecular slippage direction would be perpendicular to the applied stress [17,28,31]. More importantly, the pile-up occurs predominantly on one particular edge of the indenter, suggesting that the plastic transformation of THPa may also be highly directional on the surface (Fig. 2b). For THPm, the absence of “pop-in” event on P-h curve implies an elastic deformation throughout the whole course of the more profound indentation. At a peak load of 1000 µN, the maximum penetration of the probe in THPm (1318.86 nm) is 4.22 times deeper than that in THPa (312.36 nm). The AFM results of THPm show no pile-up of material, indicating that displacement of molecules proceeds along the same direction of the applied stress [28]. More interestingly, the surface of THPm remains smooth with no sign of previous indentation, after the stress has been released (Fig. 2b). Therefore, THPm demonstrates a more superior ability in recovering its crystal shape.

X-ray crystallography reveals that the inclusion of water molecule significantly alters the packing architecture (Table 1 and Fig. S6 in Supporting information). THPa crystallizes in an orthorhombic space group Pna2₁ with only one THP molecule in the asymmetric unit. The THP molecules form a 1-D hydrogen bond chain along the [011] direction. Via the screw symmetry, these THP molecules result in wavy molecular sheets that stack on one another by π-π interaction along the b-axis. (Fig. 3a). Meanwhile, THPm crystallizes in a monoclinic space group P112₁/b instead, with one THP and one water molecule in the asymmetric unit. Thermal analysis also supports equimolar of water and THP in the lattice (Fig. S7 in Supporting information) [32]. Water molecules presumably form a 1-D hydrogen bond chain along the a-axis. THP molecules no longer adopt a 1-D hydrogen bond chain but 0-D dimers, which alternatingly interact with water along the c-axis and result in a 3-D hydrogen bond network. Noticeably, these THP dimers also stack on one another by π–π interaction, and the putative slip plane appears along [120] direction (Fig. 3b). More plausible slip planes can be identified by further visual inspection, including the (200) plane in THPa at the interface of two 1-D hydrogen bond chains of THP and the (002) plane in THPm along the water channel (Figs. 3c and d). One may conceptualize plastic deformation as a consequence of applied strain along a slip plane, and elastic deformation as a process during which a material undergoes simultaneously compression and dilation along the direction of applied strain (Figs. 3e and f) [2]. However, how exactly water harnesses elasticity in THPm and the determining factor for plasticity of THPa remain unclear. Atomistic modeling was therefore performed to elucidate the underlying reasons.

When subject to an anisotropic strain, atoms inside the deformed lattice are displaced and attain a new configuration that corresponds to a local energy minimum. By probing the structural changes and the relating potential energy surface, the comprehensive interplay between macroscopic mechanical properties and the corresponding atomistic events could be better understood. For a given crystal lattice R matrix, its deformed lattice R' can be expressed as a transformation via a distortion matrix D [33-35]:

(1)

where δ_i corresponds to an anisotropic strain. The internal energy of a crystal E under a set of general strain δ_i can be expanded in form of a Taylor's series [36]:

(2)

where E₀ and V are the ground state energy and the volume of the unstrained crystal respectively, σ_i is the element of stress tensor with respect to δ_i, and c_ij is the anisotropic elastic modulus. c₁₁, c₂₂ and c₃₃ account for the axial compression (or dilation) along the a-, b- and c-axis respectively, whereas c₄₄, c₅₅ and c₆₆ account for the shear across the bc-, ca-, and ab-plane respectively [33-35,37]. Note that some shear directions may break the symmetry of the original space group [34]. Other remaining elastic moduli are symmetrical, i.e., c_ij = c_ji, and correspond to a combination of any two of the abovementioned deformations [37]. At small strains, the internal energy varies almost like a harmonic function, and each c_ij can be obtained from the second-order coefficient of the Taylor's series with respect to the deforming strain component (Table 1). Regarding the ground state energy of THPm, the water protons are under the inversion symmetry constraint of the P112₁/b space group and result in a suboptimal “head-to-head” contact with one another after DFT-XDM optimization (Fig. 4a) [38-42]. The removal of inversion constraint effectively lowers the space group to P11b and the water protons instead form a “zig-zag” hydrogen bond chain with neighboring water molecules (Fig. 4b). This alternative proton configuration yields an enthalpy excess of nearly 2.4 kcal/mol and inflicts a minimal structural change with a RMSD₃₀ of 0.176 Å only (Table 1). Therefore, the P11b structure was taken as the unstrained state of THPm for calculating the anisotropic elastic moduli.

The calculated anisotropic stiffness predominantly depends on the impact of deformation direction on hydrogen bonding. The high axial stiffness in THPm, for example, is attributed to the distortion of the water hydrogen bond chain (for c₁₁), the THP dimer (for c₂₂) or both (for c₃₃) caused by each axial strain (Table 1). Similarly, in THPa, the strains along the a- or c-axis, particularly the former, directly disrupt the hydrogen bond distance and angle, whereas the strain along the b-axis does not significantly alter the relative position of THPs involved in the hydrogen bond. Thus, c₁₁ and c₃₃ are nearly 2 times larger than c₂₂ (Table 1). Nevertheless, THPm in general has larger axial moduli than THPa does. Specifically, c₂₂ and c₃₃ of THPm are around 1.5–2.0 times greater than the largest axial modulus c₁₁ of THPa. The 3-D hydrogen bond network in THPm offers a more robust architecture against axial strains. This prediction result appears paradoxical to the nanoindentation experiment, where THPm shows good compliance to an axial stress. However, its shear moduli are noticeably smaller than those of THPa (Table 1). Its two smallest moduli, c₄₄ and c₅₅, are even less than 50% of the second largest shear modulus c₅₅ of THPa, because these shear directions preserve the quasi-inversion relationship between the THP dimer molecules and inflict minimal effect to the hydrogen bonding (Fig. 5a). Consequently, the generally low shear stiffness of THPm may be the major reason for its high compliance in the nanoindentation experiment, provided that the shear direction would be parallel to the applied stress. It is interesting to note that, when modeling the stiffest shear modulus c₆₆ of THPm, the increasing strain changes the tilting angle of the THP stack, which then squashes against the neighboring water chain and eventually dislodges it slightly along the c-axis. This structural event is reflected from the discontinuity of the energy-strain curve (Fig. 5b). The results from our Hirshfeld surface and energy framework calculations indicate a potential slip plane on the (002) surface here (Figs. S14 and S15b in Supporting information). In the highly strained lattices, the displacement of the water chain occurs exactly on this predicted slip plane, even though no slipping of the THP layer is observed.

Abrupt structural change induced by an applied strain was also predicted in THPa. At a moderate axial strain on a-axis, each strained 1-D hydrogen bond chain is restored to normal at the expense of the weaker dispersion between the chains. This molecular event is also reflected from the discontinuity on the energy-strain curve. Eventually isolated THP chains are observed in the highly strained structure (Fig. 6a). Our Hirshfeld surface and energy framework calculations also indicate this (200) surface as a potential slip plane (Figs. S13 and S15a in Supporting information). Likewise, the simulated shear across the ab-plane unexpectedly results in the largest modulus in THPa, despite its presumably limited impact on hydrogen bonding. One explanation would be the wavy layout of molecular sheet largely hinders its translocation in this particular direction. More strikingly, our simulation shows that, with an increasing strain (δ₆/2), the original wavy architecture (up-down-down-up) undergoes rearrangements twice: (1) Beyond 2.5% strain, the architecture gradually turns into an up-up-up-up layout; (2) At 7.0% strain, there is a slight slippage on the (020) plane and the THP in the same hydrogen bond chain becomes almost coplanar (Fig. 6b). Each rearrangement is associated with a huge energy barrier, followed by a gradual convergence to a local energy minimum. This energy landscape may well explain the plastic behavior of THPa: Firstly, each metastable strained state is now restricted from returning to its preceding state, and eventually to the unstrained equilibrium state. Secondly, the multiple “pop-in” events during nanoindentation could be related to the stepwise plastic transformations induced by this shear strain. Thirdly, the highly directional pile-up in AFM may be caused by this directional plastic transformation. Finally, the average intermolecular HO distance between the C—H on the heterocylic ring and the O═C group decreases with an increasing strain (δ₆/2), i.e.2.21 Å at 0.0%, 2.18 Å at 6.5% and 2.16 Å at 10.0%. This could have weakened the interactions between these two functional groups. Micro-Raman spectroscopy (Fig. S8 in Supporting information, top) indeed shows a stronger blue shift of the C—H stretch in the inner arc of the bent THPa crystal (2966.19 cm^–1), when compared to the outer arc (2964.53 cm^–1). This experimental finding may further support this shear direction would be highly related to the plasticity of THPa.

To investigate whether the irreversible deformation in plastic crystal is linked with a significant energy barrier during structural rearrangement, similar calculations were extended to two more well-known examples of plastic crystals. Caffeine monohydrate is noteworthy for its plastic behavior [12], but shares similar packing features with theophylline monohydrate, including space group, 1-D water chain and the xanthine dimer. In our extended calculations, shearing across each of the three faces in caffeine monohydrate crystal does not result in a significant discontinuity on the energy-strain curve as well as rearrangements of the caffeine stacking. However, a compression on the b-axis alongside a dilation on the a-axis does (See supporting information). The weak hydrogen bond between caffeine dimer, instead of the strong hydrogen bond between theophylline dimer, offers more flexibility for rearranging the caffeine stacks. Similarly, in the P2₁/c setting of hexachlorobenzene [45], another notable example of plastic crystal, our extended calculations also show that molecules overcome an energy barrier and adopt a different tilting angle in response to a simultaneous dilation on the c-axis and an increasing β angle (see Supporting information). Our results possibly suggest the energy barrier hampering the strain-induced π-stacking rearrangement could be a general mechanism for plasticity in crystals of small aromatic molecules.

In this study, water molecule abolishes completely the plastic property of anhydrous theophylline and converts it into a flexible material with excellent resilience against deformation. Under a small applied strain, our predicted anisotropic stiffness of the materials shows strong connections to the directional hydrogen bond network inside the crystal lattice. Axial strain along a 1-D hydrogen bond chain or a 0-D dimer faces huge resistance, whereas shearing about the latter tends to experience less resistance. When a higher strain is further applied, molecular slip planes or related structural events predicted from Hirshfeld surface analysis and energy framework calculations are also captured in our anisotropic stiffness predictions, reflected by a discontinuity on the energy-strain curve. This outcome indicates that all these methods can be complimentary to one another. Both our experimental and predicted results suggest aromatic stacking and molecular slip plane are prerequisites but may not be the only factors for the plasticity of an organic crystal. Additional mechanisms may be required for hindering the recovery of the material back to its unstrained state – here the energy penalty for rearranging the wavy molecular sheets in anhydrous theophylline is a good example. Meanwhile, in the monohydrate form, water molecules form a 1-D chain, break down the wavy sheets of theophylline into isolated columns of theophylline dimers, and offers a robust 3-D hydrogen bonded architecture to undermine the molecular slippage for plasticity. Moreover, the internal energy surface remains as a continuous concave function for a wide range of strain, indicating that a deformed THPm would be capable of returning to its unstrained equilibrium state. Since these calculations investigate the impact of more extreme strains on the molecular packing, there may be a potential application for tracking the pathway of mechanically induced phase transformation [46]. To conclude, searching for a stable crystal form with optimal pharmaceutical properties is the ultimate objective in experiment. The ability to predict further the material properties would definitely add values to crystal structure prediction [47-49].

Declaration of competing interest

Less

The authors declare no competing financial interest.

Acknowledgments

Less

This study was supported by the Young Scientists Promotion Fund of Natural Science Foundation of Guangdong Province, National Natural Science Foundation of China (No. 81703438), the Key R&D Program of Shandong Province (No. 2021ZDSYS26) National Science and Technology Major Project of China (No. 2018ZX09711001-001-013), CAMS Innovation Fund for Medical Sciences (No. 2017-I2M-3-010) and Student Innovation Training Program (No. 201812121093). Authors would like to thank Dr Grahame R. Woollam (Novartis Pharma AG, Switzerland), Dr Hafsa S. Javed (University of Malta) and Dr Sten Nilsson (AstraZeneca, Sweden) for their valuable discussions.

Supplementary materials

Less

Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.cclet.2022.108057.

References

Less

[1]

E. Ahmed, D.P. Karothu, M. Warren, P. Naumov, Nat. Commun. 10 (2019) 3723.

[2]

S. Ghosh, M.K. Mishra, S.B. Kadambi, et al., Angew. Chem. Int. Ed. 54 (2015) 2674–2678.

[3]

S. Saha, M.K. Mishra, C.M. Reddy, G.R. Desiraju, Acc. Chem. Res. 51 (2018) 2957–2967.

[4]

B. Tang, B. Liu, H. Liu, H. Zhang, Adv. Funct. Mater. 30 (2020) 2004116.

[5]

J.R. Wang, M. Li, Q. Yu, et al., Chem. Commun. 55 (2019) 8532–8535.

[6]

H. Liu, Z. Lu, B. Tang, et al., Angew. Chem. Int. Ed. 59 (2020) 12944–12950.

[7]

J. Li, J. Li, H. Liu, et al., Chin. Chem. Lett. 33 (2022) 4069–4073.

[8]

E. Ahmed, D.P. Karothu, P. Naumov, Angew. Chem. Int. Ed. 57 (2018) 8837–8846.

[9]

S. Saha, G.R. Desiraju, J. Am. Chem. Soc. 139 (2017) 1975–1983.

[10]

S.P. Thomas, M.W. Shi, G.A. Koutsantonis, et al., Angew. Chem. Int. Ed. 56 (2017) 8468–8472.

[11]

P.R. Spackman, M.J. Turner, J.J. McKinnon, et al., J. Appl. Crystallogr. 54 (2021) 1006–1011.

[12]

S. Hu, M.K. Mishra, C.C. Sun, Chem. Mater. 31 (2019) 3818–3822.

[13]

Z. Zhou, W. Li, W.J. Sun, et al., Int. J. Pharm. 509 (2016) 391–399.

[14]

D. Guan, B. Xuan, C. Wang, et al., Pharmaceutics, 13 (2021) 2160.

[15]

L. Dezerald, D. Rodney, E. Clouet, et al., Nat. Commun. 7 (2016) 11695.

[16]

C. Wang, C.C. Sun. Mol. Pharm. 16 (2019) 1732–1741.

[17]

K. Zhang, C.C. Sun, Y. Liu, et al., Chem. Mater. 33 (2021) 1053–1060.

[18]

C.F. Mackenzie, P.R. Spackman, D. Jayatilaka, M.A. Spackman, IUCrJ 4 (2017) 575–587.

[19]

M.J. Turner, S.P. Thomas, M.W. Shi, et al., Chem. Commun. 51 (2015) 3735–3738.

[20]

C. Wei, J. Wang, Y. Zhang, et al., Chin. Chem. Lett. (2022) 107896.

[21]

H. Liu, J. Nie, H.C. Stephen Chan, et al., Int. J. Pharm. 601 (2021) 120537.

[22]

D.D. Huang, H.C.S. Chan, Y.S. Wu, et al., J. Mol. Liq. 329 (2021) 115604.

[23]

H.M. Yu, B.X. Zhang, W.H. Xing, et al., Chin. Chem. Lett. 34 (2023) 107668.

[24]

H.C. Stephen Chan, J. Kendrick, M.A. Neumann, F.J.J. Leusen, CrystEngComm 15 (2013) 3799–3807.

[25]

K. Wang, M.K. Mishra, C.C. Sun, Chem. Mater. 31 (2019) 1794–1799.

[26]

M. Telford, Mater. Today 7 (2004) 36–43.

[27]

S. Ghosh, M.K. Mishra, S.B. Kadambi, et al., Angew. Chem. Int. Ed. 54 (2015) 2674–2678.

[28]

S. Varughese, M.S. Kiran, U. Ramamurty, G.R. Desiraju, Angew. Chem. Int. Ed. 52 (2013) 2701–2712.

[29]

Y. Chen, Z. Chang, J. Zhang, J. Gong, Angew. Chem. Int. Ed. 60 (2021) 22424–22431.

[30]

P. Commins, A.B. Dippenaar, L. Li, et al., Chem. Sci. 12 (2021) 6188–6193.

[31]

P. Sanphui, M.K. Mishra, U. Ramamurty, G.R. Desiraju, Mol. Pharm. 12 (2015) 889–897.

[32]

H. Liu, H.H.Y. Tong, Z. Zhou, J. Therm. Anal. Calorim. 147 (2022) 12947–12963.

[33]

J.H. Jang, S.J. Park, T.H. Lee, H.K.D.H. Bhadeshia, Mater. Sci. Tech. 36 (2020) 615–622.

[34]

K.B. Panda, K.S.R. Chandran, Acta. Mater. 54 (2006) 1641–1657.

[35]

G. Ghosh, Aip. Adv. 5 (2015), 087102.

[36]

D.C. Wallace, Am. J. Phys. 40 (1972) 1718–1719.

[37]

A.U. Ortiz, A. Boutin, A.H. Fuchs, F.X. Coudert, Phys. Rev. Lett. 109 (2012) 195502.

[38]

P. Giannozzi, S. Baroni, N. Bonini, et al., J. Phys. Condens. Matter 21 (2009) 395502.

[39]

P. Giannozzi, O. Andreussi, T. Brumme, et al., J. Phys. Condens. Matter 29 (2017) 465901.

[40]

J.P. Perdew, K. Burke, M. Ernzerhof, Phys. Rev. Lett. 77 (1996) 3865–3868.

[41]

A. Dal Corso, Comp. Mater. Sci. 95 (2014) 337–350.

[42]

A. Otero-de-la-Roza, E.R. Johnson, J. Chem. Phys. 136 (2012) 174109.

[43]

J.A. Chisholm, S. Motherwell, J. Appl. Crystallogr. 38 (2005) 228–231.

[44]

C.F. Macrae, I. Sovago, S.J. Cottrell, et al. J. Appl. Crystallogr. 53 (2020) 226–235.

[45]

M.K. Panda, S. Ghosh, N. Yasuda, et al., Nat. Chem. 7 (2015) 65–72.

[46]

D. Hasa, W. Jones, Adv. Drug. Deliver. Rev. 117 (2017) 147–161.

[47]

D.D. Huang, D.Z. Yang, Y. Lv, et al., J. Mol. Struct. 1250 (2022) 131891.

[48]

J. Du, J. Guo, D. Kang, et al., Chin. Chem. Lett. 31 (2020) 1695–1708.

[49]

A. Halme, M.J. Quayle, S.O. Nilsson Lill, et al., Cryst. Growth. Des. 19 (2019) 3670–3680.

Appendix

Less

Year 2023 volume 34 Issue 8

PDF

Cite this Article

BibTeX

Article Info

doi: 10.1016/j.cclet.2022.108057

Receive Date：2022-11-06
Online Date：2025-11-21
Published：2023-08-15

Article Data

Affiliations

History

Received：2022-11-06
Revised：2022-12-04
Accepted：2022-12-09

Affiliations

^b Actmol AG, Hagentalerstrasse 49, Basel CH-4055, Switzerland

^c Research Center for Computer-Aided Drug Discovery, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China

^d Beijing City Key Laboratory of Polymorphic Drugs, Institute of Materia Medica, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing 100050, China

^e Department of Forensic Toxicological Analysis, Zhongshan School of Medicine, Sun Yat-sen University, Guangzhou 510080, China

References

Share

https://castjournals.cast.org.cn/joweb/ccl/EN/10.1016/j.cclet.2022.108057

Share to

Scan QR to access full text

Cite this article

BibTeX

Citations

表12种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

Scheme 1 The molecular structures of caffeine and theophylline (THP).

Fig. 1. (a) THPa crystal shows plastic deformation, when applying a load with forceps and metallic needle. (b) THPa crystals with super plasticity were able to be molded into letters “SMU” respectively. (c) THPm crystal exhibits elastic deformation, where the bent crystal restores its linear shape when the applied force is withdrawn.

Fig. 2. (a) Representative load–displacement (P-h) curves obtained from nanoindentation on THPa and THPm crystals. “Pop-in” events (three black arrows) happened when loading THPa. (b) Atomic force microscopy images for the nanoindentation on THPa with pile-up predominantly located on the upper side of the indent. (c) In contrast, a relatively smooth surface after nanoindentation on THPm crystals.

Fig. 3. (a) The wavy stacking of THPa along the b-axis. (b) Stacking of THP dimers in THPm along [120] direction. Other putative slip planes (c) (020) in THPa and (d) (002) THPm; Possible molecular arrangement during plastic or elastic deformation in (e) THPa and (f) THPm crystals.

Fig. 4. (a) The inversion symmetry in P112₁/b space group constraints the water proton positions and interrupts hydrogen bonding between neighboring water molecules; (b) Without inversion symmetry in P11b space group, water protons are free to attain an optimal hydrogen bond geometry.

Fig. 5. (a) Shear strain δ₄ in THPm has less significant effect on the hydrogen bond of THP dimers. The hydrogen bonds in the unstrained and the sheared structure (δ₄/2 = 10.0%) are colored in cyan and yellow respectively. (b) Energy-strain curve of THPm under shear strain δ₆ which eventually causes a displacement of water chain in the c-axis direction (highlighted in yellow circles) and a discontinuity on the curve (zones Ⅰ and Ⅱ).

Fig. 6. (a) Energy-strain curve of THPa with axial strain δ₁. Discontinuity is due to restoration of the 1-D hydrgeon bond chain (zone Ⅰ). Isolated chains are seen in highly strained structure (zone Ⅱ), and (b) with shear strain δ₆. The wavy molecular sheet undergoes significant rearrangements and reaches different metastable strained states (zones Ⅱ and Ⅲ respectively).

Articles: Latest Articles; Most Read; Collections

Updates: Events; News; Multimedia

About: About Us

Contact

No. 86 Xueyuan South Road, Haidian District, Beijing

100081

010-62199257

qkjq@cast.org.cn

Copyright © 2025 China Association for Science and Technology. All rights reserved. For all open access content, the relevant licensing terms apply.
Sponsored by the Office of the Leading Group for Cybersecurity and Informatization of CAST, and supported by Science and Technology Review Publishing House