• <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.

    男女国产视频网站| 波多野结衣av一区二区av| av在线app专区| 国产又爽黄色视频| 好男人电影高清在线观看| 亚洲精品中文字幕在线视频| a级毛片黄视频| 咕卡用的链子| 99久久99久久久精品蜜桃| 午夜91福利影院| av电影中文网址| 精品卡一卡二卡四卡免费| 美女视频免费永久观看网站| 人人妻人人爽人人添夜夜欢视频| 免费高清在线观看日韩| 国产精品久久久久久精品古装| 亚洲精品自拍成人| 欧美成狂野欧美在线观看| av有码第一页| videos熟女内射| 91字幕亚洲| 婷婷色综合www| 男女国产视频网站| 日本一区二区免费在线视频| 性高湖久久久久久久久免费观看| 在线观看人妻少妇| 亚洲午夜精品一区,二区,三区| 欧美大码av| 99香蕉大伊视频| 久久精品国产亚洲av高清一级| 蜜桃在线观看..| 飞空精品影院首页| 欧美国产精品va在线观看不卡| 可以免费在线观看a视频的电影网站| 黄色一级大片看看| 97在线人人人人妻| 日韩 亚洲 欧美在线| 一级毛片女人18水好多 | 精品少妇久久久久久888优播| 亚洲国产av影院在线观看| 国产精品久久久久久人妻精品电影 | 大码成人一级视频| 久久影院123| 欧美国产精品一级二级三级| 男女国产视频网站| 精品国产国语对白av| 中文字幕精品免费在线观看视频| 欧美日韩av久久| 国产高清视频在线播放一区 | 亚洲视频免费观看视频| 波野结衣二区三区在线| 麻豆av在线久日| 国产欧美日韩综合在线一区二区| 精品久久蜜臀av无| 国产成人啪精品午夜网站| 校园人妻丝袜中文字幕| 中文字幕精品免费在线观看视频| 国产精品久久久久久精品古装| 午夜影院在线不卡| 亚洲男人天堂网一区| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩成人在线一区二区| 99热网站在线观看| 丰满人妻熟妇乱又伦精品不卡| 久久国产精品人妻蜜桃| 亚洲欧美一区二区三区久久| 在线观看www视频免费| 美女中出高潮动态图| 人体艺术视频欧美日本| 国产黄色免费在线视频| 一级毛片 在线播放| av在线播放精品| 2021少妇久久久久久久久久久| 国语对白做爰xxxⅹ性视频网站| 国精品久久久久久国模美| svipshipincom国产片| 久久久久精品国产欧美久久久 | 99国产精品99久久久久| 91成人精品电影| 叶爱在线成人免费视频播放| 久久久久国产一级毛片高清牌| 婷婷色综合www| 欧美精品高潮呻吟av久久| 精品国产一区二区三区久久久樱花| 极品人妻少妇av视频| 王馨瑶露胸无遮挡在线观看| 十八禁高潮呻吟视频| 色网站视频免费| 亚洲av男天堂| 欧美国产精品va在线观看不卡| 久久精品熟女亚洲av麻豆精品| 人成视频在线观看免费观看| 久久热在线av| av天堂在线播放| 99精国产麻豆久久婷婷| 日韩av免费高清视频| 午夜激情av网站| 国产亚洲欧美精品永久| 一级黄色大片毛片| 99国产精品免费福利视频| 波多野结衣一区麻豆| 国产1区2区3区精品| 在线亚洲精品国产二区图片欧美| 色网站视频免费| 精品久久蜜臀av无| 日韩中文字幕欧美一区二区 | 9191精品国产免费久久| 无遮挡黄片免费观看| 国产伦人伦偷精品视频| 国产亚洲欧美精品永久| 亚洲熟女毛片儿| 国产伦理片在线播放av一区| 精品国产一区二区久久| 老司机亚洲免费影院| avwww免费| 十八禁人妻一区二区| 青春草亚洲视频在线观看| 韩国高清视频一区二区三区| 岛国毛片在线播放| 国产一区亚洲一区在线观看| 亚洲中文字幕日韩| 亚洲人成77777在线视频| 日本vs欧美在线观看视频| 午夜福利影视在线免费观看| 久久久欧美国产精品| 免费高清在线观看视频在线观看| 国产成人a∨麻豆精品| 丝袜在线中文字幕| 久久精品国产亚洲av高清一级| 男女高潮啪啪啪动态图| 久久久久久免费高清国产稀缺| 国产在线视频一区二区| 精品国产一区二区久久| 久久久久国产精品人妻一区二区| 亚洲久久久国产精品| 亚洲av成人精品一二三区| 99九九在线精品视频| 搡老岳熟女国产| 国产精品久久久久久人妻精品电影 | 男的添女的下面高潮视频| 男女高潮啪啪啪动态图| 人人妻,人人澡人人爽秒播 | 国产精品久久久久成人av| 国产免费又黄又爽又色| 日本av免费视频播放| 国产成人精品久久二区二区免费| 久久亚洲精品不卡| 制服诱惑二区| 亚洲国产欧美网| 美女中出高潮动态图| 久久精品亚洲av国产电影网| 丰满少妇做爰视频| h视频一区二区三区| 汤姆久久久久久久影院中文字幕| 色播在线永久视频| 午夜影院在线不卡| 在线观看免费日韩欧美大片| 一本综合久久免费| 欧美黄色片欧美黄色片| 各种免费的搞黄视频| 97人妻天天添夜夜摸| 午夜福利,免费看| 日本wwww免费看| 女性被躁到高潮视频| 亚洲视频免费观看视频| 18禁裸乳无遮挡动漫免费视频| 97精品久久久久久久久久精品| 成年人免费黄色播放视频| 极品人妻少妇av视频| 9热在线视频观看99| 免费在线观看完整版高清| 国产国语露脸激情在线看| 婷婷色综合www| 久久天堂一区二区三区四区| 日本wwww免费看| 亚洲欧美中文字幕日韩二区| 亚洲中文av在线| 亚洲国产欧美在线一区| 亚洲成人免费电影在线观看 | a级毛片黄视频| 国产成人一区二区在线| e午夜精品久久久久久久| 精品福利观看| 国产精品一区二区在线观看99| 中文精品一卡2卡3卡4更新| 麻豆乱淫一区二区| 亚洲成色77777| √禁漫天堂资源中文www| 欧美日韩黄片免| 亚洲av电影在线进入| av网站在线播放免费| 丰满饥渴人妻一区二区三| 色网站视频免费| av一本久久久久| 欧美+亚洲+日韩+国产| 满18在线观看网站| 免费一级毛片在线播放高清视频 | 丝瓜视频免费看黄片| 纵有疾风起免费观看全集完整版| 91精品国产国语对白视频| 巨乳人妻的诱惑在线观看| 亚洲美女黄色视频免费看| 欧美日韩国产mv在线观看视频| 国产av一区二区精品久久| 精品熟女少妇八av免费久了| 久久久久国产精品人妻一区二区| 性高湖久久久久久久久免费观看| 国产男女内射视频| 乱人伦中国视频| 日韩人妻精品一区2区三区| 亚洲 欧美一区二区三区| 国产一区有黄有色的免费视频| 亚洲人成电影观看| 热99国产精品久久久久久7| 亚洲男人天堂网一区| 首页视频小说图片口味搜索 | 一区二区av电影网| 亚洲精品一卡2卡三卡4卡5卡 | 最近最新中文字幕大全免费视频 | 亚洲欧美一区二区三区黑人| 国产成人免费无遮挡视频| 欧美97在线视频| 欧美日韩一级在线毛片| 久久国产精品影院| 天堂中文最新版在线下载| 亚洲精品久久午夜乱码| 两个人免费观看高清视频| 在线观看免费午夜福利视频| a级毛片黄视频| 五月开心婷婷网| www.精华液| av线在线观看网站| 国产一区二区三区av在线| 日韩免费高清中文字幕av| 欧美老熟妇乱子伦牲交| 啦啦啦在线观看免费高清www| 999精品在线视频| 欧美xxⅹ黑人| 亚洲欧美一区二区三区久久| 一二三四社区在线视频社区8| 日韩中文字幕视频在线看片| av在线老鸭窝| 麻豆国产av国片精品| 在现免费观看毛片| 亚洲人成77777在线视频| 久久热在线av| 午夜福利一区二区在线看| 99久久精品国产亚洲精品| 在线观看www视频免费| 久久中文字幕一级| 啦啦啦视频在线资源免费观看| 亚洲国产看品久久| 99热网站在线观看| 久久人人97超碰香蕉20202| 一级,二级,三级黄色视频| 精品久久久精品久久久| 精品少妇黑人巨大在线播放| 男女国产视频网站| 一区福利在线观看| 久久久久视频综合| 91成人精品电影| 69精品国产乱码久久久| 一区二区三区激情视频| 亚洲欧美一区二区三区黑人| 女人爽到高潮嗷嗷叫在线视频| 久久av网站| 国产精品人妻久久久影院| 午夜视频精品福利| 超色免费av| 国产精品久久久久成人av| av有码第一页| 亚洲色图综合在线观看| 两个人免费观看高清视频| 一级毛片 在线播放| 免费在线观看影片大全网站 | 欧美 亚洲 国产 日韩一| 日本黄色日本黄色录像| 久久影院123| 精品国产一区二区三区四区第35| 中文字幕色久视频| 97在线人人人人妻| 婷婷色综合大香蕉| 久久精品久久久久久久性| 在线看a的网站| 9191精品国产免费久久| 日韩熟女老妇一区二区性免费视频| 男女国产视频网站| 成人午夜精彩视频在线观看| 亚洲综合色网址| 97在线人人人人妻| 精品久久久久久电影网| 狠狠精品人妻久久久久久综合| 麻豆av在线久日| 欧美性长视频在线观看| 亚洲九九香蕉| 大片免费播放器 马上看| 久久av网站| 国产精品 欧美亚洲| 涩涩av久久男人的天堂| 99热网站在线观看| 亚洲少妇的诱惑av| 少妇人妻 视频| 婷婷成人精品国产| 午夜av观看不卡| 97人妻天天添夜夜摸| 免费女性裸体啪啪无遮挡网站| 亚洲精品av麻豆狂野| 国产一卡二卡三卡精品| 国产精品一区二区在线不卡| 国产av一区二区精品久久| 91麻豆av在线| 亚洲国产毛片av蜜桃av| 99热网站在线观看| 国产亚洲精品第一综合不卡| 女人精品久久久久毛片| videos熟女内射| 视频区图区小说| 免费日韩欧美在线观看| 美女中出高潮动态图| 欧美国产精品va在线观看不卡| 丝袜喷水一区| 久久人妻熟女aⅴ| 最近中文字幕2019免费版| 国产精品人妻久久久影院| 一边摸一边抽搐一进一出视频| 欧美在线黄色| 性色av乱码一区二区三区2| 色播在线永久视频| 欧美日韩综合久久久久久| 精品高清国产在线一区| 女人久久www免费人成看片| 免费在线观看日本一区| 黄色视频不卡| 真人做人爱边吃奶动态| 国产成人精品在线电影| 国产成人a∨麻豆精品| 黄色视频不卡| 91精品伊人久久大香线蕉| 91精品三级在线观看| 中文字幕最新亚洲高清| 免费在线观看完整版高清| 欧美精品av麻豆av| 亚洲视频免费观看视频| 欧美久久黑人一区二区| 人人妻人人爽人人添夜夜欢视频| 久久久久久久精品精品| 欧美激情极品国产一区二区三区| 午夜福利视频精品| 久久精品人人爽人人爽视色| 午夜91福利影院| 人人妻人人澡人人爽人人夜夜| 啦啦啦视频在线资源免费观看| 在线天堂中文资源库| 国产精品一区二区免费欧美 | 看免费成人av毛片| 亚洲 欧美一区二区三区| 久久久久精品人妻al黑| 亚洲精品久久午夜乱码| 五月开心婷婷网| 免费一级毛片在线播放高清视频 | 尾随美女入室| 久久综合国产亚洲精品| 91精品三级在线观看| 久久久久久久精品精品| 女性生殖器流出的白浆| 午夜91福利影院| 成在线人永久免费视频| 一区二区三区乱码不卡18| 叶爱在线成人免费视频播放| 热99国产精品久久久久久7| 亚洲欧美一区二区三区黑人| 午夜激情av网站| a级毛片在线看网站| 亚洲专区国产一区二区| 成人手机av| 一区二区三区激情视频| 99香蕉大伊视频| 狠狠精品人妻久久久久久综合| 国产极品粉嫩免费观看在线| 99久久人妻综合| 国产成人免费观看mmmm| 欧美av亚洲av综合av国产av| 日本猛色少妇xxxxx猛交久久| 久久亚洲精品不卡| 亚洲国产精品一区三区| 中文字幕人妻丝袜一区二区| 亚洲 欧美一区二区三区| 国产精品九九99| 亚洲精品av麻豆狂野| 看十八女毛片水多多多| 男男h啪啪无遮挡| 国产免费现黄频在线看| 国语对白做爰xxxⅹ性视频网站| 你懂的网址亚洲精品在线观看| 超碰成人久久| 男女边摸边吃奶| 久久国产精品影院| 国产一区二区三区av在线| 女人被躁到高潮嗷嗷叫费观| 日本vs欧美在线观看视频| 久久久国产欧美日韩av| 国产男女内射视频| 欧美激情 高清一区二区三区| 人人澡人人妻人| 99re6热这里在线精品视频| 色视频在线一区二区三区| 在线天堂中文资源库| 国产成人精品久久二区二区免费| 性色av乱码一区二区三区2| 亚洲一区中文字幕在线| 满18在线观看网站| 精品少妇一区二区三区视频日本电影| 国产色视频综合| 久久久久久久大尺度免费视频| 99香蕉大伊视频| 黑人巨大精品欧美一区二区蜜桃| 一级毛片 在线播放| 欧美在线黄色| 制服人妻中文乱码| 捣出白浆h1v1| 色播在线永久视频| 国产成人精品无人区| 97人妻天天添夜夜摸| 免费高清在线观看视频在线观看| 国产欧美日韩一区二区三区在线| 在线av久久热| 日韩 欧美 亚洲 中文字幕| 色综合欧美亚洲国产小说| 80岁老熟妇乱子伦牲交| 久久人人爽av亚洲精品天堂| 色94色欧美一区二区| 午夜91福利影院| 男女免费视频国产| 久久国产精品影院| 中文字幕人妻熟女乱码| 午夜福利视频在线观看免费| 大片电影免费在线观看免费| 国产不卡av网站在线观看| 久久久久精品国产欧美久久久 | 看十八女毛片水多多多| 成人影院久久| 婷婷色av中文字幕| 亚洲中文字幕日韩| 午夜福利免费观看在线| 99热网站在线观看| 91字幕亚洲| 视频区图区小说| www.自偷自拍.com| 国产成人影院久久av| 天天躁日日躁夜夜躁夜夜| 啦啦啦 在线观看视频| 丰满人妻熟妇乱又伦精品不卡| 天天影视国产精品| 一本—道久久a久久精品蜜桃钙片| 欧美成人午夜精品| 欧美精品一区二区大全| 91九色精品人成在线观看| 国产三级黄色录像| 国产高清国产精品国产三级| 亚洲第一av免费看| 别揉我奶头~嗯~啊~动态视频 | 亚洲精品av麻豆狂野| 国产91精品成人一区二区三区 | 亚洲综合色网址| 国产精品久久久av美女十八| 国产高清国产精品国产三级| 亚洲中文av在线| 又黄又粗又硬又大视频| 又大又爽又粗| 纯流量卡能插随身wifi吗| 亚洲欧美成人综合另类久久久| 精品亚洲成a人片在线观看| 日本猛色少妇xxxxx猛交久久| 人人妻人人爽人人添夜夜欢视频| 精品国产国语对白av| cao死你这个sao货| 性色av一级| 亚洲国产av影院在线观看| 一个人免费看片子| 一级片'在线观看视频| 欧美黄色片欧美黄色片| 精品少妇内射三级| 久久99一区二区三区| 一边摸一边做爽爽视频免费| 美女主播在线视频| 亚洲第一av免费看| 天堂8中文在线网| 久久精品人人爽人人爽视色| a 毛片基地| 日本色播在线视频| 精品亚洲成a人片在线观看| 欧美少妇被猛烈插入视频| 九草在线视频观看| 免费在线观看视频国产中文字幕亚洲 | 大香蕉久久成人网| 女人高潮潮喷娇喘18禁视频| √禁漫天堂资源中文www| 久久鲁丝午夜福利片| 人妻人人澡人人爽人人| 国产精品国产三级专区第一集| 久久人人97超碰香蕉20202| 久久久国产一区二区| 亚洲久久久国产精品| 天天躁日日躁夜夜躁夜夜| 日韩av在线免费看完整版不卡| 悠悠久久av| 十分钟在线观看高清视频www| 热99久久久久精品小说推荐| 久久 成人 亚洲| 国产伦人伦偷精品视频| 热re99久久精品国产66热6| 七月丁香在线播放| 亚洲精品久久成人aⅴ小说| 可以免费在线观看a视频的电影网站| 午夜福利视频在线观看免费| 亚洲成人免费av在线播放| 男人操女人黄网站| 欧美成人午夜精品| 一级,二级,三级黄色视频| 曰老女人黄片| 精品福利永久在线观看| 高清欧美精品videossex| 国产精品一区二区在线观看99| 99久久综合免费| 大陆偷拍与自拍| 9191精品国产免费久久| 久久中文字幕一级| 又黄又粗又硬又大视频| 天堂中文最新版在线下载| 女性被躁到高潮视频| 99国产精品一区二区蜜桃av | 欧美亚洲日本最大视频资源| 日本午夜av视频| 少妇粗大呻吟视频| 日本午夜av视频| 国产爽快片一区二区三区| 日本vs欧美在线观看视频| 国产野战对白在线观看| 可以免费在线观看a视频的电影网站| 亚洲男人天堂网一区| 无限看片的www在线观看| 最新的欧美精品一区二区| 韩国精品一区二区三区| 欧美精品av麻豆av| 丝袜美腿诱惑在线| 国产精品国产三级国产专区5o| 国产成人免费观看mmmm| 一本一本久久a久久精品综合妖精| 一本大道久久a久久精品| 香蕉国产在线看| 无遮挡黄片免费观看| 50天的宝宝边吃奶边哭怎么回事| 精品国产一区二区三区久久久樱花| 午夜两性在线视频| 各种免费的搞黄视频| 十八禁网站网址无遮挡| 亚洲午夜精品一区,二区,三区| 日本一区二区免费在线视频| 国产av一区二区精品久久| 国产免费一区二区三区四区乱码| 欧美亚洲 丝袜 人妻 在线| 免费黄频网站在线观看国产| 亚洲欧美日韩高清在线视频 | 欧美日韩综合久久久久久| 午夜两性在线视频| 久久久久视频综合| 成年人免费黄色播放视频| 黄色视频在线播放观看不卡| 又大又黄又爽视频免费| 国产av一区二区精品久久| 欧美日韩国产mv在线观看视频| 在线观看免费日韩欧美大片| 日本欧美视频一区| 啦啦啦视频在线资源免费观看| a级毛片在线看网站| 十八禁高潮呻吟视频| 人人澡人人妻人| 性色av乱码一区二区三区2| 晚上一个人看的免费电影| √禁漫天堂资源中文www| 午夜福利在线免费观看网站| 久久久精品区二区三区| 色综合欧美亚洲国产小说| 国产欧美亚洲国产| 男人爽女人下面视频在线观看| 亚洲欧美成人综合另类久久久| 1024香蕉在线观看| 日本欧美国产在线视频| 视频区图区小说| 亚洲精品一二三| 亚洲精品国产av成人精品| 老鸭窝网址在线观看| 亚洲中文字幕日韩| av在线老鸭窝| 最新在线观看一区二区三区 | 欧美激情 高清一区二区三区| 9色porny在线观看| 日韩熟女老妇一区二区性免费视频| 国产在视频线精品| 亚洲av欧美aⅴ国产| 热re99久久精品国产66热6| 一区二区av电影网| 一个人免费看片子| 丰满人妻熟妇乱又伦精品不卡| 免费看不卡的av| 日日爽夜夜爽网站| 成在线人永久免费视频| 欧美精品啪啪一区二区三区 | 波野结衣二区三区在线| 亚洲精品日本国产第一区| 日韩人妻精品一区2区三区| 午夜两性在线视频| 观看av在线不卡| 国产成人系列免费观看| 国产视频首页在线观看|