2D periodicity in lammps

thank you axel for reply.

  i will do it in a 3D system that is periodic in 2 directions and non
periodic in z direction.
   but i am having a problem with tip4p model.
    i already made a program for 2d ewald but i have to match my results
with lammps.
  when i am doing it in lammps , i am using following files,

  #FORCEFIELD.TIP4P
   mass 1 1.0008 # H
   mass 2 15.9994 # O

pair_style lj/cut/tip4p/long 2 1 1 1 0.15 14.9
pair_modify tail no
kspace_style ewald/tip4p 1.0e-5

pair_coeff * * 0.0 0.0
pair_coeff 2 2 0.0 3.1589

bond_style harmonic
bond_coeff 1 0.0 0.9572

angle_style harmonic
angle_coeff 1 0.0 109.47

#lmpfile
  #LAMMPS Input file generated by Avogadro
# tip4p

# Intialization
units real
dimension 2
boundary p p fm

atom_style full

# Atom Definition
read_data wat16.lmpdat
include forcefield.TIP4P
replicate 1 1 1
variable watMoleMass equal 18.0153 # /(g/mol)
variable nAvog equal 6.0221415e23 # Avogadro's number
variable watMoleculeMass equal (\{watMoleMass\}/{nAvog}) # /(g/molecule)
fix 1 all store/force
fix constraint all shake 0.0001 100 0 b 1 a 1

# Settings
#velocity all create 298.15 4928459 rot yes mom yes dist gaussian
fix ensemble all nve

timestep 1.0

# Output
dump 1 all custom 1 forces.atom fx fy fz
dump dumpXYZ all xyz 1 wat.trj.xyz
thermo_style multi
thermo 1

# Run the simulation
run 1

  but as i go through the lammps documentation. there is no command for
ewald method with tip4p model. so i have a question that can i do the 2d
simulation ( 3d system but 2d periodicity) with tip4p model in lammps

thank you

thank you axel for reply.

[…]

but as i go through the lammps documentation. there is no command for

ewald method with tip4p model. so i have a question that can i do the 2d
simulation ( 3d system but 2d periodicity) with tip4p model in lammps

what about kspace_style pppm/tip4p?

axel.