import three;

import x11colors;

real Rbond=206;

triple atom1=(0,0,0);

triple atom2=atom1+scale_len*Rbond*Z;

surface xe_atm=scale3(70*scale_atm)*unitsphere;

draw(shift(atom1)*xe_atm,lightblue);

draw(atom1--atom2, mediumgray+opacity(c_bdtransp)+linewidth(scale_bnd*4mm),currentlight);

label("$Xe1$",atom1 ,c_align);

draw(atom1--atom1+lmod1,red+size_shaft,Arrow3(size_arrow,size_arrow*1.2),currentlight);

import three;
import x11colors;

// initialization
currentprojection=perspective(200,350,100);
//currentprojection=orthographic(200,350,100);
//currentlight.background=palecyan;

real scale_atm=3/10;
real scale_bnd=0.7;
real scale_len=0.7;
real scale_nmd=200;

real c_bdtransp=0.4;
real size_arrow=12;
real size_shaft=3;
triple c_align=(0,8,8);

// Atom balls

surface xe_atm=scale3(70*scale_atm)*unitsphere;
surface f_atm=scale3(45*scale_atm)*unitsphere;

// Cartesian coordinates
real Rbond=206; // 2.06 Angstrom

triple atom1=(0,0,0);
triple atom2=atom1+scale_len*Rbond*Z;
triple atom3=atom1+scale_len*Rbond*X;
triple atom4=atom1+scale_len*Rbond*Y;
triple atom5=atom1-scale_len*Rbond*X;
triple atom6=atom1-scale_len*Rbond*Y;
triple atom7=atom1-scale_len*Rbond*Z;

// Normal modes (Oh/eg)

triple lmod1=( 0.000, 0.000, 0.000)*scale_nmd;
triple lmod2=( 0.000, 0.000, 0.577)*scale_nmd;
triple lmod3=(-0.289, 0.000, 0.000)*scale_nmd;
triple lmod4=( 0.000,-0.289, 0.000)*scale_nmd;
triple lmod5=( 0.289, 0.000, 0.000)*scale_nmd;
triple lmod6=( 0.000, 0.289, 0.000)*scale_nmd;
triple lmod7=( 0.000, 0.000,-0.577)*scale_nmd;

defaultrender=render(compression=Zero,merge=true);

// draw atoms
draw(shift(atom1)*xe_atm,lightblue);
draw(shift(atom2)*f_atm,palecyan);
draw(shift(atom3)*f_atm,palecyan);
draw(shift(atom4)*f_atm,palecyan);
draw(shift(atom5)*f_atm,palecyan);
draw(shift(atom6)*f_atm,palecyan);
draw(shift(atom7)*f_atm,palecyan);

// draw bonds
draw(atom1--atom2, mediumgray+opacity(c_bdtransp)+linewidth(scale_bnd*4mm),currentlight);
draw(atom1--atom3, mediumgray+opacity(c_bdtransp)+linewidth(scale_bnd*4mm),currentlight);
draw(atom1--atom4, mediumgray+opacity(c_bdtransp)+linewidth(scale_bnd*4mm),currentlight);
draw(atom1--atom5, mediumgray+opacity(c_bdtransp)+linewidth(scale_bnd*4mm),currentlight);
draw(atom1--atom6, mediumgray+opacity(c_bdtransp)+linewidth(scale_bnd*4mm),currentlight);
draw(atom1--atom7, mediumgray+opacity(c_bdtransp)+linewidth(scale_bnd*4mm),currentlight);

// labels
label("$Xe1$",atom1 ,c_align);
label("$F2 $",atom2 ,c_align);
label("$F3 $",atom3 ,c_align);
label("$F4 $",atom4 ,c_align);
label("$F5 $",atom5 ,c_align);
label("$F6 $",atom6 ,c_align);
label("$F7 $",atom7 ,c_align);

// draw arrows
draw(atom1--atom1+lmod1,red+size_shaft,Arrow3(size_arrow,size_arrow*1.2),currentlight);
draw(atom2--atom2+lmod2,red+size_shaft,Arrow3(size_arrow,size_arrow*1.2),currentlight);
draw(atom3--atom3+lmod3,red+size_shaft,Arrow3(size_arrow,size_arrow*1.2),currentlight);
draw(atom4--atom4+lmod4,red+size_shaft,Arrow3(size_arrow,size_arrow*1.2),currentlight);
draw(atom5--atom5+lmod5,red+size_shaft,Arrow3(size_arrow,size_arrow*1.2),currentlight);
draw(atom6--atom6+lmod6,red+size_shaft,Arrow3(size_arrow,size_arrow*1.2),currentlight);
draw(atom7--atom7+lmod7,red+size_shaft,Arrow3(size_arrow,size_arrow*1.2),currentlight);

8Â¥: Originally posted by neweroica at 2016-06-07 21:26:20
Ö§³Ö˜ÇÖ÷¡£ ÎÒÔںܾÃÒÔÇ°Ò²Ôø‡LÔ‡ÓÃAsymptoteºÍMetaPost®‹ß^·Ö×Ó½Y˜‹£¬Œ‘ÁË‚€½ÐNeatMolµÄºê°ü£ºhttp://bbs.ctex.org/forum.php?mod=viewthread&tid=46540