• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Preparation of reduced sensitivity co-crystals of cyclic nitramines using spray flash evaporation

    2020-04-09 18:41:52MrinlGhoshSikderShiblBnerjeeTlwrSikder
    Defence Technology 2020年1期

    Mrinl Ghosh ,A.K.Sikder *,Shibl Bnerjee ,M.B.Tlwr ,N.Sikder

    a High Energy Materials Research Laboratory,Pune 411 021,India

    b Defence Institute of Advanced Technology,Pune 411 025,India

    Keywords:Cyclic nitramine Solubility Spray evaporation Co-crystal Reduced sensitivity

    ABSTRACT The present day weapon technology demands novel energetic materials that exhibit simultaneous high explosive yield and reduced sensitivity.This article demonstrates application of spray evaporation to prepare reduced sensitive co-crystals of high performance nitramine explosives like HMX and CL-20 with a relatively less insensitive explosive 1,1-diamino-2,2-dinitroethylene or FOX-7.Stronger intermolecular hydrogen bonding in FOX-7 is responsible for limited solubility in most of organic solvents.Large solubility differences of FOX-7 with HMX and CL-20 restricts it's co-crystallization through classical methods that yields thermodynamically favorable product.Spray flash evaporation,a kinetic crystallization method,has been therefore adopted and could successfully produce CL-20/FOX-7(2:1)and HMX/FOX-7(4:1)co-crystals.The fine powdered materials obtained were characterized by SEM,powder XRD,Raman spectroscopy,DSC-TGA etc.Multipoint Raman spectra showed consistent occurrence of spectral features indicating stoichiometric co-existence of ingredients in the crystal lattices.DSC analysis showed absence of all thermally assisted solidsolid phase transformation in the co-crystals as they were observed in pristine materials.The thermal stability calculated in terms of activation barrier for decomposition,revealed the CL-20/FOX-7 co-crystal to be intermediately stable on comparison to their constituents while,the HMX/FOX-7 co-crystal is more stable.Compared to pure HMX and CL-20,both the co-crystals have shown higher insensitivity to impact force,suggesting them to be suitable for future generation insensitive munitions.?2020 China Ordnance Society.Production and hosting by Elsevier B.V.on behalf of KeAi Communications Co.This is an open access article under the CCBY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    1.Introduction

    The most pronounced nitramine class of explosives such as hexahydro-1,3,5-trinitro-s-triazinane (RDX),tetrahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine(HMX)and 2,4,6,8,10.16-hexanitro-2,4,6,8,10,12-hexaazatetracyclo[5.5.0.0[3,11].0[5,9]]dodecane(CL-20)are characterized by velocity of detonation(VoD)well above 8.5 km/s.How ever,suffer from inherent limitations of their high sensitivity to mechanical stimuli.Explosive power and safety,often contradicting each other,are the two most concerned parameters in the field of high energy materials.The safest energetic materials(EMs) are least powerful such as trinitrotoluene (TNT),dinitrotoluene(DNT),trinitrotriaminobenzene(TATB)etc.,have VoD only near to 7.0 km/s[1].This has been a most apprehending circumstance w hile developing new energetic materials[2,3].Thus,there exists a grow ing demand for high performance and less sensitive(HPLS)energetic materials.As a matter of fact,only three energetic materials namely, nitrotriazolone (NTO), dinitrodiaminopyrazinoxide (LLM-105)and dinitrodiaminoethane(FOX-7)fall into this category with VoD above 8.5 km/s and impact insensitivity(h50%)more than 50 cm.Two pronounced approaches being adopted to reduce sensitivity of energetic materials are based on manipulating chemical stability[4-16](inclusion of-NH2,-OH,-COOH,conjugation,aromatic functionality etc.),and particle characteristics[17-19](reduced crystalline defects,polymorphism,spherocity,etc.).For existing energetic materials,how ever,chemical manipulation is limited and thus physical manipulation becomes the only option.A lot of efforts are being made to desensitize new age materials to improve their safety parameters however augmentation acquired in this direction is still inadequate.In recent days,co-crystallization has been widely accepted for solid modification in pharmaceuticals[20-22].The same concept is being experimented to manipulate hazard characteristics energetic materials.Co-crystals are built up owing to intermolecular nonbonding interactions such as hydrogen bonds,ionic bonds,π-π stacking and van der Waals forces.Presence of polar functional groups like-NO2,-NH2and a variety of acidic protons in nitramine explosive molecules greatly assist to form co-crystal.Progressively,new energetic co-crystals being produced opened up a potential route to high performance safer explosives and thus relaxing the contradiction[23-30].How ever,only a few of them such as CL-20/TNT(1:1);CL-20/HMX(2:1)and CL-20/BTF(1:1)have appreciable density and calculated Vo D above 8.5 km/s.Hence,further experiment to realize more powerful and better insensitive co-crystals of less sensitive EMs like FOX-7,TATBetc.,has become area of current interest.Among these,FOX-7,have most attractive performance indexes with density of 1.885 g/cm3,impact insensitivity upto h50%=126-159 cm,friction insensitivity>350 N and VoD of 9.1 km/s[31-34].Under this study,attempts were thus made to prepare co-crystals of CL-20 and HMX with FOX-7 as co-crystal former.Previously,Zhang et al.[29]proposed much attractive figures for densities of CL-20/FOX-7 and HMX/FOX-7 to be 1.993 g/cm3and 1.895 g/cm3,respectively.In a theoretical study published by Gao et al.[35]proposed high forming probability of 1:1 cocrystal of CL-20/FOX-7 system with a derived density of 2.001 g/cm3.As experimental studies are essential to practically realize such potential co-crystals of FOX-7,hardly any work is published in the recent literature.Using triamminotrinitrobenzene(TATB)as cocrystal former,a similar molecule with equivalent amino&nitro functionality like FOX-7,a 3:1 CL-20/TATB co-crystal is only reported till date[36].How ever,due to poor energetic performance of TATB,the overall activity of the co-crystal is not appreciable even though higher insensitivity could be expected.A previous study also reported HMX/TATBco-crystal with 8:1 stoichiometry[37]but,was later proved to be a physical mixture rather than co-crystal by morphological analysis[38].Attempt was thus made to prepare cocrystals of CL-20 and HMX with FOX-7 to impart a balance between insensitivity and performance.

    Initial attempts were made to obtain desired co-crystals through classical methods as followed in preparation of most pharmaceutical co-crystals.A brief account of screening of method is described in experimental section.All these thermodynamically favoured methods were unsuccessful to yield expected co-crystals.Hence,kinetic crystallization method,spray flash evaporation(SFE)was attempted to enable forced co-crystallization.Fine crystalline powders of desired co-crystals could be obtained successfully with this method.The powdered materials thus obtained were well characterized using powder X-ray diffraction,Raman spectrometry and SEM for chemical composition and particle morphology.The insensitivity parameters pertaining to thermal initiation and initiation by mechanical stimuli measured as,improved activation energy of decomposition and minimum height of impact were obtained using differential scanning calorimetry(DSC)and BAM Impact tester.

    2.Experimental

    2.1.Materials

    The energetic materials used for this study were obtained from various sources.ε-CL-20,HMX and FOX-7 were obtained from indigenous sources.All these materials were repeatedly recrystallized prior to their use to obtain chemical purity above 98.5%.All other chemical used were obtained from Thermo Fisher.UHP grade Nitrogen was used as pressurized gas for spray.

    2.2.Screening of co-crystallization method

    Screening out a suitable co-crystallization methodology has been the hard-hitting task under this study.Most effectively used techniques in different arena include solvent evaporation,grinding[39-46],slurry crystallization[47],melt extrusion[48],sonocrystallization[49,50]etc.Most of the energetic co-crystals realized so far have been prepared through solvent mediated process.High temperature melt extrusion how ever,not recommended for crystallization of thermally labile hazardous explosives.Wet grinding[51,52],a proven and safe method for preparation of fine energetic particulates was also expected to be useful.Thus,safer and thermodynamically favoured methods like solvent evaporation(fast and slow),w et grinding,ultrasound assisted cooling[53],vapour diffusion[54]etc.,have been attempted to co-crystallize FOX-7 with CL-20 and HMX.All these experiments were carried out using variety of solvents such as acetone,methanol,ethyl acetate,ethyl methyl ketone etc.Saturated solutions of explosive mixtures at room temperature were employed in each method.Thus,solvent volume was varied in all cases.However,each method yielded particulate mixtures rather than co-crystals as revealed by micro Raman analysis(Table 1).

    Relative solubility of the constituents in the solvent system plays major role during co-crystallization.FOX-7 exhibits strong intra and intermolecular hydrogen bonding with average distance of 2.31?[55]resulting in poor solubility in even polar solvents while,nitramines are appreciably soluble in the same medium.The homo intermolecular interaction in FOX-7 competed over to that in FOX-7/solvent and FOX-7/nitramine interactions resulting in separate homo molecular nucleation rather than forming hetero molecular clusters.In slow and fast evaporation,separated crystals were observed almost simultaneously.During vapour diffusion and ultra sound assisted cooling,in the first set of experiments(Table 1),FOX-7 partially crystallized out but CL-20 remained in solution. How ever, almost simultaneous but partial recrystallization of ingredients was observed for HMX/FOX-7 pair.These observations led to exclusion of thermodynamic methods for co-crystallization.

    Rapid expansion of supercritical solutions(RESS)was developed as an alternate method to prepare nano-explosives[56].This specific technique how ever,offer constrains in the form of poor solubility of explosives in supercritical CO2.Spray flash evaporation(SFE)or commonly know n as spray drying has been successfully implicated to generate pharmaceutical nano materials[57-59]and also nano energetic materials[60-62].Feasibility of this methodology has been demonstrated vide continuous preparation of nanoco-crystals of 1:1 CL-20/TNT and 2:1 CL-20/HMX[63,64].In recent study by Yang et al.[65],fine co-crystals of 1:1 BTF/TNT has been realized through spray drying.The study suggested spray drying as a fast,simple and continuous methodology for preparation and scale-up of energetic co-crystal.SFE facilitates extraordinarily fast crystallization compared to classical techniques through rapid solvent evaporation,which favours enhanced interaction between chemical entities at molecular level to produce crystalline hybrids.Hence,this method is expected to provide an opportunity to cocrystallize solutes with in-concomitant solubility.How ever,in flash evaporation,sometimes depending on the nature of intermolecular interactions between the chemical pairs,different crystalline composites like separated particulate mix and core-shells and semi crystalline composites may also form.

    Owing to advantages of rapid kinetic crystallization,SFE was attempted in this case and successfully yielded the desired cocrystals.The detailed experimental procedure is described in the section 2.3.

    Table 1 The observations from the classical methods of crystallization.

    2.3.Co-crystal preparation by SFE method

    The sub-micron size co-crystals were prepared by means of SFE that allow s a continuous preparative process for spherical submicron particulates of sensitive energetic materials within a single processing step.The basis of this process is to sudden drop in spray liquid pressure in a preheated large chamber.Solvents with boiling point below 60°C were used to suit the process.The substances under study were taken as solution in acetone in stoichiometric ratio.Unlike the commercial equipments,in the present set up,neither a heated nozzle nor did any pre heated drying gas was used.Instead,the chamber is externally heated to a temperature of 70°C.This makes the air inside the chamber hotter and partially creates a low pressure zone.The solution was then nebulized inside the chamber using an UHP nitrogen overpressure of 75 kg/cm2through a single fluid hollow cone nozzle of 0.3 mm diameter.Under a state of sudden liquid pressure drop in a heated atmosphere,the thermodynamic equilibrium is displaced,making the liquid droplets unstable.To regain stability,the excess thermal energy is converted into latent heat,enabling ultra fast evaporation of the solvent.Ultra fast removal of solvent forces the solute molecules to arrange themselves into most favorable lattice positions to built a crystalline material.Such a crystal would thus contain stoichiometric ratio of both the solutes within same lattice,thus forming a co-crystal.A schematic representation of the method is shown in Fig.1.Only those free flowing particles settled at the bottom of drying chamber were collected for characterization.

    Experiments were carried out in varied molar stoichiometry such as 6:1,4:1,2:1,1:1,1:2,1:4 and 1:6.The product from each experiment for both the explosive combinations i.e.CL-20/FOX-7 and HMX/FOX-7 were characterized by micro Raman spectrometer.Particle to particle probing helped to identify residual uncocrystallized materials.Multi point mapping has been extremely helpful to establish molar stoichiometry of complete cocrystallization.In this study,complete co-crystallization was obtained at 2:1 for CL-20/FOX-7 and 4:1 for HMX/FOX-7 combinations.Other stoichiometric combinations yielded mixture of cocrystal and pure crystals of the component that is in excess.Since,this study was aimed only to establish feasibility of spray technique to prepare co-crystals of non-concomitantly soluble energetic materials,quantitative calculation of yield has not been derived.

    2.4.Co-crystal characterization

    Fig.1.Experimental laboratory set-up and process scheme of Spray Flash Evaporation.

    Particle morphology were recorded using Quanta 200,FEI Netherlands make ESEM.All the imaging was performed under low vacuum and beam accelerating potential of 12.5-15 k V.Raman spectra were collected by using a‘in Via Reflex’dispersive micro Raman spectrometer from Renishaw,UK with 514 nm Argon ion laser.The spectrometer is attached with a‘Leica’DM 250 research grade optical microscope with motorized stage having low est translational step size of 1μm.Spectra were collected in the Raman shift region of 100-3500 cm-1at 1 cm-1spectral resolution under 10%of laser pow er input to prevent sample damage and minimize thermal noise.FTIR spectrums were collected using a Nicolet 5700 FTIR spectrometer in the spectral range 400-4000 cm-1at 1 cm-1spectral resolution and 32 scans per spectrum.In order to avoid grinding of the co-crystal particles,which could alter the lattice properties,FTIR spectrums were collected in ATR mode using a single bounce Diamond ATR accessory.Powder X-ray diffraction data of pure and co-crystal samples were collected on Bruker D8 Advance diffractometer with nickel filtered Cu Kα radiation(λ=1.5406?),run at 40 k V voltage and 40 n A current.Data were collected over diffraction angle range of 2θ=10°-70°at a scan rate of 0.1°/min.Thermal analysis were performed in a DSC 7,differential scanning calorimeter from Perkin Elmer,USA.Sample sizes were restricted below 0.5 mg for decomposition studies and>2 mg for phase transformation.Samples were taken in hermetically sealed Aluminum pans and were subjected to thermal profiling.A heating rate of 10°C/min was employed for solid-solid phase transition studies,while decomposition profiles were obtained at different heating rates viz.5°C/min,10°C/min,15°C/min and 20°C/min so as to derive kinetic parameters.

    2.5.Sensitivity analysis

    Insensitivity to impact force was measured by Bruceton Staircase method using a fall hammer weighing 2 kg.About 20 mg of the pure ingredients and their co-crystals were taken in stainless steel anvils and subjected to impact of the falling hammer from a certain height.Every impact test was noted for probable exposition.A positive result was marked to have occurred if there evidenced a loud explosion,ignition(flash and/or smoke),crackling,sparking and even brow n spots with burning smell from sample.In the event of no reaction,the test was repeated with a fresh sample from height 5 cm greater than previous and if a reaction occurred,fresh samples were tested at successively lower height.By using an up and down method and analyzing the data statistically using Bruceton staircase method,a height for 50%ignition probability designated as‘h50%’was determined.Ahigher value of h50%signifies better insensitivity.The experimental value of‘h50%’is obtained in‘cm’.How ever,for universal acceptability,the insensitivity to impact is expressed in Joules and is calculated as

    In the present measurement,fall hammer weight is 2 kg and value of g=9.8 m/s2.

    3.Results and discussion

    To demonstrate feasibility of to generate insensitive co-crystals of CL-20 and HMX with co-crystal former that differ largely in solubility,CL-20/FOX-7 and HMX/FOX-7 pairs were chosen.Selection of FOX-7 was made due to its energetic performance and insensitivity.FOX-7 with density 1.88 g/cm3and VoD 8870 m/sec[66]exhibited performance close to RDX,while,being significantly insensitive.Density has almost a linear co-relation to Vo D and to achieve VoD above 8500 m/sec,a density value of>1.78 g/cm3is expected for any energetic material[30].A co-crystal usually possesses density mutually compromised by the pure components but,is not always true[24,67].In cases,where significantly low dense co-crystal former is used,e.g.TNT in CL-20/TNT 1:1 co-crystal,both density and performance has been limited by the co-crystal former[68].FOX-7 as a co-crystal former thus would limit the extent of compromise in performance and impart superior insensitivity.

    The proposed co-crystals were characterized by SEM for morphology,Raman&FTIR for conformational analysis and DSC for thermal stability&solid-solid phase transitions.It is well conceived that,spectroscopic and thermal analysis may not authenticate cocrystals formation as these analyses do not explain the lattice structure how ever,powder X-ray diffraction can support to some extent.Thus,it was attempted to carry out single-crystal X-ray diffraction(SC-XRD)analysis of both the proposed co-crystals.It was how ever,not a successful attempt owing very poor or no diffraction intensity from the tiny particulates.Further,smaller size and agglomeration characteristics of the particulates posed significant challenge to mount a single particulate and their stability under prolonged X-ray exposure.SC-XRD measurements is know n to provide explicit evidence towards formation of a co-crystal,how ever,this technique may not be suitable for all co-crystals prepared through different techniques.In absence of SC-XRD analysis,it would be difficult for accurate characterization of the proposed co-crystals.Therefore,extreme care has been taken to differentiate the proposed co-crystals from their pure constituents while analyzing by the conventional tools.

    3.1.Solubility analysis

    Solubility difference of CL-20 and HMX with FOX-7 restricted our initial attempt to co-crystallize by solvent evaporation method.All the three ingredients had their maximum solubility in acetone with values 75 g/100 ml,1.85 g/100 ml and 0.247 g/100 ml,respectively.Both,slow and rapid evaporation were futile to produce cocrystals,rather,components crystallized out separately.This can be understood by considering two main factors-(i)solubility difference and(ii)strong intermolecular hydrogen bonding in FOX-7.The measured solubility ratios in acetone at 298 K for CL-20/FOX-7 and HMX/FOX-7 pairs are 300:1 and 8:1,respectively.Such large solubility differences pushes FOX-7 to attain super saturation preferentially and hence crystallized out priory.The crystal structure of FOX-7,has been reported to exhibit wave shaped lattice layers[69]with intra-layer hydrogen bonding and inter-layer weaker van der Waals interactions.The inter-layer distance(~3?)is shorter than observed in graphite(~3.5?)and TATB(~6?).Apart from that,two intra molecular hydrogen bonds exist within the molecule measuring bond distance 1.77?.Such strong intermolecular forces in solid state,essentially lowers solubility of FOX-7 in moderately polar solvents like acetone.The same reason also explains extremely insensitive nature of FOX-7 towards mechanical stimuli.

    3.2.Morphology and surface analysis

    SEM micrographs of spray dried materials,of CL-20/FOX-7(Figs.2(a)-(c))and HMX/FOX-7 (Figs.2(d)-(e))co-crystals,respectively,show ed spherical particulates with smooth surface texture.However,very few particulates show ed non geometric morphology other than the majority of spherical particles.This is understood as possible amorphous content.The mean particle size measured for CL-20/FOX-7(2:1)and HMX/FOX-7(4:1)co-crystal systems are Xmean=3.3 and 4.6μm,respectively,with narrow size distribution.The SEM images were analyzed using available image analysis software“Nmage J1.52a”developed by National Institute of Health,USA.Being a two dimensional image processing tool,it could not generate spherocity,rather obtained the values for circularity and roundness,which are equivalent parameters to spherocity.The results are summarized in Table 2.

    The results of image analysis suggest the CL-20/FOX-7 to be more spherical in shape than the HMX/FOX-7 particulates.However,the material collected from the wall of drying chamber show ed agglomeration with non-spherical particles as presented in Fig.2(g).This may be due to early particle adhesion facilitated by incomplete solvent evaporation before striking with the inner wall of drying chamber.This would thus be very important point to note for usefulness of this technique to generate spherical particles for improved the solid loading of composites drastically.

    Fig.2.SEM micrographs of(a)-(c)CL-20/FOX-7;(d)-(f)HMX/FOX-7 co-crystals collected from collection chamber and(g)HMX/FOX-7 particulates collected form drying chamber wall.

    Table 2 Summary of results from image analysis.

    In absence of SC-XRD analysis,careful analysis of crystal morphology by microscopic techniques may provide evidence of changes in unit cell geometries.For instance,by comparing the reported morphology of HMX/TATB co-crystal[37]and HMX morphology predicted using BFDH methodology,Wiscon et al.,could unambiguously prove the former to be re-crystallized HMX particles rather than co-crystals.How ever,the particles of proposed co-crystals generated under this study show much varied morphology than their pristine constituents.Fig.3 shows a schematic representation of morphological changes that occur during co-crystallization process.

    ε-CL-20 is obtained as bi-pyramidal crystals while,FOX-7 crystallizes as prismatic or cuboid shaped particles and often show stacking of multiple plates.How ever,the particulates obtained through spray method are mostly spherical in nature.Similarly,the particulates obtained from HMX/FOX-7 system are spherical in shape,largely varied form shape of pure HMX and FOX-7.This indeed serves as useful evidence towards lattice alteration and thus formation of co-crystals.

    3.3.Raman spectral analysis

    The Raman spectra of the energetic/energetic co-crystals are given in Figs.4(a)-(b).FOX-7 is identified by a strong-ONO bending vibration band at 856 cm-1,two moderately strong bands at 1205 and 1351 cm-1for C-N and-ONO symmetric stretching,and NH stretching vibrations at 3332&3425 cm-1.The spectra of co-crystal show a band near 281 cm-1,a feature observed in α and β phases of CL-20,differentiates the conformation from pristine ε phase.A careful look into spectral pattern suggests the CL-20 conformation in co-crystal to be close to β phase.The-ONO bending vibrations in β-CL-20 and β-HMX are noted at 832&833 cm-1,respectively,while,the O-N-O symmetric stretching were marked near 1324 and 1309 cm-1.In both CL-20 and HMX,multiple bands originating from C-H bending vibrations,which are absent in case of FOX-7,are observed nearly at 1250 cm-1.Another set of distinguishing vibrational band,the N-N stretching vibrations in CL-20 and HMX appear in the region of 950 cm-1.Other distinguishing features in CL-20/FOX-7 co-crystal are C-H bending vibrations near 1270 cm-1,a single strong band at 859 cm-1showing a slight shift,a doublet at 318 cm-1,ring deformation band at 281 cm-1,C-H stretching at 3044 cm-1and N-H stretching at 3332 &3425 cm-1.All these,Raman peaks leads to probable conclusion of co-crystal formation with CL-20 being in β-phase conformation.As the SFE set up had not been isolated from atmosphere,presence of atmospheric moisture is unavoidable.Being metastable,an isolated particulate of β-CL-20 is highly susceptible for conversion to α phase.Retention of the conformation during the process and later is indicative of sufficient stabilization of the phase induced by intermolecular interactions with FOX-7 molecules in the co-crystal lattice.Similarly,the HMX/FOX-7 co-crystal can be identified by distinguished features like,a broad band at 1348 cm-1originating from O-N-O symmetric stretching with a blue shift compared to FOX-7 and HMX,a very sharp and strongest band at 866 cm-1,a doublet at 464&490 cm-1,singlet at 631 cm-1etc.Moreover,appearance of C-H stretching vibrations at 2977&3063 cm-1and N-H stretching at 3337&3431 cm-1,with slight blue shift are indicative of co-crystal formation.

    In order to investigate homogeneity and therefore identifying isolated non-cocrystals,Raman mapping of both the proposed cocrystals were performed.An area of 60μm×60μm over the sample surface was mapped at 5μm step size thus generating 169 spectrums in totality.Otherwise,it can be understood that,Raman spectra were collected from 169 different locations confined within a 60μm×60μm area.The mappings were carried out using a 50×(NA 0.95)microscope eyepiece and 514 nm Ar ion laser source with approximate beam diameter of 330 nm at focal point,much low er than the observed particle size.This ensured collection of Raman spectrum even from smallest particles present in the sample.Collections were enabled with‘a(chǎn)uto focus track’to optimize the working distance matching to focal distance before every collection.The Raman maps of sample surface thus generated are presented in Fig.5.

    Fig.3.Morphological changes observed during co-crystallization of proposed co-crystals.

    Color gradient from red to dark red represents low to high peak intensity of the vibrational band under consideration and a white spot would represent no scattered intensity.In case of CL-20/FOX-7 2:1 co-crystal,the Raman mapping(Fig.5(a))was obtained for the peak at 281 cm-1,ring deformation of β-CL-20 conformer.The color gradient shows change in intensity of the band at different location on the mapped surface.This may be due to random orientation of the particles and uneven surface as weak Raman scattered lights are prone to be absorbed by surrounding particles as well.Overall,all the positions on mapped surface show ed presence of vibrational band at 281 cm-1,signifying homogeneous presence of co-crystal particles.A corresponding map(Fig.5(b))for HMX/FOX-7 cocrystal was obtained for the vibrational band at 631 cm-1.This also showed similar kind of color gradient and thus overall homogeneity of co-crystal particles.These observations strongly corroborate to the condition of co-crystal formation.How ever,the results could not produce any clear indication of isolated particles.

    3.4.Powder XRD analysis

    To con f i rm the formation of co-crystals,powder X-ray diffraction pattern in the range of 2θ=10°-50°was considered for all the co-crystals and the pristine components.In each case,a comparative analysis is presented in Figs.6(a)-(b).The XRD pattern of FOX-7 show very few diffraction intensities with distinguishing features at 2θ=14.9°,20.0°,26.8°,27.8°,31.0°etc.with the band at 2θ=27.8°be the most intense.While,β-CL-20 is characterized by numbers of diffraction peaks with most of them are of low intensity.The most prominent peaks are at 2θ=13.7°(most intense),18.0°,19.0°,24.1°,28.3°,30.3°etc.In contrast,the co-crystal shows XRD pattern with bands at 2θ=13.6°,15.2°,18.1°,20.2°,24.0°,28.2°,30.2°,32.6°and 36.6°etc.The diffraction band at 28.2°and 30.2°is comparatively broadened,something not very similar to pristine materials.Interestingly,the diffraction pattern of co-crystal does not show bands matching directly to that pristine FOX-7.How ever,a closure look reveals a few bands with significant broadening such as at 2θ=15.2°,20.2°and 30.2°nearly matches to FOX-7.This is indicative of inclusion of FOX-7 in co-crystal lattice rather than formation of any amorphous powder.The XRD profile of pristine β-HMX shows prominent features at 2θ=15.0°,16.0°,18.7°,20.9°,23.4°,27.7°,30.1°,32.3°,37.5°41.7°and 42.7°with band at 2θ=32.3°to be the most intense.Whereas,the HMX/FOX-7 co-crystal exhibits new distinguishing features at 2θ=14.6°,16.9°,20.4°,22.9°,24.5°,26.0°,27.9°,29.5°,31.8°,35.0°,37.0°etc.Moreover,the entire pattern also shows specific alterations.In case of FOX-7 and β-HMX,most intense bands noted at 2θ=27.8°and 32.3°respectively,whereas,in the co-crystal,three almost equally intense peaks are observed at 2θ=20.4°,27.9°,31.8°.The peaks at 20.4°,22.9°and 29.5°are some new features observed in the cocrystal.FOX-7,however,show a medium intense doublet in this region.Moreover,the FWHM of all the features in the co-crystal are comparatively wider than the diffraction bands of pristine components.

    Although,there may be several other factors that may influence peak broadening,formation of co-crystals with nano crystallite sizes cannot be ruled out.It is therefore likely that,the observed cocrystals particles may be of poly crystalline nature.Also,wide background is evident in the diffraction pattern of both the cocrystals,especially in case of HMX/FOX-7 co-crystal.It is therefore possible that,both the co-crystals may contain certain amount of amorphous content.The SEM images also show presence of particles with non-geometric shape,likely to be the amorphous phases.How ever,a quantitative analysis of amorphous content is beyond the scope of present study.

    3.5.FTIR spectroscopic analysis

    Fig.4.Comparative Raman spectra of CL-20/FOX-7 co-crystal(a),and HMX/FOX-7 co-crystal(b)in the range 100-1600 cm-1.

    As the information obtained from Raman spectrometry and powder X-ray diffraction analysis of HMX/FOX-7 co-crystal were not sufficient to establish the claim,FTIR spectral analysis were carried out in support.FTIR spectra of the pow der material obtained from spray drying of 4:1 stoichiometric mixture of HMX/FOX-7,as presented in Fig.7,show some features distinguishable from the pristine materials.The new material essentially show no hydrate peak in the region of 3500 cm-1and above which nullify any possibility of formation separate HMX phases.

    Furthermore,the new material contains vibrational bands at positions similar to that of β-HMX.The N-H stretching vibrations of FOX-7 at 3409 cm-1was observed to split into two distinguishable peaks at 3423 cm-1&3403 cm-1in the new material.Also,the strong N-H bending vibration of FOX-7 at 1623 cm-1has been manifested as two well separated peaks in the new material located at 1609 cm-1and 1637 cm-1.This indicates occurrence of two types of N-H stretching vibrations supposed to be originated from differently oriented intermolecular forces in the co-crystal.In a possible co-crystal lattice,the-NH hydrogen atoms of FOX-7 are believed to experience two types of intermolecular hydrogen bonds.One set of hydrogen bonds probably involved N-H&N-O of FOX-7 molecules(intermolecular)and other set involved N-H of FOX-7&N-O of neighboring HMX molecules in the same crystal lattice.The proposed intermolecular interactions are presentation in Scheme 1.

    Similarly,the C-H stretching vibration of HMX at 3031 cm-1split into two new peaks at 3036&3027 cm-1.A minor broadening of peak at 2985 cm-1(C-H stretching,HMX)was also noted in the new material.Further,noticeable changes are observed in the regions of 1400-1500 cm-1,a multiplet structure at 1241 cm-1,two weaker but new peaks at 1110&1169 cm-1,broadening of the peak at 1022 cm-1,a strong and completely new vibrational peak at 908 cm-1and splitting of vibrational peaks in the region of 700-770 cm-1.All these distinguishable features supports to the possibility of formation of a co-crystal lattice rather than an amorphous physical mixture or a composite.

    Fig.5.The Raman maps of sample surface.

    Fig.6.Comparative powder X-ray diffraction pattern of(a)2:1 CL-20/FOX-7 co-crystal and(b)4:1 HMX/FOX-7 co-crystal.

    Fig.7.Comparative FTIR spectrum of HMX/FOX co-crystal and pristine materials(a)3800-2800 cm-1,(b)1950-400 cm-1.

    3.6.Thermal analysis

    3.6.1.Thermally assisted phase transformation

    Fig.8.Comparative phase transition profiles collected at heating rate 10°C/min for(a)CL-20/FOX-7 co-crystal,and(b)HMX/FOX-7 co-crystal.

    Fig.9.Comparative decomposition profiles recorded at 10°C/min of(a)CL-20/:FOX-7 co-crystal,and(b)HMX/FOX-7 co-crystal.

    The thermally assisted endothermic solid-solid phase transformation and decomposition profiles of pristine explosives and their likely co-crystals are presented in Fig.8 and Fig.9,respectively.FOX-7 shows two endothermic transitions at 115°C and 168°C,corresponding to α→β and β→γ solid-solid phase transformations(Fig.8(a)).This is follow ed by two exothermic transitions at 232°C and 272°C,corresponding to two step thermal decomposition of FOX-7(Fig.9(a)).These values are well corroborating with literature reports[34,70].The solid-solid endothermic phase transformation in CL-20 is noted at 172°C representing the ε→γpolymorphic change.Further,in the context of observed CL-20 conformation in spray dried material,as described based on Raman spectroscopic studies,the temperature of thermally assisted β→γCL-20 phase transformation on freshly prepared β-CL-20 was also measured and found to be at 151°C.These values matched well with the literature reports[71,72].In case of the CL-20/FOX-7 2:1 spray dried material,however,no endothermic transition observed up to 200°C(Fig.8(a)).This clearly indicates absence of any isolated lattice structures for both FOX-7 and CL-20,conforming to the conclusion of co-crystal formation.CL-20 shows a wide decomposition exotherm with Tmax=255°C while,an entirely different decomposition pattern was observed for the co-crystal.A closer look into the pattern presented in Fig.9(a)reveals three exothermic peaks overlapping substantially with decomposition maxim as at 240°C,247°C and 253°C.The onset of first step of FOX-7 degradation might have been shifted towards a higher temperature closer to that of CL-20 decomposition.Strong intermolecular interactions,more precisely,the enhanced hydrogen bonding in the co-crystal is believed to be responsible for such alterations.Both inter and intra molecular hydrogen bonding contributes to improved thermal stability in many explosives like TATB,FOX-7 etc.The first maxima thus correspond to first step of FOX-7 decomposition;the second one is possibly due to decomposition of intermediates formed through high temperature reaction of evolved gases from FOX-7 and CL-20,and the last one corresponding to combined phenomena of CL-20 decomposition and 2nd step of FOX-7 degradation.It is important to note that,no separate exotherm corresponding to second step of FOX-7 decomposition was observed.This combined decomposition thus also contributes towards concluding the formation of co-crystal.

    β-HMX,show endothermic transition at 185°Crelated to β→δ solid phase change as presented in Fig.8(b).Eisenreich[73]reported this transition to occur at 193°C based on powder X-ray diffraction and Brill et al.reported the same at 185.5°C[74].In the HMX/FOX-7 co-crystal,how ever,no endothermic transitions are observed corresponding to those of FOX-7 and HMX solid-solid phase transformation.Moreover,the first step decomposition of FOX-7 is shifted to 247°C(Fig.9(b)),15°Chigher than that observed in pure FOX-7.The second decomposition is however merged with that of HMX decomposition with maximum temperature of the combined effect at 278°C.The initial decomposition products of FOX-7 instigated HMX to decompose at a lower temperature through secondary reaction in the crystal lattice.The reactive gases from thermolysis of FOX-7 gets into reaction with HMX before could leave the lattice.The thermal profiling is thus supporting our findings saying co-crystal formation in the lattice.

    3.6.2.Non-isothermal decomposition kinetics

    The kinetic parameters for decomposition of the CL-20/FOX-7 and HMX/FOX-7 co-crystal and the pure ingredients were derived from results of DSC analysis carried out at multiple heating rates following the Kissinger[75]and Ozawa[76,77]models.Fig.8 shows the comparative thermal decomposition profile of both the cocrystals along with their ingredients at heating rate of 10°C/min.

    The mathematical forms of Kissinger and Ozawa models are given as follow s-

    where,β-linear heating rate(K/min),R-gasconstant 8.314 J/(K·mol),A-pre exponential or frequency factor,E-activation energy,and T-peak decomposition temperature.

    These models are based on the fact that,with increasing heating rates during dynamic scanning,the temperature of maximum reaction shifts towards higher temperature[78].Many transitions such as evaporation,crystallization,curing reaction,decomposition,etc.are kinetic events and are function of both time and temperature.This implicates that such transitions(peak maxima)will shift to a higher temperature while heated under an elevatedrate.The fact underlying is that,a particular event would have less time at any specific temperature.Thus,the sensitivity lag to heating rate can be employed to identify kinetic transitions in the sample and thereby deriving their kinetic parameters.Under this study,the decomposition profiles were recorded at four different heating rates 5°C/min,10°C/min,15°C/min&20°C/min.The values of kinetic parameters,the apparent activation energies(E),linear corelation co-efficient(R2),pre exponential factor(A),rate constants(k)calculated using Kissinger and Ozawa methods are listed in Table 3.

    Table 3 Experimental Kinetic parameters for ε-CL-20,β-HMX and the Co-crystals.

    The kinetic parameters derived from DSC analysis of the CL-20/FOX-7 co-crystal suggest thermal stability of the co-crystals to be intermediate to that of CL-20 and FOX-7 with energy of activation 56.28 kcal/mol,almost 13 kcal higher than that of CL-20.FOX-7 exhibits a two step decomposition path.The first step decomposition data only has been considered to assess the stability of FOX-7 and its co-crystal and comparison thereof.The decomposition rate constant suggested a faster decomposition rate of CL-20/FOX-7 cocrystal over to that of CL-20.The second decomposition of the cocrystal with E=35.615 kcal/mol,how ever has been more favoured in comparison to both CL-20 and FOX-7.A self catalysis by the decomposition products of CL-20 is expected to be responsible for this low activation barrier.The‘E’of 1st decomposition in HMX/FOX-7 co-crystal has been significantly enhanced compared to that of pristine components with value of E=136.85 kcal/mol.This indicated improved thermal stability of the co-crystal.Such a higher activation barrier can be attributed to enhanced intermolecular interactions as compared to those in neat compounds.A significant rise in decomposition maxima by about 15°C at β=10°C/min is also supportive enhanced intermolecular interaction.A close packing lattice arrangement may have facilitated the heat dissipation within the particle bulk through a non contact energy transfer mechanism thus improving thermal stability.Although,there has been a rise in decomposition maxima,the HMX/FOX-7 co-crystal exhibited much faster decomposition rate.Enhanced surface area owing to spherical particle morphology and possible auto catalytic effects are suspected to be active in this case.From the derived results using Kissinger method,the Arrhenius equation of exothermic decomposition process of the cocrystals may be expressed as-

    3.7.Impact sensitivity analysis

    The sensitivity parameter in terms of fall weight Impact sensitivity was determined using a Fall Hammer apparatus.The measured impact insensitivity values for ε-CL-20,βHMX and FOX-7 were determined as 5.88 J,9.6 J and 18.23 J respectively.Whereas,the co-crystal were found suitably insensitive compared to CL-20 and HMX with impact insensitivity up to 13.13 J and 14.9 J,respectively,for CL-20/FOX-7(2:1)and HMX/FOX-7(4:1)cocrystals.These values indicate the co-crystals to be safer than the pure nitramine explosives.Inclusion of FOX-7 molecules in the cocrystal lattice through preferential hetero intermolecular hydrogen bodings attributed to the improvement of insensitive characteristics of the novel hybrids.

    The co-crystals thus generated through spray flash evaporation carry two very important characteristics(i)reduced sensitivity and(ii)spherical particle morphologies.Coexistence of both the qualities in the same material would make them significantly reliable for futuristic insensitive munitions with high explosive yield.

    4.Conclusion

    Feasibility of spray flash evaporation has been demonstrated for preparation of co-crystals of differently soluble explosive molecules.CL-20/FOX-7(2:1)and HMX/FOX-7(4:1)co-crystals were prepared and characterized for structure and thermal properties.Xray diffraction indicated polycrystalline nature of the particulates with crystallite sizes possibly in the nano domain.Raman and FTIR spectral analysis indicated possible co-crystal formation through alteration in vibrational pattern.Absence of thermally assisted phase transformations in the co-crystals,corresponding to that of pristine ingredients,demonstrated complete co-crystallization.In both the co-crystals,the first decomposition step of FOX-7 has been shifted towards higher temperature making the co-crystals thermally more stable.The sensitivity to impact figures also suggested the co-crystals to be safer than pure nitramine ingredients.Further exploration of this methodology would be helpful to produce cocrystals in continuous manner to support requisites of future Insensitive Munitions.

    Author contributions

    The manuscript was written through contributions of all authors.All authors have given approval to the final version of the manuscript.

    Notes

    The authors declare no competing financial interest.

    Acknowledgment

    The authors of this manuscript are highly obliged to the Director,HEMRL for his support to carry out the work.The authors also acknowledge Defense Research&Development Organization(DRDO),India for financial support towards this study.

    老女人水多毛片| 色综合亚洲欧美另类图片| 欧美+日韩+精品| 99热这里只有是精品50| 国产一级毛片在线| 亚洲av免费高清在线观看| 激情 狠狠 欧美| 97超视频在线观看视频| 中文欧美无线码| 高清午夜精品一区二区三区| 亚洲最大成人中文| 午夜免费男女啪啪视频观看| 成人美女网站在线观看视频| 三级国产精品片| 午夜福利视频1000在线观看| 亚洲欧美日韩卡通动漫| 欧美一区二区国产精品久久精品| 精品人妻视频免费看| 精品一区二区三区人妻视频| 色5月婷婷丁香| 丝袜喷水一区| 如何舔出高潮| 2021少妇久久久久久久久久久| 久久久精品欧美日韩精品| 国产精品人妻久久久影院| 男人舔女人下体高潮全视频| 日本五十路高清| av在线蜜桃| h日本视频在线播放| 中文字幕av在线有码专区| 国产综合懂色| 亚洲av不卡在线观看| 秋霞在线观看毛片| 日韩欧美精品免费久久| 日本免费a在线| 亚洲美女视频黄频| 国产真实伦视频高清在线观看| 在线观看av片永久免费下载| 亚洲av二区三区四区| 欧美日韩国产亚洲二区| 真实男女啪啪啪动态图| av线在线观看网站| 亚洲精品456在线播放app| 日日干狠狠操夜夜爽| 国产一区二区在线观看日韩| 国产视频内射| 久久综合国产亚洲精品| 国产成人freesex在线| 99热这里只有精品一区| 99热这里只有是精品在线观看| 久久精品夜色国产| 国产人妻一区二区三区在| 亚洲高清免费不卡视频| 日本爱情动作片www.在线观看| 国产高清国产精品国产三级 | 免费看av在线观看网站| 日本欧美国产在线视频| 丰满少妇做爰视频| av视频在线观看入口| 日本av手机在线免费观看| 九九热线精品视视频播放| 日韩三级伦理在线观看| 久久久成人免费电影| 亚洲av成人精品一二三区| 国产精品av视频在线免费观看| eeuss影院久久| 狂野欧美激情性xxxx在线观看| 成人国产麻豆网| 欧美日韩精品成人综合77777| 久久国产乱子免费精品| 久久久午夜欧美精品| 日韩精品青青久久久久久| 欧美一区二区精品小视频在线| 身体一侧抽搐| 中文天堂在线官网| 成人毛片a级毛片在线播放| 国产女主播在线喷水免费视频网站 | 七月丁香在线播放| 婷婷六月久久综合丁香| 国产成人免费观看mmmm| 亚洲乱码一区二区免费版| 日本色播在线视频| 久久国产乱子免费精品| 欧美日韩精品成人综合77777| av国产免费在线观看| 国产三级在线视频| 欧美日韩精品成人综合77777| 热99在线观看视频| 汤姆久久久久久久影院中文字幕 | 久久久久久久亚洲中文字幕| 亚洲精华国产精华液的使用体验| 国产在视频线在精品| 在线免费观看不下载黄p国产| 乱系列少妇在线播放| 国产成人aa在线观看| 非洲黑人性xxxx精品又粗又长| av在线播放精品| 白带黄色成豆腐渣| 最近中文字幕2019免费版| 中文亚洲av片在线观看爽| 国产高潮美女av| 成年女人看的毛片在线观看| 99视频精品全部免费 在线| 伦精品一区二区三区| 卡戴珊不雅视频在线播放| 国产精品久久久久久精品电影| 国产 一区精品| 亚洲乱码一区二区免费版| www日本黄色视频网| 日日啪夜夜撸| av播播在线观看一区| 少妇丰满av| 久久婷婷人人爽人人干人人爱| 国产一区二区在线观看日韩| 少妇人妻一区二区三区视频| 亚洲真实伦在线观看| av在线观看视频网站免费| 国产三级在线视频| www日本黄色视频网| 九九久久精品国产亚洲av麻豆| 精品国产三级普通话版| 99热网站在线观看| 免费电影在线观看免费观看| 亚洲一级一片aⅴ在线观看| 亚洲一级一片aⅴ在线观看| 亚洲国产欧美在线一区| 欧美一区二区国产精品久久精品| 天美传媒精品一区二区| 插阴视频在线观看视频| 好男人视频免费观看在线| 国产精品久久视频播放| 99久久精品国产国产毛片| 欧美丝袜亚洲另类| 在线免费观看不下载黄p国产| 国产黄色视频一区二区在线观看 | 久久精品国产99精品国产亚洲性色| 日本黄色片子视频| 国产淫片久久久久久久久| 大香蕉97超碰在线| 久久国内精品自在自线图片| 麻豆精品久久久久久蜜桃| 欧美成人午夜免费资源| 观看免费一级毛片| 麻豆一二三区av精品| 国产黄片美女视频| 国产免费福利视频在线观看| 91午夜精品亚洲一区二区三区| 亚洲乱码一区二区免费版| 国产不卡一卡二| 精品午夜福利在线看| 狂野欧美激情性xxxx在线观看| 天美传媒精品一区二区| 亚洲五月天丁香| av国产久精品久网站免费入址| 最近2019中文字幕mv第一页| 久久久国产成人免费| 麻豆av噜噜一区二区三区| 欧美另类亚洲清纯唯美| 久久久精品大字幕| 亚洲无线观看免费| 搞女人的毛片| 欧美最新免费一区二区三区| 97超碰精品成人国产| 国产精品1区2区在线观看.| 汤姆久久久久久久影院中文字幕 | 乱系列少妇在线播放| 免费黄色在线免费观看| 麻豆av噜噜一区二区三区| 国产乱来视频区| av在线老鸭窝| 中文资源天堂在线| 一个人看视频在线观看www免费| 九草在线视频观看| 国产午夜福利久久久久久| 亚洲av日韩在线播放| 午夜日本视频在线| 一区二区三区高清视频在线| 成人毛片a级毛片在线播放| 自拍偷自拍亚洲精品老妇| 22中文网久久字幕| 中文天堂在线官网| 久久久久精品久久久久真实原创| 亚洲av免费在线观看| 偷拍熟女少妇极品色| 丰满少妇做爰视频| 99在线视频只有这里精品首页| 成人二区视频| 久久久久久久亚洲中文字幕| 久久久久免费精品人妻一区二区| 国内精品一区二区在线观看| 熟女电影av网| 国产精品综合久久久久久久免费| 91午夜精品亚洲一区二区三区| 国产一区二区三区av在线| 亚洲精品色激情综合| 欧美性猛交╳xxx乱大交人| 亚洲av熟女| 日韩亚洲欧美综合| 韩国高清视频一区二区三区| 99热这里只有精品一区| 成人一区二区视频在线观看| 啦啦啦观看免费观看视频高清| 少妇的逼水好多| 欧美激情国产日韩精品一区| 九九在线视频观看精品| 美女大奶头视频| 亚洲国产精品久久男人天堂| 晚上一个人看的免费电影| 国产又黄又爽又无遮挡在线| 国产午夜精品一二区理论片| 人体艺术视频欧美日本| 久久欧美精品欧美久久欧美| 国产精品国产三级国产av玫瑰| 舔av片在线| 国产成人a∨麻豆精品| 久久精品国产亚洲网站| 日韩欧美国产在线观看| 少妇的逼水好多| 夫妻性生交免费视频一级片| av在线观看视频网站免费| 五月玫瑰六月丁香| 两个人视频免费观看高清| 男人和女人高潮做爰伦理| 亚洲国产成人一精品久久久| 国产午夜精品一二区理论片| 中文亚洲av片在线观看爽| 五月伊人婷婷丁香| 国产亚洲午夜精品一区二区久久 | 又黄又爽又刺激的免费视频.| www.av在线官网国产| 国产极品天堂在线| 婷婷六月久久综合丁香| 我要看日韩黄色一级片| 免费观看的影片在线观看| 国产精品久久电影中文字幕| 国产成人午夜福利电影在线观看| 国产乱人视频| 在线免费观看的www视频| 青春草亚洲视频在线观看| 日本一二三区视频观看| 亚洲国产精品久久男人天堂| 91狼人影院| 国产成人精品久久久久久| 九九久久精品国产亚洲av麻豆| 最近中文字幕高清免费大全6| 在线免费观看不下载黄p国产| 国产综合懂色| 亚洲av二区三区四区| 中文欧美无线码| 床上黄色一级片| 一本一本综合久久| 高清毛片免费看| 国产探花在线观看一区二区| 网址你懂的国产日韩在线| 亚洲精品久久久久久婷婷小说 | 国产一区二区在线av高清观看| 亚洲国产最新在线播放| 亚洲人成网站在线播| 高清视频免费观看一区二区 | 中文亚洲av片在线观看爽| 国产一区二区在线观看日韩| 听说在线观看完整版免费高清| 国产真实乱freesex| 亚洲国产精品久久男人天堂| 97超碰精品成人国产| 午夜福利高清视频| 亚洲av男天堂| 卡戴珊不雅视频在线播放| 男人狂女人下面高潮的视频| 18+在线观看网站| 日本-黄色视频高清免费观看| 久久婷婷人人爽人人干人人爱| 噜噜噜噜噜久久久久久91| 大香蕉97超碰在线| 精品人妻一区二区三区麻豆| 黄色配什么色好看| 2021天堂中文幕一二区在线观| 国产男人的电影天堂91| 精品久久久久久久末码| 午夜福利成人在线免费观看| 欧美+日韩+精品| 久久热精品热| 久久人人爽人人爽人人片va| 久久国内精品自在自线图片| 国产精品久久电影中文字幕| 国产伦在线观看视频一区| 国产精品综合久久久久久久免费| 亚洲美女视频黄频| 我要看日韩黄色一级片| 亚洲国产色片| 亚洲不卡免费看| 免费看光身美女| 国产精品无大码| 嫩草影院精品99| 男人和女人高潮做爰伦理| 99久久成人亚洲精品观看| 看免费成人av毛片| 国产午夜精品久久久久久一区二区三区| 日韩强制内射视频| 亚洲欧美一区二区三区国产| 大又大粗又爽又黄少妇毛片口| 亚洲电影在线观看av| 亚洲av男天堂| 一区二区三区免费毛片| av线在线观看网站| 少妇高潮的动态图| 日产精品乱码卡一卡2卡三| 亚州av有码| 久久国内精品自在自线图片| 老女人水多毛片| .国产精品久久| 亚洲av成人精品一二三区| 亚洲欧美日韩卡通动漫| 国产一级毛片在线| 麻豆一二三区av精品| 纵有疾风起免费观看全集完整版 | 国产午夜精品论理片| 1000部很黄的大片| 久久精品影院6| 久久久色成人| 别揉我奶头 嗯啊视频| 亚洲精品影视一区二区三区av| 国产黄色视频一区二区在线观看 | 少妇被粗大猛烈的视频| 又粗又硬又长又爽又黄的视频| 国产亚洲av嫩草精品影院| 久久精品熟女亚洲av麻豆精品 | 级片在线观看| 一级黄色大片毛片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 最后的刺客免费高清国语| 美女高潮的动态| 亚洲丝袜综合中文字幕| 国内精品美女久久久久久| 亚洲精品乱久久久久久| 禁无遮挡网站| 乱码一卡2卡4卡精品| 精品国产露脸久久av麻豆 | 欧美变态另类bdsm刘玥| 欧美日韩综合久久久久久| 国产精品一区二区三区四区久久| 午夜精品在线福利| 亚洲图色成人| 女人久久www免费人成看片 | 精品无人区乱码1区二区| kizo精华| 欧美日韩在线观看h| 高清日韩中文字幕在线| 久久这里只有精品中国| 搡女人真爽免费视频火全软件| 国产黄a三级三级三级人| 97超视频在线观看视频| 亚洲国产精品sss在线观看| videos熟女内射| 又粗又硬又长又爽又黄的视频| 亚洲av免费在线观看| 97热精品久久久久久| 天堂影院成人在线观看| 午夜福利在线观看免费完整高清在| 五月玫瑰六月丁香| 在线免费十八禁| 国产激情偷乱视频一区二区| 久久久午夜欧美精品| 国产精品女同一区二区软件| 欧美成人免费av一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 黄片无遮挡物在线观看| 超碰av人人做人人爽久久| 久热久热在线精品观看| 老司机影院毛片| 亚洲18禁久久av| 国产淫片久久久久久久久| 日日摸夜夜添夜夜添av毛片| 97人妻精品一区二区三区麻豆| 欧美成人一区二区免费高清观看| 国产亚洲av片在线观看秒播厂 | 免费一级毛片在线播放高清视频| 免费看a级黄色片| 久久6这里有精品| 久久久久久久久久久丰满| 欧美日韩一区二区视频在线观看视频在线 | av在线蜜桃| 亚洲国产欧洲综合997久久,| 国产真实乱freesex| 两性午夜刺激爽爽歪歪视频在线观看| www.色视频.com| 国产淫片久久久久久久久| 亚洲人成网站在线观看播放| av播播在线观看一区| 国产精品三级大全| 亚洲va在线va天堂va国产| 国产大屁股一区二区在线视频| 亚洲,欧美,日韩| 免费播放大片免费观看视频在线观看 | 深爱激情五月婷婷| 99久国产av精品| 中文亚洲av片在线观看爽| 麻豆国产97在线/欧美| 美女高潮的动态| 99久久九九国产精品国产免费| 麻豆国产97在线/欧美| 99久国产av精品| 欧美三级亚洲精品| 国国产精品蜜臀av免费| 免费看av在线观看网站| 欧美三级亚洲精品| 免费看美女性在线毛片视频| 亚洲图色成人| 成人毛片a级毛片在线播放| 视频中文字幕在线观看| 国产日韩欧美在线精品| 亚洲欧美精品专区久久| 日韩高清综合在线| 一级毛片aaaaaa免费看小| 91aial.com中文字幕在线观看| 亚洲欧美清纯卡通| 日韩在线高清观看一区二区三区| 少妇被粗大猛烈的视频| 一边摸一边抽搐一进一小说| 国产成人午夜福利电影在线观看| 国产精品一及| 欧美3d第一页| 国产免费福利视频在线观看| 69人妻影院| 国产 一区 欧美 日韩| 久久99蜜桃精品久久| 亚洲av免费高清在线观看| kizo精华| 蜜臀久久99精品久久宅男| 精品免费久久久久久久清纯| 国产成人a区在线观看| 亚洲精品乱码久久久久久按摩| 亚洲av电影在线观看一区二区三区 | 国产午夜精品论理片| 3wmmmm亚洲av在线观看| 免费观看的影片在线观看| 97超碰精品成人国产| 日韩中字成人| av国产免费在线观看| 亚洲精品国产av成人精品| 国产午夜精品一二区理论片| 久久99热这里只有精品18| 最近的中文字幕免费完整| 欧美日韩综合久久久久久| 欧美一区二区亚洲| 观看免费一级毛片| 蜜臀久久99精品久久宅男| 毛片一级片免费看久久久久| 亚洲成av人片在线播放无| 精品一区二区免费观看| 热99re8久久精品国产| 黄色日韩在线| 综合色av麻豆| 色尼玛亚洲综合影院| 日韩欧美 国产精品| 国产精品久久久久久av不卡| 欧美高清成人免费视频www| 亚州av有码| 国产精品蜜桃在线观看| 我的老师免费观看完整版| 91精品伊人久久大香线蕉| 国产成人午夜福利电影在线观看| 亚洲av免费在线观看| 国产91av在线免费观看| 1000部很黄的大片| 少妇猛男粗大的猛烈进出视频 | 最近视频中文字幕2019在线8| 国产亚洲av片在线观看秒播厂 | 国产又色又爽无遮挡免| 亚洲av成人精品一二三区| 色视频www国产| 小说图片视频综合网站| 在线观看美女被高潮喷水网站| 成人一区二区视频在线观看| 午夜a级毛片| 建设人人有责人人尽责人人享有的 | 精品人妻熟女av久视频| 搞女人的毛片| 淫秽高清视频在线观看| 精品国内亚洲2022精品成人| 欧美日韩精品成人综合77777| 亚洲人与动物交配视频| 国产熟女欧美一区二区| 天美传媒精品一区二区| 成人一区二区视频在线观看| 日日撸夜夜添| 中文字幕亚洲精品专区| 村上凉子中文字幕在线| 好男人在线观看高清免费视频| 老司机影院毛片| 国产高清有码在线观看视频| kizo精华| 长腿黑丝高跟| 成人亚洲精品av一区二区| 亚洲真实伦在线观看| 久久久国产成人免费| 亚洲中文字幕日韩| 一个人免费在线观看电影| av专区在线播放| 日日撸夜夜添| 老司机福利观看| 村上凉子中文字幕在线| 国产精品美女特级片免费视频播放器| 精品不卡国产一区二区三区| 亚洲av熟女| 国产伦一二天堂av在线观看| 淫秽高清视频在线观看| 国产伦精品一区二区三区四那| 国产成人免费观看mmmm| 如何舔出高潮| 日本与韩国留学比较| 狂野欧美白嫩少妇大欣赏| 边亲边吃奶的免费视频| 日韩亚洲欧美综合| 夫妻性生交免费视频一级片| 嫩草影院新地址| 欧美日本视频| www.色视频.com| 别揉我奶头 嗯啊视频| 亚洲国产精品sss在线观看| 99热这里只有是精品在线观看| 男的添女的下面高潮视频| 国产精品嫩草影院av在线观看| 国语自产精品视频在线第100页| 听说在线观看完整版免费高清| 99久久精品一区二区三区| 亚洲精品乱久久久久久| 高清毛片免费看| 亚洲精品日韩av片在线观看| 免费黄色在线免费观看| 日本免费一区二区三区高清不卡| 老司机影院成人| 大话2 男鬼变身卡| 最近手机中文字幕大全| 国产黄片美女视频| 久久6这里有精品| 日韩国内少妇激情av| 91久久精品国产一区二区三区| 国产精品野战在线观看| 中文精品一卡2卡3卡4更新| 亚洲av成人精品一区久久| av免费在线看不卡| 直男gayav资源| 美女大奶头视频| 国产精品国产三级国产专区5o | 2021少妇久久久久久久久久久| 国产av码专区亚洲av| 男人的好看免费观看在线视频| 亚洲国产高清在线一区二区三| 一区二区三区免费毛片| 亚洲av二区三区四区| av卡一久久| 免费av不卡在线播放| 久久韩国三级中文字幕| h日本视频在线播放| 嫩草影院入口| 黄色日韩在线| 一区二区三区免费毛片| 欧美区成人在线视频| videossex国产| 色综合色国产| 国产高清有码在线观看视频| 亚洲国产精品合色在线| 最近中文字幕高清免费大全6| 亚洲,欧美,日韩| av在线播放精品| 国内精品一区二区在线观看| 一二三四中文在线观看免费高清| 亚洲色图av天堂| 国产精品久久视频播放| 高清午夜精品一区二区三区| 99久国产av精品国产电影| 身体一侧抽搐| 日韩 亚洲 欧美在线| 汤姆久久久久久久影院中文字幕 | 国产精品精品国产色婷婷| 18+在线观看网站| 免费人成在线观看视频色| 淫秽高清视频在线观看| 91精品一卡2卡3卡4卡| 中国国产av一级| 伦精品一区二区三区| 我要搜黄色片| 欧美激情久久久久久爽电影| 嫩草影院入口| 成年女人看的毛片在线观看| 中文乱码字字幕精品一区二区三区 | 国产综合懂色| 亚洲国产欧洲综合997久久,| 99久久精品一区二区三区| 久久久精品欧美日韩精品| 麻豆成人av视频| 69人妻影院| 在线观看66精品国产| 久久久久九九精品影院| 国产精品久久久久久久电影| 午夜免费激情av| 少妇高潮的动态图| 亚洲欧美日韩卡通动漫| 日本熟妇午夜| 97热精品久久久久久| 国产精品久久久久久久电影| 精品欧美国产一区二区三| 边亲边吃奶的免费视频| 99热6这里只有精品| 91av网一区二区| 国产91av在线免费观看| 美女黄网站色视频| 久久这里有精品视频免费| 男女国产视频网站| 内射极品少妇av片p| 日韩av在线大香蕉| 美女脱内裤让男人舔精品视频| 男人和女人高潮做爰伦理| 成人毛片60女人毛片免费|